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ABSTRACT 

All  Navy  aircraft  are  required  to  accommodate, 
anthropometrically ,  ninety  percent  of  the  user  population. 
Some  designs  have  been  criticized  for  their  low  accommo- 
dations but  those  accommodations  have  never  been  quantified. 
The  purpose  of  this  thesis  was  to  quantify  the  accommoda- 
tion, by  each  type  of  operational  Naval  aircraft,  of 
populations  of  Naval  aviation  personnel  of  1964,  1969,  and 
1975.   The  Computerized  Accommodation  Percentage  Evaluation 
(CAPE)  model  was  used  to  generate  data  points  since  only 
summary  statistics  were  available  for  two  of  the  popula- 
tions.  Each  subject  of  every  population  was  checked  against 
the  requirements  of  the  design  specification,  and  against 
the  limitations  of  each  aircraft.   All  aircraft  were  found 
to  accommodate  more  than  ninety  percent  of  the  1975  popu- 
lation.  Time  related  changes  in  the  populations  were  noted 
and  unexplained  inconsistencies  in  the  data  were  discovered. 
Possible  sources  of  error  were  discussed  and  potential 
solutions  proposed. 
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I.   INTRODUCTION 

The  Human  Engineering  Design  Criteria  for  Military 

Systems,  Equipment  and  Facilities  [MIL-STD-1472B]  states 

the  requirement  for  anthropometric  accommodation  of 

equipment  users: 

Design  and  sizing  shall  insure  accommodation,  com- 
patibility, and  maintainability  by  at  least  90 
percent  of  the  user  population.   Generally,  design 
limits  shall  be  based  upon  a  range  from  the  5th 
percentile  to  the  95th  percentile  values  for 
critical  body  dimensions. 

For  the  special  case  of  Naval  Aircraft  cockpits,  the 
accommodation  criterion  was  enlarged  in  1973  to  include  the 
central  95  percent  of  the  user  population  CNAVAIR  SD-24K, 
1973] .   The  anthropometric  description  of  the  population  is 
based  on  data  collected  in  1964.   Although  a  future  survey, 
scheduled  for  1981,  will  provide  a  more  current  description 
of  the  anthropometric  features  of  Naval  aviation  personnel, 
aircraft  now  being  designed  for  production  and  use  well  into 
the  future  are  being  designed  to  accommodate  the  1964  popu- 
lation.  The  F-18,  projected  to  become  operational  in  1982 
will  be  built  to  accommodate  the  pilots  of  1964.   Other 
designs,  if  they  have  advanced  much  beyond  the  concept 
stage,  are  committed  to  the  same  aircrew  accommodation 
standards. 
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The  anthropometric  characteristics  of  the  population 
are  changing  but  the  design  specification  is  not.   That  the 
population  is  changing  is  reflected  in  efforts  by  the 
National  Aeronautics  and  Space  Administration  and  by  the 
United  States  Air  Force  to  predict  the  typical  pilot  size 
of  the  future.   A  1978  study  conducted  by  NASA  concluded 
that  stature  (height)  increases  among  pilot  and  potential 
astronaut  males  could  be  expected  to  amount  to  about  8  mm 
(1/3  inch)  per  decade  [NASA  1024  Vol.  I,  1978].   The 
U.S.  Air  Force  predicted  an  increase  in  height  of  0.418 
inch  per  decade   Roebuck,  Kroemer,  and  Thompson,  1975  . 
These  changes,  since  1964  and  into  the  future,  affect 
pilot  accommodation  in  operational  Naval  aircraft,  but  the 
exact  nature  of  the  effect  is  not  known. 

In  1977  the  Navy  implemented  the  Anthropometric 
Compatibility  Assignment  Program  [OPNAVINST  3710.36],  a 
system  aimed  at  preventing  the  assignment  of  Naval  Aviators 
and  Naval  Flight  Officers  to  aircraft  which  are  incompatible 
with  their  anthropometric  features.   Under  the  provisions  of 
this  program,  pilots'  and  NFOs '  anthropometric  measurements 
are  taken  when  they  enter  flight  training  CBUMEDINST  3710.1, 
1977].   Values  for  sitting  height,  functional  reach, 
buttock-knee  length,  and  leg  length  (for  descriptions  of 
anthropometric  terms,  see  Appendix  A)  are  encoded  and 
become  part  of  their  permanent  records.   For  example,  a 
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man  coded 

6  4  2  3 
would  have  a  sitting  height  measurement  between  38.5  and 
38.9  inches,  a  functional  reach  of  30.5  to  30.9  inches, 
buttock-knee  length  of  23.0  to  23.9  inches  and  a  leg  length 
of  40.0  to  42.9  inches. 

The  cockpits  of  each  operational  aircraft  type  are 
measured  to  determine  the  limits  of  pilot  size  compatible 
with  the  cockpit  [NAVAIRINST  3710.9,  1979].   The  maximum 
and  minimum  acceptable  dimensions  for  sitting  height, 
functional  reach,  buttock-knee  length,  and  leg  length  are 
determined  for  each  cockpit.   These  maxima  and  minima  are 
then  converted  to  four  coded  descriptors,  which  signify 
dimensions  of  exclusion  for  pilots.   For  example,  an  A-4E 
is  coded: 

089  9  789  0 
As  with  the  personal  anthropometric  codes,  the  first  of  the 
four  groups  of  numbers  corresponds  to  sitting  height.   The 
second  group  corresponds  to  functional  reach;  the  third  to 
buttock-knee  length,  and  the  fourth  to  leg  length.   The 
code  089  9  789  0  indicates  that  a  pilot  with  a  sitting 
height  code  0,  8,  or  9  would  not  be  eligible  for  assignment 
to  an  A-4E.   Similarly  any  pilot  with  a  functional  reach 
code  9,  buttock-knee  length  code  7,  8,  or  9,  or  a  leg 
length  code  0  would  be  excluded. 

The  Chief  of  Naval  Aviation  Training  (CNATRA)  is 
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required  to  utilize  these  codes  early  in  each  aviation 
student's  training  to  exclude  the  student  from  training 
sequences  leading  to  aircraft  with  which  he  is  not  com- 
patible [CNATRAINST  13520.1,  1980} .   The  Navy  Military 
Personnel  Command  uses  the  anthropometric  codes  in  the 
process  of  assigning  pilots  and  NFOs  to  aircraft  squadrons. 

Currently,  a  similar  coding  scheme  is  under  consider- 
ation for  use  internationally  among  free  world  nations 
(Clauser,  1980] .   The  need  for  such  a  program  was  indi- 
cated by  difficulties  encountered  by  some  foreign  countries 
operating  aircraft  acquired  from  the  United  States  and 
other  foreign  sources. 

In  Figure  1,  it  is  apparent  that  the  height  of  U.S. 
Air  Force  flight  personnel  is  generally  greater  than  that 
of  other  countries.   The  fifth  percentile  in  U.S.  height 
is  greater  than  the  50th  percentile  of  the  samples  from 
Japan,  Thailand,  and  Vietnam.   Figure  2  shows  a  similar 
relationship  for  sitting  height. 
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Figure  1 
Heights  of  Seven  Military  Populations 
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Figure  2 
Sitting  Heights  of  Seven  Military  Populations 
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In  this  case  the  5th  percentile  is  larger  than  the  50th 
percentile  of  Thailand  and  Vietnam.   Aircraft  designed  to 
accommodate  people  as  small  as  the  5th  percentile  of  the  U.S 
population  can  be  expected  to  exclude  about  half  of  certain 
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foreign  populations  and  as  much  as  75  percent  of  the 
Vietnam  aviation  personnel  [Kennedy,  19733. 

Although  design  criteria  require  accommodation  of  the 
central  ninety  percent  of  the  Naval  aviation  population, 
the  existing  population  is  imperfectly  described  by  the  1964 
sample.   Changes  in  the  selection  process,  trends  in  the 
United  States  general  population  were  cited  as  factors  con- 
tributing to  the  differences  between  the  1964  and  the  1969 
populations  TMoroney,  Kennedy,  Gifford,  and  Provost,  197l]. 
In  addition,  designing  to  accommodate  the  fifth  and  ninety- 
fifth  percentiles  in  each  variable  does  not  yield  an 
accommodation  of  ninety  percent  of  the  total  population. 
Since  exclusion  results  from  any  one  measurement  being  out- 
side limits,  accommodation  is  as  much  a  function  of  the 
correlation  between  measurements  as  it  is  a  function  of 
each  measurement.   The  correlation  between  measurements  is 
appreciably  less  than  1.0  and  accommodation  is  less  than 
that  indicated  by  any  univariate  estimate.   When  twelve 
measurements  were  considered,  the  exclusion  of  subjects  who 
were  outside  the  5th  and  95th  percentile  range  in  any 
measurement  yielded  a  cumulative  exclusion  of  48  percent 
[Moroney  and  Smith,  1972].   As  the  existing  population  is 
altered  further  by  expanding  the  acceptable  range  to 
include  a  greater  percentage  of  the  female  population,  user 
accommodation  will  be  adversely  affected.   For  example,  a 
generalized  cockpit  built  to  male  anthropometric  standards 
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excluded  22  percent  of  the  male  user  population.  The  same 
cockpit  excluded  90  percent  of  the  females  in  the  U.S.  Air 
Force  [Ketcham-Weidl  and  Bittner,  19761. 

Before  considering  future  changes  in  the  population, 
the  current  anthropometric  exclusion  situation  should  be 
determined.   However,  the  proportion  of  the  Naval  aviation 
population  excluded  by  the  existing  design  specification 
has  not  been  quantitatively  determined.   Similarly  the 
proportion  excluded  by  each  specific  aircraft  design  is 
unknown.   Some  aircraft  can  be  expected  to  exclude  more 
than  the  ten  percent  designed  out  by  the  specification. 
Other  aircraft  will  exclude  less  than  the  maximum  allowable 
ten  percent  by  virtue  of  design  allowance  or  perhaps  by 
chance. 

The  methods  employed  in  measuring  aircrew  personnel 
have  been  found  to  be  rather  unreliable  [Moroney,  et .  al . , 
1971].   These  measurements,  however,  form  the  basis  of  the 
Anthropometric  Compatibility  Assignment  Program  and  their 
accuracy  is  critical  since  they  constitute  a  determining 
factor  in  an  officer's  career  path  starting  point.   Addition- 
ally, these  same  data  will  be  used  to  help  identify  (or 
eliminate)  subsequent  aviation  duty  assignments.   Consider- 
able attention  has  been  placed  on  determining  the  proper 
use  of  the  data  but  the  actual  collection  of  the  data 
appears  to  be  an  area  where  reliability  could  be  improved. 
Currently,  there  is  no  systematic  method  of  checking 
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anthropometric  measurements.   Each  student  is  measured  once 
without  any  attempt  at  remeasurement  or  checking  of  the 
accuracy  of  the  measurements.   If  a  measurement  is  recorded 
in  error,  the  error  becomes  part  of  the  student's  record, 
undetected  until  some  independent  event  brings  the  error  to 
light.   Errors  in  measurement  of  as  much  as  three  inches 
have  been  found  Co'Leary,  1980 J  subsequent  to  the  erroneous 
measurement  being  used  to  determine  anthropometric  suit- 
ability and  assignment  of  a  student  to  an  aircraft  type. 
Some  students  have  completed  flight  training  before  anthro- 
pometric measurement  errors  were  detected,  resulting  in 
reassignment  and  expensive  retraining. 

The  purpose  of  this  thesis  is  to  quantify  the  amount 
of  incompatibility  between  aircraft  and  aircrewmen  and  to 
propose  some  alternatives  to  improve  the  reliability  of  the 
personal  anthropometric  measurements  now  being  taken.   By 
examining  samples  drawn  from  three  Naval  aviation  popula- 
tions, each  of  a  different  year,  proportions  of  those 
populations  excluded  from  Naval  aircraft  will  be  determined. 
The  exclusions  resulting  from  changes  in  the  population  will 
be  quantified  and  the  specific  proportions  excluded  from 
each  aircraft  type,  as  a  function  of  sample  source,  will  be 
calculated.   Those  results  will  lead  to  a  discussion  of  some 
inconsistent  patterns  in  the  data  and  recommended  changes  in 
the  data  collection  procedure  now  used  by  the  Navy. 
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Chapter  2,  DATA  SOURCES,  describes  the  three  data  sets 
analyzed.   For  each  of  the  three  samples:   the  subjects  were 
described,  previous  statistical  analyses  discussed,  and 
methods  of  data  collection  compared.   Chapter  3,  DATA 
PREPARATION,  describes  the  preparation  of  each  data  set 
to  facilitate  analysis.   Preparations  consisted  primarily 
of  the  use  of  a  segment  of  the  Computerized  Accommodation 
Percentage  Evaluation  (CAPE)  model  to  generate  data  points 
from  a  statistical  description  of  the  population.   Chapter 
4,  ANALYSIS  OF  DATA,  describes  the  screening  of  each  data 
set  to  determine  the  proportion  of  the  sample  excluded. 
Screening  involves  the  elimination  of  any  subject  with  an 
anthropometric  measurement  outside  prescribed  limits.   The 
proportions  excluded  from  each  sample  are  then  compared. 
Observations  on  the  apparent  shortcomings  in  data  collection 
methods,  errors  in  measurements,  and  unexplained  inconsis- 
tencies were  included.   Chapter  5,  SUMMARY  AND  CONCLUSIONS, 
includes  inferences  concerning  interpretations  of  anthropo- 
metric data  and  recommendations  in  the  area  of  data 
collection  methods  for  use  by  the  Navy. 


20 


II.   DATA  SOURCES 

The  data  to  be  analyzed  consisted  of  samples  from  three 
populations.   Samples  of  Naval  aviation  personnel  were  taken 
in  1964,  1969,  and  1975.   Differences  between  populations 
other  than  the  year  of  the  sample  exist  and  they  will  be 
discussed  below. 

A.   THE  1964  SURVEY  DATA 

In  1964,  a  survey  was  conducted  to  determine  the  anthro- 
pometric characteristics  of  personnel  flying  in  Naval 
aircraft.   The  survey  was  conducted  in  a  manner  to  yield 
results  sufficiently  precise  to  be  utilized  subsequently  for 
design  of  aircraft.   The  Bureau  of  Medicine  and  Surgery, 
tasked  to  collect  data  for  use  by  the  Bureau  of  Weapons, 
sent  a  trained  anthropometric  measuring  team  to  ten  different 
Naval  and  Marine  Corps  air  stations  to  take  measurements. 

Included  in  the  sample  were  pilots,  Naval  Aviation 
Observers,  bombardier-navigators,  radar  observers,  flight 
surgeons  who  were  designated  Naval  Aviators,  and  enlisted 
personnel  who  were  permanently  designated  aeronautical 
personnel  and  who  were  still  eligible  for  assignment  to 
subsequent  duty  involving  flying.   Naval  Aviators  included 
both  U.S.  Navy  and  the  U.S.  Marine  Corps  flyers. 

The  subjects  in  this  survey  comprised  approximately  ten 
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percent  of  the  Naval  aviation  population  on  active  duty  in 
1964.   The  sample  included  personnel  of  various  ages,  ranks, 
ethnic  backgrounds  and  sources  (Naval  Academy,  NROTC,  etc.). 
In  that  respect  the  sample  undoubtedly  provided  an  accurate 
representation  of  the  Naval  aviation  population.   The  com- 
position of  the  sample  is  somewhat  different  from  the  later 
samples  in  the  respect  that  enlisted  men  and  older  men  were 
surveyed  in  1964  but  not  in  the  later  samples.   The  relia- 
bility of  the  measurements  received  special  attention  and 
was  enhanced  by  the  repeat  measurement  of  randomly  selected 
subjects,  insuring  rapid  detection  of  measurement  errors  and 
encouraging  the  data  collectors  to  be  careful  in  their  task. 
Each  measurement  was  taken  by  personnel  specifically  trained 
for  the  job,  using  a  technique  of  measurement  recognized  as 
the  most  accurate  available. 

Published  results  of  the  survey  include  means  and 
standard  deviations,  as  well  as  percentile  listings  for  each 
variable  CGifford,  Provost,  and  Lazo ,  19651.   Subsequent 
analysis  of  the  data  led  to  publication  of  correlations 
between  pairs  of  variables  CMoroney,  et .  al . ,  1971]. 

The  results  of  this  survey  were  adopted  by  the  Navy  as 
the  basis  for  aircraft  design  specifications.   MIL-STD-1333A, 
Military  Standard  Aircrew  Station  Geometry  for  Military 
Aircraft,  includes  NAEC-ACEL-553  as  part  of  the  standard, 
thus  specifying  the  1964  data  as  the  guideline  for  aircraft 
being  built  today  (NAEC-ACEL-553  reported  the  results  of  the 
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1964  survey).  In  that  sense,  these  results  are  recognized 
as  the  standard  and  they  have  become  the  yardstick  against 
which  other  data  can  be  compared. 

This  data  set  has  two  shortcomings  with  respect  to  this 
thesis,  however.   No  measurements  of  leg  length  (or  buttock- 
leg  length)  were  taken.   The  components  of  leg  length 
(buttock-knee  length,  knee  height  sitting  and  popliteal 
height  sitting)  were  recorded  but  leg  length  had  not  yet 
been  recognized  as  a  critical  selection  dimension.   Also, 
the  results  of  data  analysis  are  available,  but  the  raw 
data  are  not.   The  importance  of  these  two  limitations  will 
become  apparent  when  procedures  are  discussed. 

B.   THE  1969  DATA 

Subsequent  to  the  introduction  of  an  integrated  anthro- 
pometric measuring  device  in  1964,  the  Navy  required  that 
the  device  be  used  to  measure  all  aviation  training  candi- 
dates during  their  aviation  training  entrance  physical 
examination,  which  was  administered  at  the  School  of  Aviation 
Medicine  (now  the  Naval  Aerospace  Medical  Institute).   The 
collection  of  anthropometric  data  was  thus  institutionalized 
and  became  a  routine  process.   Eight  measurements  were  taken 
on  each  entering  Student  Naval  Aviator  and  Student  Naval 
Flight  Officer.   These  measurements  were: 

weight  trunk  height 

height  buttock-knee  length 

sitting  height  buttock-leg  length 

shoulder  width  functional  reach 
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Between  January  1966  and  August  1969,  6,534  students 
were  measured.   Their  measurements  comprise  the  data  set 
herein  called  the  1969  data.   The  description  of  this  data 
set  was  published  in  1971  LMoroney,  et .  al . ,  1971 1   and 
includes  the  mean  and  standard  deviation  of  six  measurements 
(shoulder  height  and  functional  reach  were  deleted  from  the 
analysis),  percentile  listings  for  those  measurements  and 
correlations  between  variables. 

No  information  on  functional  reach  or  shoulder  height 
was  published  because  of  the  discovery  of  measurement  tech- 
nique errors  resulting  in  measurements  uncorrectible  by 
any  means  short  of  repeat  measurement.   This  deviation  from 
ideal  procedure  should  not  have  been  altogether  surprising 
since  measurements  were  taken  by  personnel  whose  level  of 
training  did  not  match  that  of  the  1964  data  collectors. 
In  addition,  measurements  were  taken  and  recorded  with  no 
subsequent  checking,  leaving  procedural  errors  undetected 
and  biasing  the  measurements  of  many  subjects. 

The  size  of  the  1969  data  set  is  quite  large,  consisting 
of  virtually  all  of  the  student  aviation  input  for  over  three 
years.   As  with  the  1964  data  set,  no  raw  data  is  available, 
only  analytical  results.   The  absence  of  functional  reach 
from  the  data  analysis  will  present  an  obstacle  to  be 
discussed  in  Chapter  3. 
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C.   THE  1975  DATA 

Following  the  publication  of  the  results  of  the  1969 
data  analysis,  an  improved  integrated  anthropometric 
measuring  device  was  installed  at  NAM I  and  it  was  subse- 
quently used  to  measure  incoming  students.   That  device  is 
still  in  use  today.   It  is  operated  by  NAMI  staff  Hospital 
Corpsmen  and  it  was  used  to  measure  the  same  eight  dimen- 
sions as  was  its  predecessor.   A  sample  of  968  students 
who  entered  training  in  1975  comprise  the  subjects  of  the 
1975  data.   The  raw  data  is  available  and  from  an  examina- 
tion of  the  data  points  thirty-two  subjects  were  found  to 
be  deficient  in  some  respect  and  did  not  meet  the  entry 
requirements  delineated  by  the  Bureau  of  Medicine  and  Surgery 
LMANMED,  1978].   Most  of  these  ineligibilities  were  based  on 
excessive  weight.   However,  all  subjects  were  screened  for 
minimum  and  maximum  values  of  weight,  height,  sitting  height, 
and  buttock-leg  length,  as  well  as  for  height-weight  rela- 
tionship.  The  remaining  936  subjects  comprise  the  acceptable 
sample. 

The  1975  data  were  collected  in  a  manner  similar  to  that 
used  for  the  1969  data  set :   a  routine  execution  of  estab- 
lished procedures  by  Navy  enlisted  personnel.   The  1975 
data  also  manifest  data  collection  deviations  similar  to 
those  of  the  1969  data.   The  number  of  subjects  is  smaller 
than  either  of  the  other  samples  but  still  large  enough  to 
lend  itself  to  statistical  analysis. 
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III.   DATA  PREPARATION 

The  goal  of  the  analysis  was  to  determine  the  proportions 
of  the  1964,  1969,  and  1975  samples  excluded  by  the  four 
anthropometric  limitations  (sitting  height,  functional 
reach,  buttock-knee  length,  and  leg  length)  required  by 
OPNAVINST  3710. 36A.   This  goal  was  achieved  by  screening 
each  subject  four  times,  once  for  each  of  his  measurements, 
and  eliminating  every  subject  whose  measurements  were  not 
all  within  limits. 

This  procedure  required  knowledge  of  the  four  individual 
measurements  for  each  subject.   In  the  1975  data  set,  the 
four  measurements  (and  more)  were  readily  available.   For 
the  earlier  data  sets,  only  means,  standard  deviations,  and 
correlations  were  available  and  no  individual  measurements 
were  known.   It  was  necessary,  therefor,  to  generate  data 
to  simulate  those  measurements. 

The  Computerized  Accommodated  Percentage  Evaluation 
(CAPE)  model  (Bittner,  19751  was  used  to  generate  the 
required  data  points.   This  model  required  as  entering 
arguments  the  mean  and  standard  deviation  for  each  of  the 
four  variables,  as  well  as  a  matrix  of  correlations  between 
variables.   The  model  used  a  Monte  Carlo  process  to  generate 
data  points.   A  program  listing  is  included  in  Appendix  C 
(Appendix  D  presents  the  same  program  modified  for  use  with 
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the  IBM  360-67  computer  and  other  machines  using  Fortran  IV) . 
The  CAPE  model  has  been  incorporated  in  the  Computerized 
Accommodated  Reach  Model  and  is  available  at  the  Naval 
Aviation  Development  Center,  Warminster,  Pennsylvania. 

Unfortunately,  some  of  the  values  required  for  generating 
the  data  points  were  missing  from  the  survey  reports. 
Specifically,  the  1964  data  included  no  information  on  leg 
length  and  the  1969  data  omitted  functional  reach.   It  was 
necessary  to  find  values  in  similar  populations  and  substitute 
them  for  those  missing.   Taking  values  from  other  populations 
required  the  assumption  that  concurrent  populations,  similar 
in  cultural,  racial,  and  demographic  composition,  and  having 
passed  similar  anthropometric  entrance  requirements,  will 
have  comparable  anthropometric  features.   That  assumption 
was  made  and  surrogate  values  for  the  missing  means  and  stan- 
dard deviations  were  taken  from  populations  which  matched 
the  descriptions  of  the  populations  under  analysis  (1964  and 
1975),  both  subjectively  and  numerically.   Inter-variable 
correlations  were  determined  on  much  the  same  basis:   substi- 
tuting correlations  from  similar  populations.   Ketcham-Weidl 
and  Bittner  L19762    had  substituted  correlations  from  a  data 
set  of  male  subjects  into  a  set  of  female  data  and  exper- 
ienced little  loss  of  precision. 

For  the  1964  data,  the  mean  and  standard  deviation  for 
leg  length  were  assumed  to  be  43.82  and  2.01  inches,  respec- 
tively.  Those  values  were  established  in  a  1968  survey  of 
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U.S.  Army  helicopter  students  CSchane,  Littell  and  Moultree, 
1969].   Of  the  information  available,  the  population  of  that 
survey,  better  than  any  other,  approximated  the  1964  Naval 
aviation  population  with  respect  to  the  factors  which  Moroney 
Cl97ll  reported  as  influencing  differences  between  the  1964 
and  1969  Navy  data  sets.   The  Sitting  height-Leg  length  cor- 
relation maintained  the  same  value  (0.45)  for  the  1969  and 
1975  data  sets  so  the  same  value  was  used  for  1964.   Func- 
tional reach-Leg  length  and  Buttock  knee-Leg  length 
correlations  (0.639  and  0.860)  were  extracted  from  a  1970 
survey  of  U.S.  Army  aviators  [Churchill,  McConville,  Laubach, 
and  White,  1971],  again  because  of  the  similarities  in 
descriptions  of  the  1970  Army  and  1964  Navy  populations. 
For  the  1969  data,  Functional  reach  values  for  USAF 
flight  personnel  reported  by  Churchill  Cl97l]  were  used  for 
the  mean  and  standard  deviation  (31.24  and  1.62  inches). 
The  sitting  height-functional  reach  and  the  buttock  knee- 
functional  reach  correlations  were  assumed  to  be  the  same  as 
those  found  in  the  1975  data,  0.44  and  0.63,  respectively. 
These  values  were  greater  than  those  of  the  1964  data  and 
their  use,  rather  than  the  1964  values,  would  influence 
error  in  the  conservative  direction,  that  is,  in  the  direc- 
tion of  lesser  difference  between  univariate  and  multi- 
variate exclusion.   The  leg  length-functional  reach  corre- 
lation was  taken  from  Churchill's  1970  survey  and  is  the 
same  value  assumed  for  the  1964  data. 
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Table  I 

CAPE  INPUTS  AND  OUTPUTS 
FOR  THE  1964  DATA 

Means,  Standard  Deviations  and  Ranges 


Variable 


Observed 
Data 


Generated 
Data 


Sitting  Height 
mean 
std.dev. 
range 

Functional  Reach 
mean 
std.dev. 
range 

Buttock-Knee  Length 
mean 
std. dev. 
range 

Leg  Length2- 
mean 
std. dev. 
range 


36.28 
1.25 
32.19-41.62 


31.51 
1.42 
27.26-36.31 


24.09 
1.00 
20.73-27.81 


43.82 
2.01 


36.29 
1.26 
31.81-40.25 


31.47 
1.40 
26.35-35.38 


24.10 
9.94 
21.06-27.11 


43.82 
1.99 
37.97-49.87 


aSchane,  W.P.,  Littell,  D.E.',  and  Moultree,  C.G., 
Selected  Anthropometric  Measurements  of  1,640  U.S.  Army 
Warrant  Officer  Candidate  Flight  Trainees,  U.S.  Army 
Aeromedical  Research  Laboratory  Report  69-2,  as  reported 
in :   Anthropometric  Source  Book  Volume  II:   A  Handbook  of 
Anthropometric  Data,  NASA  Reference  Publication  1024, 
1978,  p.  156. 
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Table  II 

CAPE  INPUTS  AND  OUTPUTS 
FOR  THE  1964  DATA 

Correlations 


Sitting    Functional  Buttock-    Leg 
Height     Reach      Knee  Length  Length 


Observed  Data 


Sitting  Height 
Functional  Reach 
Buttock-Knee  Length 
Leg  Length 


1.0 


.376 
.0 

.382 
.586 
1.0 

.45a 
,639b 
.86b 
1.0 

Generated  Data 


Sitting  Height 
Functional  Reach 
Buttock-Knee  Length 
Leg  Length 


1.0 


.410 

.374 

.446 

1.0 

.587 

.650 

1.0 

.867 
1.0 

aSame  as  observed  value  for  1969  and  1975  samples. 

^Churchill,  E. ,  McConville,  J.T.,  Laubach,  L. ,  and 
White,  R.M.,  Anthropometry  of  U.S.  Army  Aviators-1970 ,  U.S. 
Army  Natick  Laboratories  Technical  Report  72-52-CE,  December 
1971,  p.  278. 
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Table  III 

CAPE  INPUTS  AND  OUTPUTS 
FOR  THE  1969  DATA 

Means,  Standard  Deviations,  and  Ranges 


Observed 


Generated 


Variable 

Data 

Data 

Sitting  Height 
mean 
std.dev. 
range 

36.76 
1.21 
31.70-41. 

,0 

36.77 
1.22 
32.43-40, 

,61 

Functional  Reac 
mean 
std. dev. 
range 

:ha 

31.24 
1.62 

31.20 
1.60 
20.02-35. 

,73 

Buttock-Knee  Length 
mean 
std. dev. 
range 

24.54 
1.26 
20.40-29. 

.90 

24.55 
1.25 
20.67-28. 

,31 

Leg  Length 
mean 
std. dev . 
range 

43.86 
2.08 
32.50-50. 

.70 

43.86 
2.05 
38.14-50, 

,02 

aChurchill,  E. ,  McConville,  J.T.,  Laubach,  L. ,  and 
White,  R.M.,  Anthropometry  of  U.S.  Army  Aviators-1970 ,  U.S. 
Army  Natick  Laboratories  Technical  Report  72-52-CE,  December 
1971,  p.  91. 
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Table  IV 

CAPE  INPUTS  AND  OUTPUTS 
FOR  THE  1969  DATA 

Correlations 


Sitting    Functional  Buttock-    Leg 
Height     Reach      Knee-Length  Length 


Observed  Data 


Sitting  Height 
Functional  Reach 
Buttock-Knee  Length 
Leg  Length 


1.0 


.44a 

.40 

.45 

1.0 

.63a 

.639b 

1.0 

.86 
1.0 

Generated  Data 


Sitting  Height 
Functional  Reach 
Buttock-Knee  Length 
Leg  Length 


1.0 


.473 

.393 

.449 

1.0 

.629 

.652 

1.0 

.801 
1.0 

aSame  as  observed  value  in  1975. 

bChurchill,  E.,  McConville,  J.T.,  Laubach,  L. ,  and 
White,  R.M.,  Anthropometry  of  U.S.  Army  Aviators-1970 ,  U.S. 
Army  Natick  Laboratories  Technical  Report  72-52-CE,  December 
1971,  p.  278. 
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Table  V 
DATA  SET  COMPARISON 

Means,  Standard  Deviations,  and  Ranges 


Variable  1964  1969  1975 

Sitting  Height 

mean  36.28  36.76  36.61 

std.dev.        1.25         1.21  1.15 

range  32.19-41.62  32.43-40.61  32.90-40.00 

Functional  Reach 

mean  31.51  31.24  30.89 

std.dev.        1.42         1.62  1.16 

range  27.26-36.31  26.02-35.73  27.10-34.60 

Buttock-Knee  Length 

mean  24.09  24.54  24.05 

std.dev.         1.00         1.26  1.07 

range  20.73-27.81  20.67-28.31  20.70-26.80 

Leg  Length 

mean  43.82  43.86  42.89 

std.dev.        2.01         2.08  1.89 

range  37.97-49.87  38.15-50.02  37.10-48.20 
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Table  VI 
DATA  SET  COMPARISON 
Correlations 


1964 

Sitting 

Functional 

Buttock- 

Leg 

Height 

Reach 

Knee  Length 

Length 

Sitting  Height 

1.0 

.376 

.382 

.45 

Functional  Reach 

1.0 

.586 

.639 

Buttock-Knee  Length 

1.0 

.86 

Leg  Length 

1.0 

1969 

Sitting  Height 

1.0 

.44 

.40 

.45 

Functional  Reach 

1.0 

.63 

.639 

Buttock-Knee  Length 

1.0 

.86 

Leg  Length 

1.0 

1975 

Sitting  Height 

1.0 

.44 

.40 

.45 

Functional  Reach 

1.0 

.63 

.61 

Buttock-Knee  Length 

1.0 

.79 

Leg  Length 

1.0 
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When  the  assumed  values  described  above  were  included 
in  the  data  analysis,  the  requirements  for  CAPE  input  were 
complete.   The  inputs  to  CAPE  are  listed  in  Tables  I 
through  IV  under  Observed  Data.   One  thousand  subjects  were 
generated  for  the  1964  and  for  the  1969  data  sets.   Together 
with  the  1975  data  in  its  original,  as  observed  form,  the 
result  was  three  similar  data  sets  with  sample  sizes  of 
1000,  1000,  and  936.   The  validity  of  this  model's  perfor- 
mance has  been  demonstrated  with  a  correlation  between 
empirical  and  generated  results  of  0.997  to  0. 999  tBittner , 
19743.   For  the  1964  and  1969  data  sets,  the  CAPE  generated 
output  was  analyzed  to  insure  that  the  means,  standard 
deviations,  and  correlations  of  the  generated  data  were 
actually  those  sought.   The  values  used  as  inputs  to  CAPE 
and  those  resulting  from  analysis  of  the  generated  data  are 
shown  in  Tables  I  through  IV.   Tables  V  and  VI  compare  the 
summary  statistics  for  all  three  data  sets. 
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IV.   ANALYSIS  OF  DATA 

As  a  result  of  the  data  generating  procedure,  all  three 
data  sets  included  a  value  for  each  of  the  four  relevant 
measurements  for  each  subject.   The  samples  were  then 
examined  more  closely  to  determine  which  subjects  in  each 
sample  would  be  excluded,  either  by  design  constraints  (5th 
to  95th  percentiles  or  3rd  to  98th  percentiles)  or  by  air- 
craft constraints  specified  by  NAVAIRINST  3710.9. 

A.   PERCENTAGES  EXCLUDED 

Fortran  program  SCREEN  was  used  to  count  the  subjects 
who  were  too  large  or  too  small  to  satisfy  anthropometric 
requirements.   For  each  of  the  three  samples,  identical 
screening  operations  were  conducted.   Each  subject  was 
examined  and  excluded  if  any  one  or  more  of  the  four 
measurements  under  consideration  was  outside  prescribed 
limits.   Each  subject's  anthropometric  features  were  tested 
against  each  cutoff  limit.   Thus  an  individual  might  be 
excluded  on  both  functional  height  and  sitting  height  but 
not  on  buttock-knee  length  and  leg  length.   For  each  elimin- 
ation encountered,  the  subject  was  counted  as  being 
eliminated  and  each  variable  causing  an  elimination  was 
tallied.   At  the  end  of  each  sample,  the  resultant  totals 
were  converted  to  percentages  of  the  total  sample,  yielding 
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a  percentage  eliminated  by  each  of  the  four  variables  and 
a  total  exclusion  resulting  from  the  net  effect  of  the  four 
variables'  eliminations. 

The  prescribed  limits  mentioned  above  were  set  equal  to 
the  percentile  values  in  the  design  specification  (5th  and 
95th  or  3rd  and  98th  percentiles).   Then  the  limits  were 
set,  successively,  to  the  anthropometric  limits  of  each 
aircraft  for  which  anthropometric  coding  had  been  estab- 
lished CNAVAIRINST  3710.9]. 

1 .   Design  Constraints 

To  screen  the  samples  by  design  specification,  the 
subjects  were  screened  to  eliminate  those  with  any  of  the 
four  measurements  less  than  the  5th  percentile  or  above  the 
95th  percentile  of  the  1964  data,  the  design  standard  CMIL- 
STD-1472BJ.   Then  the  screening  cutoff  points  were  changed 
to  correspond  to  the  3rd  and  98th  percentiles  CNAVAIR  SD-24K, 
1973^,  and  the  screening  process  was  repeated. 

It  is  interesting  to  note  at  this  point  that  percen- 
tile values  for  leg  length  of  the  1964  data  do  not  exist 
since  that  measurement  was  not  taken  in  the  1964  survey. 
Still,  the  Navy  considers  leg  length  to  be  a  critical  dimen- 
sion in  that  it  is  one  of  the  measurements  included  in  the 
Anthropometric  Compatibility  Assignment  Program.   Critical 
dimensions  are  also  the  basis  of  the  design  specification 
since  '...design  limits  shall  be  based  upon  a  range  from  the 
5th  percentile  to  the  95th  percentile  values  for  critical 
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body  dimensions'  CMIL-STD-1472B3 .   Since  the  percentile 
values  were  not  available,  values  used  for  leg  length 
percentiles  were  determined  from  the  CAPE  generated  sample 
of  1000  subjects,  described  in  Chapter  III,  DATA  PREPARATION. 
The  design  specification  percentile  screens  resulted 
in  the  maximum  exclusion  percentage  to  be  expected  if  air- 
craft were  designed  to  accommodate  only  the  specified 
percentiles,  and  no  more,  for  each  dimension.   This  result 
amounts  to  a  worst  case  exclusion  and  could  occur  only  if 
a  design  exactly  met  the  bare  minimum  accommodation  in  each 
variable,  a  very  unlikely  event,  and  one  which  would  still 
not  meet  the  requirement  to  accommodate  90  percent  of  the 
user  population  as  specified  by  MIL-STD-1472B.   The  percen- 
tile screen  results,  shown  in  Table  VII,  show  that  cumulative 
exclusion  can  greatly  exceed  the  maximum  single  variable 
elimination.   For  example,  the  1969  data  screened  on  the 
5th  to  95th  percentiles  shows  a  total  exclusion  of  40.4 
percent  but  the  highest  single  variable  elimination  is  the 
21  percent  excluded  by  buttock-knee  length. 
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Table  VII 

Design  Specification  Exclusions  for 
the  5-95  Percentiles  and  the  3-98  Percentiles 


Variable 


1964 


1969 


1975 


5th  to  95th  Percentile 
Percentage  Excluded 


Sitting  Height 

11.1 

11.5 

7.7 

Functional  Reach 

8.7 

14.6 

8.4 

Buttock-Knee  Length 

9.1 

21.0 

13.2 

Leg  Length 

9.8 

11.3 

11.8 

Cumulative 

26 

4 

40 

4 

26.7 

3rd  to  98th  Percentile 
Percentage  Excluded 


Sitting  Height 

5.0 

4.7 

2.7 

Functional  Reach 

4.8 

9.4 

5.9 

Buttock-Knee  Length 

5.4 

11.6 

6.4 

Leg  Length 

4.9 

4.7 

6.9 

Cumulative 

14 

3 

22 

7 

14.3 
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The  difference  between  cumulative  exclusion  and  the 
largest  univariate  exclusion  figure,  which  could  be  taken 
falsely  as  an  upper  bound,  is  due  to  the  non-unity  corre- 
lation between  variables.   Although  some  anthropometric 
models  have  described  all  anthropometric  dimensions  as  a 
specific  multiple  of  height  £Roozbazar,  19793  with  the 
implicit  assumption  of  intervariable  correlation  equal  to 
1.00,  those  models  fail  to  recognize  the  variable  propor- 
tions of  people  and  would  lead  to  the  acceptance  of  the 
largest  univariate  exclusion  as  the  net  cumulative  exclu- 
sion.  The  assumption  of  a  correlation  of  1.00  amounts  to 
treating  the  population  as  if  the  man  with  95th  percentile 
height  must  also  have  95th  percentile  leg  length,  95th 
percentile  shoulder  width,  and  95th  percentile  in  all 
measurements,  i.e.,  'the  95th  percentile  man.'   If  that  were 
the  case,  the  percentile  screens  would  have  eliminated  the 
same  subjects  on  all  variables,  as  well  as  on  the  cumulative 
exclusion.   If,  at  the  other  extreme,  intervariable  corre- 
lation were  zero,  the  single  variable  eliminations  would  be 
independent  and  their  arithmetic  sum  would  appear  as  the 
cumulative  exclusion. 

In  reality,  some  subjects  with  leg  lengths  shorter 
than  the  5th  percentile  also  have  sitting  heights  less  than 
the  5th  percentile.   Others  have  short  legs  but  sitting 
heights  nearer  the  mean.   The  evidence  of  variable  propor- 
tions along  with  variable  sizes  is  found  in  the  cumulative 
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exclusion  values.   Generally,  the  cumulative  exclusion  is 
less  than  the  sum  of  the  single  variable  eliminations  but  is 
greater  than  any  of  the  single  value  eliminations.   Only 
rounding  error  could  cause  the  cumulative  exclusion  to 
exceed  the  sum  of  the  single  variable  eliminations.   The 
cumulative  exclusion  is  generally  as  small  as  the  largest 
single  variable  when  only  one  variable  resulted  in  any 
elimination. 

2.   Aircraft  Constraints 

The  specific  aircraft  screens  were  accomplished  with 
the  same  technique  as  used  on  the  design  specification 
screens.   The  cut  points  for  each  variable  were  the  values 
corresponding  to  each  aircraft's  anthropometric  restrictions 
specified  by  NAVAIRINST  3710.9.   Then,  for  each  aircraft, 
the  three  samples  were  examined  to  eliminate  subjects  whose 
anthropometric  measurements  were  large  enough,  or  small 
enough,  to  be  excluded  by  the  aircraft's  restrictions.   The 
number  of  subjects  eliminated  by  each  measurement  and  the 
cumulative  exclusions  were  totaled  and  converted  to  a  percen- 
tage of  the  sample.   Different  aircraft  with  the  same 
anthropometric  restrictions  were  consolidated  into  eighteen 
groups  and  the  screens  were  executed,  resulting  in  the  1964, 
1969,  and  1975  samples  being  matched  to  all  aircraft  types. 
The  percentage  exclusions  are  shown  in  Table  VIII,  which 
lists  aircraft  types  in  alphabetical  order. 
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Table  VIII 
Exclusion  by  Aircraft  Type 


Code 

Variable 

Excl 

usion  Percent 

age 

Aircraft 

1964 

1969 

1975 

A-3 

sit .ht . 

9 

fnct . rch. 

0.6 

2.6 

0.3 

- 

b-knee  lnth. 

- 

- 

- 

0 

leg  lnth. 

0.1 

0.0 

0.1 

cumulative 

0. 

6 

2.6 

0.4 

A-4E 

089 

sit .ht . 

0.9 

1.7 

0.4 

A-4F 

9 

fnct . rch. 

0.6 

2.6 

0.3 

789 

b-knee  lnth. 

1.0 

6.6 

0.6 

0 

leg  length 

0.1 

0.0 

0.1 

cumulative 

2. 

5 

10.3 

1.5 

A-4M 

0 

sit . ht . 

0.2 

0.1 

0.0 

9 

fnct . rch. 

0.6 

2.6 

0.3 

789 

b-knee  lnth. 

1.0 

6.6 

0.6 

0 

leg  length 

0.1 

0.0 

0.1 

cumulative 

1. 

8 

9.3 

1.1 

A-6A 

06789 

sit . ht . 

5.0 

9.6 

4.0 

A-6B 

- 

fnct . rch. 

- 

- 

- 

A-6E 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

- 

leg  length 

- 

- 

- 

cumulative 

5. 

0 

9.9 

4.0 

A-7A 

089 

sit . ht . 

0.9 

1.7 

0.4 

A-7C 

6789 

fnct . rch. 

7.4 

13.6 

8.8 

A-7E 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

0 

leg  length 

0.1 

0.0 

0.1 

cumulative 

8. 

2 

15.6 

9.2 

AV-8A 

09 

sit . ht . 

0.4 

0.8 

0.1 

9 

fnct . rch. 

0.6 

2.6 

0.3 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

— 

leg  length 

- 

- 

- 

cumulative 

1. 

0 

3.8 

0.4 

C-118 

_ 

sit . ht . 

_ 

_ 

_ 

C-131 

- 

fnct . rch. 

- 

- 

- 

- 

b-knee  lnth. 

— 

- 

- 

0 

leg  length 

0.1 

0.0 

0.1 

cumulative 

0. 

1 

0.7 

1.6 
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Table  VIII  (Continued) 
Exclusion  by  Aircraft  Type 


Code 

Variable 

Ex 

:cl 

usion  Percent 

age 

Aircraft 

1964 

1969 

1975 

CH-53 

_ 

sit . ht . 

^ 

_ 

_ 

— 

fnct . rch. 

- 

- 

- 

- 

b-knee  lnth. 

— 

- 

- 

01 

leg  length 

0.7 

0.7 

1.6 

cumulative 

0. 

7 

0.7 

1.6 

EA-3 

_ 

sit . ht . 

_ 

_ 

_ 

9 

fnct . rch. 

0.6 

2.6 

0.3 

— 

b-knee  lnth. 

- 

- 

- 

0 

leg  length 

0.1 

0.0 

0.1 

cumulative 

0. 

6 

2.6 

0.4 

EA-6B 

06789 

sit . ht . 

5.0 

9.6 

4.0 

- 

fnct . rch. 

- 

- 

- 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

— 

leg  length 

- 

- 

- 

cumulative 

5. 

0 

9.9 

4.0 

EA-6B 

06789 

sit . ht . 

5.0 

9.6 

4.0 

(LF) 

9 

fnct . rch. 

0.6 

2.6 

0.3 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

— 

leg  length 

- 

- 

- 

cumulative 

5. 

6 

12.5 

4.3 

F-4B(F) 

0789 

sit . ht . 

2.0 

3.9 

1.7 

F-4J(F) 

89 

fnct . rch. 

1.8 

4.3 

1.9 

F-4N(F) 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

F-4S(F) 

0 

leg  length 

0.1 

0.0 

0.1 

cumulative 

3. 

8 

8.4 

3.6 

F-4B(R) 

89 

sit . ht . 

0.7 

1.6 

0.4 

F-4J(R) 

- 

fnct . rch. 

- 

- 

- 

F-4N(R) 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

F-4S(R) 

- 

leg  length 

- 

- 

- 

cumulative 

0. 

7 

2.0 

0.4 
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Table  VIII  (Continued) 
Exclusion  by  Aircraft  Type 


Code 

Variable 

Excl 

usion  Percent 

age 

Aircraft 

1964 

1969 

1975 

F-8D 

0789 

sit .ht . 

2.0 

3.9 

1.7 

F-8J 

- 

fnct .rch. 

— 

- 

- 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

- 

leg  length 

- 

- 

- 

cumulative 

2. 

0 

4.2 

1.7 

F-14(F) 

0 

sit . ht . 

0.2 

0.1 

0.0 

9 

fnct . rch. 

0.6 

2.6 

0.3 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

- 

leg  length 

- 

- 

- 

cumulative 

0. 

8 

3.1 

0.3 

F-14(R) 

_ 

sit . ht . 

_ 

_ 

_ 

— 

fnct . rch. 

- 

- 

- 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

- 

leg  length 

- 

- 

- 

cumulative 

0. 

0 

0.4 

0.0 

H-2 

_ 

sit . ht . 

_ 

_ 

_ 

89 

fnct . rch. 

1.8 

4.3 

1.9 

- 

b-knee  lnth. 

- 

- 

- 

0 

leg  length 

0.1 

0.0 

0.1 

cumulative 

1. 

8 

4.3 

1.9 

OV-10 



sit . ht . 

_ 

_ 

_ 

9 

fnct . rch . 

0.6 

2.6 

0.3 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

- 

leg  length 

- 

- 

- 

cumulative 

0. 

6 

3.0 

0.3 

RA-5C 

09 

sit . ht . 

0.4 

0.8 

0.1 

(F) 

- 

fnct . rch. 

- 

- 

- 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

0 

leg  length 

0.1 

0.0 

0.1 

cumulative 

0. 

5 

1.2 

0.2 

RF-4B 

0789 

sit . ht . 

2.0 

3.9 

1.7 

(F) 

89 

fnct . rch. 

1.8 

4.3 

1.9 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

0 

leg  length 

0.1 

0.0 

0.0 

cumulative 

3. 

8 

8.4 

3.6 
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Table  VIII  (Continued) 
Exclusion  by  Aircraft  Type 


ift  Code 

Variable 

Excl 

usion  Percent 

age 

Aircrs 

1964 

1969 

1975 

RF-4B 

89 

sit . ht . 

0.7 

1.6 

0.4 

(R) 

- 

fnct . rch. 

- 

- 

- 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

- 

leg  length 

- 

- 

- 

cumulative 

0.7 

2.0 

0.4 

RF-8 

0789 

sit . ht . 

2.0 

3.9 

1.7 

- 

fnct . rch. 

- 

- 

- 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

- 

leg  length 

- 

- 

- 

cumulative 

2.0 

4.2 

1.7 

S-3 

_ 

sit .ht . 

_ 

_ 

_ 

- 

fnct . rch. 

- 

- 

— 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

- 

leg  length 

- 

- 

- 

cumulative 

0.0 

0.4 

0.0 

T-2C 

0 

sit . ht . 

0.2 

0.1 

0.0 

TA-4 

9 

fnct . rch. 

0.6 

2.6 

0.3 

TAV-8 

9 

b-knee  lnth. 

0.0 

0.4 

0.0 

0 

leg  length 

0.1 

0.0 

0.1 

cumulative 

0.8 

3.1 

0.4 

U-16 

_ 

sit . ht . 

_ 

_ 

_ 

- 

fnct . rch. 

- 

- 

- 

— 

b-knee  lnth. 

- 

- 

- 

0 

leg  length 

0.1 

0.0 

0.1 

cumulative 

0.1 

0.0 

0.1 

NOTE: 

(F)  - 

(R)  - 
(LF)  - 

Front  seat 

Rear  seat 

Left  front  seat 
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B .   RESULTS 

The  design  criterion  requiring  accommodation  of  the  5th 
and  95th  percentiles  results  in  exclusion  of  over  26  percent 
of  the  population.   That  result  fails  to  meet  the  goal  of 
90  percent  accommodation.   Designing  to  the  3rd  and  98th 
percentiles,  instead  of  the  5th  and  95th,  reduced  exclusion 
to  14  percent.   That  result  still  does  not  meet  the  90  per- 
cent accommodation  goal. 

No  aircraft  type's  coding  eliminated  as  much  as  ten  per- 
cent of  the  1964  sample.   The  A-7  exclusion  was  8.2  percent 
and  was  the  worst  case.   The  accommodation  by  all  aircraft 
of  more  than  90  percent  would  seem  to  indicate  that  the 
cockpit  designs  included  a  margin  of  accommodation  allowance 
in  addition  to  the  minimum  required.   Such  a  safety  factor 
in  design  is  not  unexpected  but  is  still  technically  an 
unrequired  bonus. 

The  exclusion  percentage  was  not  constant,  but  changed 
from  one  sample  to  the  next.   The  greatest  exclusions  were 
found  in  the  1969  data  while  the  1964  and  1975  cumulative 
exclusions  were  about  the  same.   Although  cumulative  exclu- 
sions were  nearly  equal,  the  contributing  factors  had 
varying  proportions.   Sitting  height  eliminations  decreased 
and  eliminations  due  to  buttock-knee  length  and  leg  length 
increased. 

The  1969  data  exhibited  other  features  of  incongruity 
which  detracted  from  the  credibility  of  the  three  data  sets 
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Figure  3 

Comparison  of  Selected  Variables 
and  Standard  Deviations 


Variable 


Inches 


Year 


1964  1969  1975 


73- 

72- 

•  • 

71- 

Height 

70- 

it 

69- 

68- 

67- 

* 

46- 

45- 

Leg 

44- 

No 

Length 

43- 
42- 
41- 

38- 

Data 

Sitting 

37- 

I 

Height 

36- 
35- 

27- 
26- 

f 

Buttock- 

25- 

T 

Knee  Length 

24- 
23- 

1 

22- 

{ 


1 


NOTE:  o  -  mean 

-  one  standard  deviation  from  mean 
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taken  as  a  whole.   These  effects  reflected  the  presence  of 
unidentified  factors,  not  entirely  explained  by  the  facts 
available.   While  height  increased  steadily  from  1964  to 
1975,  its  components  did  not.   As  seen  in  Figure  3,  which 
shows  mean  values  and  one  standard  deviation  on  either  side 
of  the  mean,  mean  height  increased  from  69.94  to  70.15 
inches,  a  rate  of  0.427  inches  per  decade.   That  rate  con- 
firms the  USAF  prediction  of  0.418  inch  per  decade  t Roebuck, 
et .  al . 3.   However,  the  increased  height  does  not  appear  as 
an  increase  in  the  other  long  bone  measurements  which  are 
components  of  height.   From  1969  to  1975,  both  mean  sitting 
height  and  mean  leg  length  appear  to  have  decreased.   The 
increase  in  height  should  be  traceable  to  an  increase  in 
sitting  height,  an  increase  in  leg  length,  or  both,  but 
neither  was  noted.   Explanations  other  than  erroneous  data 
are  possible,  but  not  supported  by  knowledge  of  the  popula- 
tion.  It  could  be  hypothesized  that  an  increasing 
concentration  of  fat  in  the  buttocks  could  increase  both 
sitting  height  and  leg  length  measurements  while  not 
affecting  height.   This  effect  has  not  been  observed. 

Another  possible  explanation,  operator  induced  loss  of 
precision,  would  be  accompanied  by  increases  in  the  variation 
An  examination  of  Table  V  and  Figure  3  indicate  that  standard 
deviations  remained  relatively  constant,  indicating  the 
absence  of  excessive  random  error  by  equipment  operators. 
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The  inclusion  of  a  constant  error  in  any  of  the  measure- 
ments in  one  or  more  data  sets  could  account  for  the 
apparent  inconsistency.   This  kind  of  error  should  be 
precluded  by  the  use  of  standard  procedures  and  calibrated 
equipment.   The  1969  data  were  collected  from  measurements 
taken  with  a  device  for  which  no  calibration  was  provided. 
There  was  no  calibration  procedure  or  equipment  associated 
with  the  integrated  anthropometric  measuring  device  until 
after  the  introduction  of  the  improved  device  in  1971. 
Although  a  means  of  calibrating  the  improved  anthropometric 
measuring  device  was  provided  after  1971,  the  calibration 
equipment  is  no  longer  available  and  the  date  of  the  last 
calibration  is  unknown. 

The  largest  contributor  to  the  apparent  error  is  the 
diminished  leg  length  in  1975.   The  reduction  of  nearly  an 
inch  is  difficult  to  believe.   Leg  length  is  also  a  diffi- 
cult measurement  to  take  accurately  UMoroney,  19803 .   The 
position  of  the  subject  being  measured  (sitting  erect  with 
his  right  leg  maximally  extended  horizontally  in  front  of 
him)  is  difficult  for  some  individuals  to  attain.   The 
inflexibility  of  some  muscular  legs  prevents  the  leg  from 
being  straightened  completely.   Methods  aimed  at  straighten- 
ing the  leg  can  result  in  the  hips  rolling  forward,  causing 
the  buttocks  to  move  outward  from  the  seat  back.   The 
measured  leg  length  is  then  either  too  long  (with  the  leg 
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straight  and  buttocks  forward)  or  too  short  (with  the  leg 
slightly  bent ) . 

The  1975  data  also  appeared  to  have  been  influenced  by 
rounding.   The  frequency  of  whole  inch  and  half  inch  measure- 
ments is  greater  than  expected.   The  bias  incurred  in  the 
mental  rounding  process  is  indeterminate  and  no  effort  was 
made  to  correct  for  it.   The  existence  of  the  effect  indi- 
cates a  need  for  a  procedural  modification. 

The  exclusion  determined  here  were  based  on  the  consider- 
ation of  only  four  variables.   Although,  of  the  available 
measurements,  they  were  considered  to  be  the  most  important 
with  respect  to  accidents,  injuries,  and  aircraft  controlla- 
bility, CMoroney,  1980]  they  are  not  the  only  measurements 
relevant  to  accommodation.   Administrative  limitations  set 
by  NMPC  established  that  no  more  than  four  variables  would 
be  used.   The  practice  of  collecting  anthropometric  data  at 
the  flight  training  entry  point  and  the  specific  measurements 
taken  there  determined  which  measurements  could  be  used. 
The  inclusion  of  additional  measurements,  such  as  shoulder 
width,  could  only  increase  exclusions.   The  substitutions  of 
other  measurements  could  affect  the  exclusion  results 
similarly. 
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V.   SUMMARY  AND  CONCLUSIONS 

The  Navy  requires  airplane  cockpit  designs  to  accom- 
modate ninety  percent  of  the  Naval  aviation  population.   In 
support  of  that  goal,  designs  are  required  to  accommodate 
the  5th  to  95th  percentiles  (or  the  3rd  to  98th  percentiles) 
of  critical  body  dimensions.   The  success  or  failure  of 
aircraft  designers  in  accommodating  the  population  had  not 
been  quantified,  heretofore.   By  examining  anthropometric 
data  from  three  populations  (1964,  1969,  and  1975),  two  of 
which  were  generated  by  the  Computerized  Accommodated  Percen- 
tage Evaluation  Model  (CAPE),  exclusions  allowed  by  the 
design  specification  and  exclusions  associated  with  each 
aircraft  type  were  determined.   The  numerical  values  corres- 
ponding to  the  percentiles  named  in  the  specification  (5  to 
95  or  3  to  98)  were  compared  to  the  measurements  of  each 
subject  in  each  population  to  determine  the  proportion  of 
subjects  of  each  population  having  at  least  one  measurement 
outside  the  limits  of  accommodation  and  therefor  excluded. 
The  results  were  the  exclusions  expected  from  designing  to 
the  specified  percentiles.   The  aircraft  exclusions  were 
found  by  comparing  aircraft  anthropometric  exclusion  codes, 
assigned  by  NAVAIR,  and  the  measurements  corresponding  to 
those  codes,  to  the  three  populations.   That  process 
resulted  in  an  exclusion  percentage  for  each  aircraft  and 

each  population. 
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It  was  found  that : 

1.  The  design  specification  requiring  accommodation  of  the 
5th  to  95th  percentiles  excluded  26  percent  of  the  specified 
population. 

2.  Designing  to  the  3rd  to  98th  percentiles  excluded  14 
percent  of  the  population. 

3.  Exclusions,  both  those  allowed  by  the  specification  and 
those  peculiar  to  each  aircraft,  varied  over  time. 

4.  All  aircraft  accommodate  more  than  the  required  ninety 
percent  of  the  most  recent  population  sample. 

5.  The  proportion  of  exclusion  caused  by  a  given  variable 
changed  over  time. 

6.  There  are  inconsistencies  in  the  data. 

A.   DESIGN  SPECIFICATION 

The  design  specification  should  result  in  at  least  ninety 
percent  of  the  user  population  being  accommodated  by  each 
aircraft.   Reference  to  the  5th  and  95th  percentiles  and  to 
the  3rd  and  98th  percentiles  of  critical  body  dimensions 
serves  only  to  fog  the  issue.   The  imaginary  '95th  percentile 
man,'  consisting  of  a  95th  percentile  measurement  for  all 
variables,  is  assumed  from  the  requirement  to  design  to  his 
dimensions.   The  tacit  assumption  that  anthropometric 
features  can  be  assembled  and  generalized  to  form  a  95th 
percentile  man  is  false  and  misleading  since  it  assumes  a 
correlation  of  1.00  between  anthropometric  features.   The 
specification  would  be  more  effective  if  it  were  to  require 
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accommodation  of  ninety  percent  of  the  population  and  remain 
silent  on  the  accommodation  of  singular  dimensions.   Deletion 
of  all  reference  to  percentiles  in  the  specification  would  be 
a  productive  change.   However,  anthropometric  data  should 
still  be  reported  in  percentiles  for  use  in  univariate 
applications. 

The  Navy  is  currently  enjoying  the  benefit  of  an  unre- 
quired bonus  in  the  area  of  user  accommodation.   Ninety 
percent  accommodation  was  required,  more  than  that  is 
delivered,  at  least  for  the  four  variables  examined,  based 
on  the  1975  data.   Although  there  appears  to  be  nothing  to 
indicate  that  this  bonus  will  not  continue,  there  is  like- 
wise nothing  to  guaranty  that  it  shall.   If  a  future  aircraft 
design  were  successful  in  accommodating  exactly  the  required 
ninety  percent,  the  consequences  would  amount  to  exclusions 
greater  than  those  being  experienced  on  any  current  Navy 
aircraft  and  the  design  specification  would  still  have  been 
satisfied.   Since  there  is  feeling  in  the  Navy  that  current 
user  accommodation  is  insufficient  and  since  the  ninety 
percent  accommodation  requirement  is  being  met,  either  the 
Navy  must  change  the  ninety  percent  requirement  or  adjust  its 
attitude  toward  current  accommodations.   The  answer  may  be 
found  in  partitioning  the  user  population  into  groups,  each 
of  which  currently  flies  in  a  particular  type  aircraft. 
Then  specify  the  'user  population'  to  be  one  of  those  groups, 
rather  than  all  flying  personnel  in  the  Navy,  and  require 

53 


accommodation  of  the  entire  group.   This  solution,  while 
easy  to  conceive,  would  be  difficult  to  administer. 

B.   COCKPIT  RESTRICTIONS 

The  exclusions  enumerated  in  Table  VIII  indicate  that 
the  Navy  has  a  minimal  problem,  on  the  global  level,  with 
incompatible  planes  and  personnel.   The  attention  that  the 
problem  has  received  would  hardly  seem  justified  by  the 
miniscule  exclusions  indicated  by  the  analysis.   The 
expense  in  terms  of  money,  man-hours,  and  effort  reflected 
in  the  production  of  instructions  by  OPNAV,  BUMED ,  CNATRA, 
and  NAVAIR,  together  with  the  resultant  effort  necessary  to 
comply  with  their  requirements  must  have  been  justified  by 
a  judgement  that  a  problem  existed  but  the  numbers  in  Table 
VIII  do  not  reflect  the  magnitude  of  the  problem.   A  reason 
that  more  impressive  numbers  don't  appear  in  the  table 
might  be  in  the  aircraft  codes'  inaccuracy,  or  in  their 
failure  to  account  for  influential  variables. 

The  accuracy  of  the  aircraft  codes  is  dependent  on  a 
judgement  factor  on  the  part  of  the  people  who  measure  air- 
craft cockpits.   The  anthropometric  codes  assigned  to 
aircraft  are  based  on  cockpit  measurements  adjusted  for 
clothing  and  equipment  worn  by  the  user.   Those  adjustments 
cannot  allow  for  all  variations.   The  codes  assigned  to 
aircraft  can  also  be  affected  by  forecasts  of  the  consequences 
of  the  codes  assigned.   The  prospect  of  excluding  an  unsuit- 
ably large  portion  of  the  user  population  could  have 
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influenced  aircraft  exclusion  code  assignments  in  the  direc- 
tion of  minimal  effect. 

The  high  cost  of  error  does  not  appear  in  the  table  but 
undoubtedly  accounts  for  the  attention  given  to  such  unim- 
pressive exclusion  figures.   The  cost  of  losing  an  airplane 
or  an  aircrew  is  considerable.   That  cost  makes  any  erroneous 
assignment  of  personnel  potentially  disastrous  and  may 
account  for  the  emphasis  in  the  area  of  aircrew  accommoda- 
tion . 

Additionally,  anthropometric  restrictions  can  affect  an 
officer's  career  rather  markedly  by  denying  him  assignment 
to  the  type  aircraft  he  is  motivated  to  fly.   That  fact 
makes  the  globally  small  exclusion  problem  a  crisis  on  the 
individual  level  and  explains  why  the  attention  level 
exceeds  that  expected  from  examining  the  exclusion  figures. 

C.   DATA  INCONSISTENCIES 

The  data  indicate  that  the  values  associated  with  anthro- 
pometric features  shifted  over  time.   Although  the  exclusions 
by  design  specification  in  1964  and  those  in  1975  were  about 
the  same  (see  Table  VIII),  the  individual  variables  respon- 
sible for  the  exclusions  were  of  varying  influence.   In  1964, 
sitting  height  caused  more  exclusions  than  it  did  in  1975, 
while  buttock-knee  length  and  leg  length  caused  more  in  1975 
than  in  1964.   The  1969  data  show  higher  exclusion  than 
either  of  the  other  two  data  sets.   Sitting  height,  buttock- 
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knee  length,  and  leg  length  all  peaked  in  1969,  accounting 
for  the  corresponding  peak  in  exclusions. 

The  exclusion  results  varied  erratically  over  time. 
There  are  inconsistent  effects  in  the  data  and  the  lack  of 
continuity  cannot  be  explained  with  the  information  avail- 
able.  The  leg  length  measurement  is  the  most  likely  source 
of  the  inconsistent  effects  and  the  discontinuance  of  leg 
length  as  a  critical  measurement  might  be  beneficial  if 
another  measurement  could  be  substituted.   Buttock-knee 
length  might  be  used  as  a  substitute  for  leg  length  since 
the  two  variables  are  well  correlated  (r=.8).   Alternatively, 
a  measurement  taken  in  a  sitting  position  approximating  a 
cockpit  seating  position  could  be  substituted.   Measuring 
from  the  heel  to  the  seat  back  via  the  front  edge  of  the 
seat  would  accomplish,  more  directly,  the  purpose  presently 
being  attempted  by  the  leg  length  measurement .   The  problems 
encountered  measuring  leg  length  would  be  eliminated  and  the 
measurements  of  personnel  and  of  aircraft  cockpits  would  be 
more  comparable  and,  consequently,  more  effective  in  matching 
personnel  to  airplanes.   Such  a  measurement  of  'functional 
leg  throw'  would  be  a  special  purpose  dimension  not  included 
in  previous  surveys.   Previous  surveys,  however,  were  not 
faced  with  the  special  problem  now  at  hand  and  the  cost  of 
designing  and  implementing  this  new  measurement  could  be 
recovered  in  the  prevention  of  a  single  aircraft  accident. 
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D.   DATA  LIMITATIONS 

Of  all  the  anthropometric  dimensions  which  might  be  con- 
sidered, the  Navy  collects  data  on  eight  anthropometric 
features  (including  weight)  considered  relevant  to  flight 
safety.   Of  those  eight  variables,  only  four  are  considered 
in  the  process  of  assigning  personnel  to  aircraft .   The 
limitation  to  four  variables  is  not  necessarily  optimum  with 
respect  to  flight  safety  since  the  administratively  driven 
maximum  of  four  variables  was  the  effective  constraint.   The 
inclusion  of  additional  variables,  if  they  were  found  to  be 
relevant,  could  only  decrease  accommodation  further.   The 
decrease  could  be  dramatic  if  the  additional  variables  were 
poorly  correlated  to  those  already  considered. 

The  problems  arising  from  variables  in  combination  (but 
not  visible  in  single  variable  evaluation)  have  been  recog- 
nized [CNATRAINST  13520.1,  1980  3  but  no  description  of  the 
interaction  between  variables  is  included  in  the  Anthropo- 
metric Compatibility  Assignment  Program.   For  example,  the 
minimum  functional  reach  and  leg  length  necessary  to  reach 
controls  varies  as  the  position  of  the  seat  changes.   The 
position  of  the  seat  is  determined  largely  by  the  sitting 
height  of  its  occupant.   An  individual  with  a  short  sitting 
height  must  raise  the  seat  high  enough  to  allow  vision  over 
the  nose  of  the  aircraft.   As  the  seat  rises,  it  also  moves 
aft  slightly,  causing  the  rudder  pedals  to  be  adjusted  aft 
and  also  requiring  a  longer  functional  reach  to  enable  the 
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man  to  reach  the  controls.   Thus,  sitting  height  determines 
in  part  the  effective  minima  for  leg  length  and  functional 
reach  values,  making  required  leg  length  and  functional 
reach  distances  functions  of  sitting  height.   The  values  used 
for  anthropometric  coding,  however,  are  determined  irres- 
pective of  seat  position  changes.   An  additional 
anthropometric  code  should  be  established  to  identify  those 
personnel  who  meet  all  the  limits  but  whose  combination  of 
anthropometric  features  could  create  problems. 

The  procedures  used  to  take  measurements  in  the  1969  and 
1975  data  sets  were  deficient  in  some  respects.   Most  notably, 
the  failure  to  calibrate  the  integrated  anthropometric 
measuring  device  was  unfortunate.   The  apparent  rounding  of 
measurements  also  induced  error.   The  calibration  deficiency 
can  be  corrected  by  instituting  a  requirement  to  document 
periodic  calibration  and  adjustments,  as  is  required  by  the 
Preventive  Maintenance  System.   The  rounding  error  could  be 
eliminated  by  converting  the  measuring  device  to  provide 
digital  displays  of  measurements.   Better,  but  also  more 
expensive,  an  automated  measuring  device  could  record 
measurements  without  requiring  the  operator  to  read  and 
record  numbers.   Such  a  device  would  also  prevent  the  oper- 
ator from  rounding  measurements  or  inserting  any  bias  into 
the  data. 

The  necessity  to  draw  on  outside  sources  for  data  con- 
cerning leg  length,  a  critical  body  dimension  about  which 
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there  is  no  data  in  the  design  specification  (1964  data), 
points  out  a  basic  deficiency  in  the  specification. 
Designers  do  not  consider  leg  length  in  the  design  of 
cockpits.   Instead,  they  use  buttock-knee  length  and  pop- 
liteal height  sitting,  which  are  components  of  leg  length. 
The  assignment  process  would  be  more  consistent  if  similar 
variables  were  used  for  exclusion  codes. 

E.   CONCLUSIONS 

The  data  collection  system  which  supports  the  Anthropo- 
metric Compatibility  Assignment  Program  consists  of 
measurements  taken  on  aircraft  cockpits  as  well  as  the 
equipment  and  procedures  used  to  measure  aviation  personnel. 
Many  inconsistencies  in  personnel  measurements  have  been 
noted  even  though  those  measurements  are  an  essential  ingre- 
dient in  the  successful  operation  of  the  program.   The 
program  will  operate  no  more  effectively  than  its  weakest, 
or  limiting,  factor.   That  limiting  factor  appears  to  be  the 
measurement  of  personnel.   Inconsistencies  in  the  data  appear 
to  be  results  of  measurement  bias.   The  appearance  of  larger 
personnel  in  1969  and  corresponding  higher  exclusions  in 
that  population  are  more  likely  the  changes  of  measurement 
bias  than  any  real  effect  in  the  population. 

Measurement  error  could  be  reduced  by: 

1.  Periodically  calibrating  the  measurement  equipment. 

2.  Randomly  remeasuring  personnel  and  comparing  the  original 
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and  repeated  measurement.   If  significant  discrepancies  are 
noted,  corrective  action  should  be  taken. 

3.   Replacing  the  integrated  anthropometric  measuring  device 
by  one  using  digital  measurement  displays,  or  automated 
measurement  recording. 

The  consequences  of  failing  to  improve  the  Navy's 
personnel  anthropometric  measuring  procedure  will  amount  to 
negating  the  Anthropometric  Compatibility  Assignment  Program 
and  accepting  the  aircraft  accident  rate,  and  ejection  injury 
rate  which  were  sufficiently  undesirable  to  lead  to  the 
program's  establishment.   Additionally,  erroneous  measurement 
can  adversely  impact  on  officers'  career  development  patterns 
and  consequently  on  officer  personnel  retention. 
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APPENDIX  A 
Definitions  of  Measurements 

1.  Buttock-Heel  Length       -see  Leg  Length 

2.  Buttock-Knee  Length       -the  horizontal  distance  from 

the  rearmost  surface  of  the  buttocks  to  the  front  of 
the  kneecaps,  with  the  subject  sitting  erect. 

3.  Buttock-Heel  Length       -see  Leg  Length 

4.  Functional  Reach  -the  distance  from  the  wall  to 

the  tip  of  the  thumb  measured  with  the  subject's 
shoulders  against  the  wall,  his  arm  extended  forward, 
and  his  index  finger  touching  the  tip  of  his  thumb. 

5.  Height  -vertical  distance  from  the 

floor  to  the  top  of  the  head  with  the  subject  standing 
erect  in  bare  feet . 

6.  Knee  Height  Sitting       -the  height  from  the  footrest 

surface  to  the  musculature  just  above  the  knee. 

7.  Leg  Length  -the  distance  from  the  base 

of  the  heel  to  a  wall  against  which  the  subject  sits 
erect  with  his  leg  maximally  extended  forward  along 
the  sitting  surface. 

8.  Sitting  Height  -the  height,  from  the  sitting 

surface,  to  the  top  of  the  head,  with  the  subject 
sitting  erect. 

9.  Stature  -see  height 

10.   Popliteal  Height  Sitting   -the  height  of  the  underside  of 
the  upper  leg  above  the  footrest  surface. 
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APPENDIX  B 

Anthropometric  Codes  as  Percentiles 
of  the  1964  Population 


Variable 

Code 

Measurement 

Equivalent 

Interval 

Percentile 

SITTING   HEIGHT 

9 

40.0-41.0 

99+ 

8 

39.5-39.9 

99+ 

7 

39.0-39.4 

98-99 

6 

38.5-38.9 

96-97 

5 

38.0-38.4 

91-95 

4 

35.0-37.9 

15-90 

3 

34.0-34.9 

3-14 

2 

33.0-33.9 

<3 

1 

32.5-32.9 

<1 

0 

32.0-32.4 

<1 

FUNCTIONAL   REACH 

9 

£27.9 

<1 

8 

28.0-28.4 

<1 

7 

28.5-28.9 

1-2 

6 

29.0-29.4 

3-6 

5 

29.5-30.4 

7-25 

4 

30.5-30.9 

26-39 

3 

31.0-31.4 

40-50 

2 

31.5-32.4 

51-75 

1 

32.5-33.9 

76-94 

0 

£34.0 

>95 

BUTTOCK-KNEE   LENGTH 

9 

>28.0 

>99 

8 

27.0-28.0 

>99 

7 

26.5-26.9 

>98 

6 

26.0-26.4 

97-98 

5 

25.5-25.9 

91-96 

4 

25.0-25.4 

81-90 

3 

24.0-24.9 

50-80 

2 

23.0-23.9 

15-49 

1 

22.0-22.9 

2-14 

0 

&21.9 

<2 

LEG   LENGTH 

9 

49.0-50.0 

>99 

8 

48.0-48.9 

98-99 

7 

47.0-47.9 

95-97 

6 

46.0-46.9 

86-94 

5 

45.0-45.9 

71-85 

4 

43.0-44.9 

34-70 

3 

40.0-42.9 

3-34 

2 

39.0-39.9 

1-2 

1 

38.0-38.9 

<1 

0 

36.0-37.9 

<1 
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Appendix  C 

Program  Listing  for  CAPE  Model 
Used  to  Generate  the  1964  and  1969  Data 


PROGRAM  GENS AMP (FMf AS  .U"E AS . SAMPLE , I NPUT .OUTPUT . 

1  TAPF1=FMEA$,TAPE2=UmEAS.TAPE J= SAMPLE* 

2  TAPcS=INPUT.TAPFft=OUTPUT) 
C 

C  •••  MODULE  A  -  MONTE  CARLO  CUMULATION  ••• 

C       *»•  MCSM  READS  IN  FACTOR  MATRIX  OF  REOUlHEO  ANTHROPOMETRIC 

C       •  »•  MEASUREMENTS  FROM  EITHER  0  I  Sf  FILES  OR  INTERACTIVE 

C       •••  TFOMTNAL.   IT  ALSO  ALLOWS  USFRS  TO  PHMM.  EDIT  AND  SAVE 

C       •••  FILES.   THE  PROGRAM  THEN  COMMUTE  RaNOOm  MULTIVARIATE  NORMAL 

C       •*•  FEATURE  VECTORS.   AT  THE  ENJ  OF  THF  PRO&PAM  IT  CALLS 

C       •••»  SUBROUTINE  TRANS  TO  PFREORm  NECESSARY  DATA  Th  AN^KORmaT  ION 

C       •  •»  OF  ANTHROPOMETRIC  SAMPLES  AND  STORFS  The  RESULTANT 

c        »••  MEASUREMENTS  ON  DISK  FILFS  FOR  CPEWSTATION  EVALUATION 

C        •  »•  MODULE. 

C 

COMMON/mCSm/CORP(78) ,ANTm<2. 12) ,V(25) .NmEAS.KN. 

1  SOROOT  (79) .X  <2S>  . I U.N TmE AS. I  COUNT. 

2  INamE  («»)  .ITmfa; -NlMEAS.NTEMP 
COMMON/MISC/IOIN.IOOUT.ICARO 

C 

C      NMEAS  -   NUMBER  OF  ANTHROPOMETRIC  MEASUREMENTS.  CURRENTLY   12 

C      NSAMP  -  NUMBER  OF  OPERATOR  SAMPLES 

C      DATA  FtLE  DESCRIPTION: 

C      NAME    PROGRAM  REF.     TAPE*   FOR/UN  CONTtNTS 

C     FMFAS      1       FOR        ANT"(2, NMEAS) .CORR(KN) 

C     UMEAS     NTMEAS  2       UN         NMEAS .KN. anTm.COwP 

C     SAMPLE     NTEMP  3       UN         IN<\mE<<»>  »V.x 

C     INPUT 5       — 

C     OUTPUT      — —  6       — 

C 

DATA  IOIN.I00UT.NUCRFW,NTEMP.NTMEAS.NlMEAS/^.6.6.3.4.7/ 

DATA  maxFUNC/*/ 

DATA  MMFAS.NTMFAS.NIMEAS/12.2.^/ 
1 .  TCOUK'T/0/ 

DATA  AMTM/2<.»0/.CORR/79»0/.INAME/<.»10H  / 

DATA  MAXOPT.MAXANTH.MAXACT/3.6.3/ 
C 

C   •••COMPUTE  (KN)  THE  LENGTH  OF  THE  LINEAR  ARRAY  CORR<KN) 
C 

IOOUT=0 

ISAMP=o 
C   ••♦MCSM  MAIN  MENU 
100    CONTINUE 

I0IN=5 

WRITE  (IS. 290) 

CALL  INPT  ( 1 ,1 .NVALS7. NMEAS. XVAL. IEOF) 

GO  TO  (130. 100. 110. 120. 2*0. 120)  .  1EOF-1 
110    CONTINUE 

CALL  HELP  (1) 

GO  TO  100 
120    CONTINUE 

WPITC  (*.3<»0) 

GO  TO  100 
130    CONTINUE 

IF  (MMFAS.LT. 1 .OR. NMEAS. GT.MMFAS)  GO  TO  120 

KN=NMFAS» (NMEAS»l)/2 

IINPT=0 
C 

C  •••  GENERATE  OPERATOR  SAMPLE  FROM  MFANS.STD  UEV  h.    CORR  MATRIX 
C 
1*0    CONTINUE 

IOIN=S 
C   •••ANTHROPOMETRIC  FUNCTIONS 
lbO    CONTINUE 

WRITF  «A.300> 
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CALL  INPT  (1 ,1 .NVALSZ.IOPT.XVAL.IEOF) 

GO  TO  (180,140.160,170,280,170),  IEOF»l 
160    CONTINUE 

CALL  HELP  (2) 

GO  TO  ISO 
170    CONTINUE" 

WPITE  (6.340) 

GO  TO  ISO 
180    CONTINUE 

IF  (IOPT.LT.l.OR.IOPT.GT.MAXANTH)  GO  TO  170 
C  INPUT, EDIT, SAVE, P^INT.GENEKATF  OPER  SAME. END 

GO  TO  (190, 200. 220.2*0. 260, 280) .  IOPT 
190    CONTINUE 
C  »*«inpuT  ANTHROPOMETRIC  MEASUREMENTS 

CALL  ACCPT  (ISAMP) 

IINPT=1 

GO  TO  140 
200    CONTINUE 
C   »»«EOIT  ANTHROPOMETRIC  DATA  FILE 

IF  (IINPT.EQ.1)  GO  TO  210 

WRITE  (6.350) 

GO  TO  140 
210    CONTINUE 

CALL  EDIT 

GO  TO  140 
2?0    CONTINUE 
C   *»»SAVE  ANTHROPOMETRIC  OATA  FILE  ON  UMEAS 

IF  (ISamP.EO.1)  GO  TO  230 

WDITF  (6,320) 

GO  TO  140 
230    CONTINUE 

CALL  SAVE 

GO  TO  140 
240    CONTINUE 
C   «»«poiNT  ANTHROPOMETRIC  DATA 

IF  (IlMPT.EQ.l)  GO  TO  250 

WRITE  (6,360) 

GO  TO  1*0 
250    CONTINUE 

CALL  PRINT 

GO  TO  140 
C 

C  **»  GENERATE  OPERATOR  SAMPLE  FROM  ANTH  MEAS  DATA 
C 
260    CONTINUE 

IF  (ISAmP.se. I)  GO  TO  270 

CALL  OPSAMP  (ISAMP) 

GO  TO  140 
270    CONTINUE 

WRITE  (6.330) 

GO  TO  140 
C   «»*TFPMINATION  OF  MCSM 
280    CONTINUE 

WRITE  (6.310) 

STOP 
C 

290    FORMAT  (•'  ENTER  NUMBER  OF  MEASUREMENTS  (1-12)? — ") 
300    FORMAT  (»   GENERATE  FUNCTIONS ( 1-INPUT ;2-EDIT; 3-SAVE; 4-PPINT;". 

♦  "5-GENERATE:6-ENn) — ") 

310    FORMAT  ('•  END  OF  OPERATOR  SAMPLE") 

320    FORMAT  ('•  UNABLE  TO  SAVE  MEASUREMENT  DATA  -  FATAL  ERRORS") 

330    FORMAT  (••  UNARLE  TO  GENERATE  SAMPLE  -  FATAL  ERRORS") 

340    FORMAT  ("  INVALID  RESPONSE") 

350    FORMAT  ("  ANTHROPOMETRIC  DATA  MUST  BE  INPUT  BEF()Rt  TnE  EDIT  ". 

♦  "FUNCTION  CAN  BE  USED") 

-\*()  FORMAT.*'?  ANTHBOPOMFTPTr  nATA  MUST  RF  TNPUT  HFFORE  ThF  PRINT". 
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♦   "FUNCTION  CAN  3E  USED") 

END 

SUBROUTINE  INPT  ( ITYPE , I VALSZ.NVALSZ. I VAL .XVAL, IEOF ) 

COMMON/MISC/IOIN.IOOUT.ICAPD 
C      FPEE  FIELD  INPUT  ROUTINE  —  READS  IN  INPUT  FROM  THE  kEYmOARD  IN 
C  FREE  FIELD  INPUT 

C  IMPUT 

C  ITYPE  -  SPECIFIES  THE  TYPE  OF  DATA  TU  BE  ENTERED 

C  1=INTEGER   2=REAL   3=ALPHA\UMERIC 

C  .  IVALSZ  -  THE  NUMBER  OF  ITEMS  EXPECTED  TO  BE  ENTEREO 

C  OUTPUT 

C  NVALSZ  -  THE  NUMBER  OF  ITEMS  ENTERED  BY  THE  USER 

C  IVAL  -  APRAY  CONATIMING  INTEGER  OR 

C  ALPHANUMERIC  DATA  SU°PLIED  BY  THE  USER 

C  XVAL  -  APRAY  CONTAINING  REAL  DATA  SUHPLIED  Br  THE 

C  USER 

C  IEOF  -  END  OF  FILE  FLAG 

C  1=IF  BLANKS  Ok  a  CARRIAGE  CONTROL  IS  THE 

C  ONLY  INPUT 

C  3=  IF  A  QUESTION  MARK  IS  ENTERED 

C  3=IF  INPUT  IS  INVALID 

C  4=IF  A  $  HAS  BEEN  ENTERED 

DIMENSION  I  CHAR (80) .  IVAL ( IVALSZ) .XVAL < IVALSZ) 

ISTOP=0 

IF  (IVALSZ. GE.O)  GO  TO  100 

ISTOP=l 

IVALSZ=IABS(IVALSZ) 
100    CONTINUE 

IOUOTE=0 

NVALSZ=0 

IFOF=0 

IF  (IOOUT.EQ.l)  GO  TO  120 

DO  110  1=1 ,80 

ICHAR(I)=lH 
110    CONTINUE 
C      READ  IN  THE  INPUT 

READ  <IOIN,330)  ICHAR 

IF  (I0IN.NE.5)  WRITE  (6.330)  ICHAR 

IF  (POF(IOIN) .NE.l)  GO  TO  130 

TFOF=l 

RETURN 
120    CONTINUE 

I1=ICAPD 

GO  TO  140 
130    CONTINUE 
C      INITIALIZE  THE  CARD  COLUMN  COUNTERS 

11=1 
140    CONTINUF 

I2  =  fl0 
ISO    CONTINIJF 

IF  (ISTOP.EO. LAND. IVALSZ. EO. NVALSZ)  RETURN 
C      SEARCH  FOR  THE  NEXT  WORD 

IF  (I1.LE.I2)  GO  TO  160 
C      SEARCH  FOR  THE  FIRST  NON-BLANK  CHARACTER  (START  OF  WORD) 
C      CHECK  IF  EXCESS  HAS  BEEN  ENTERED 

IF  (NVALSZ. LE. IVALSZ)  RETURN 

TE0F=5 

RETURN 
160    CONTINUE 

DO  170  1=11,12 

IF  (ICHflR(I) .E0.1H  )  GO  TO  170 

IFI9ST=I 

GO  TO  ISO 
170    CONTINUE 

IF  (NVALSZ. EO.O)  IEOF=l 
C      CHECK  IF  EXCESS  HAS  BEEN  ENTERED 
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IF  (NVALSZ.GT.IVALSZ)  IE0F=5 

RETURN 
IPO    CONTINUE 
C      SFARCH  FOR  LAST  NON-BLANK  CHAPACTER  (END  OF  wORD) 

IF  (ICHAP(IFIRST) .EO.lH?)  IEOF=2 

IF  (ICHAR(IFIRST) .EQ.lHS)  IEOF=4 

IF  (IFOF.NE.O)  PFTURN 

IF  (ITYDE.EQ.^)  GO  TO  300 

IF  (ITYPE.NE.3)  GO  TO  190 

IF  (ICHAP(IFIPST) .NE.lH».AND.ICHAH(IFIRST) ,NE.lH«>  GO  TO  190 

IFIRST=IFIRST*1 

IOUOTE=l 
190    CONTINUF 

DO  220  I=IFIRST.I2 

IF  (IQUOTE.EQ.O)  GO  TO  200 

IF  (ICHAR(I) .FQ.lH".OP.ICHAR(I> .EQ.lHM  GO  TO  210 

GO  TO  220 
200    CONTINUE 

IF  (ICHftR(I) .NE.1H  .ANO.ICHAR(I) .NE.iH, )  GO  TO  220 
210    CONTINUF 

ILAST=I-1 

GO  TO  230 
220    CONTINUE 

ILAST=I? 
230    CONTINUE 

NVALSZ=MVALSZ»1 
C      RESET  THE  FIRST  CHARACTER  TO  HE  SEARCHED  FOR  THE  '^EXT  WORD 
C      INCREMENT  THE  NUMBER  OF  WORDS  FOUND 
*I1*ILAST*2 

IF  (NVALSZ.GT.IVALSZ)  GO  TO  150 
C      CHECK  IF  INPUT  IS  ALPHA-NUMERIC 

IF  (ITYPE.EQ.3)  GO  TO  310 
C      INPUT  TS  NUMERIC 
C      CHECK  IF  NUMBER  IS  SIGNED 

ISIGNsl 

IF     UCHAR(IFIRST) .NE.lH-.AND.ICHAft(IFlRST) .NE.lH*)    GO    TO    240 

IF     (ICHAR(IFIRST) .E0.1H-)     ISIGN=-1 

IFIPST=IFIRST*1 
2<+0         CONTINUE 

INUM=0 

mEc=o 

IDECF=0 

DO  260  I=IFIRST»ILAST 

IF  (ICmaR(I) .NE.1H.)  GO  TO  250 

IDECF=IOECF*l 

IDEC=1 

GO  TO  260 
250    CONTINUF 

CALL  ICONV  (ICHAR(I) .IADD»IE0F) 

IF  (IEOF.EQ.3)  RETURN 

INUM=INUM»10*IADD 

IF  (IDEC.NE.O)  IDEC=IDEC*1 
260    CONTINUF 

IF  (IDFCF.LT.2)  GO  TO  280 
270    CONTINUF 

IFOF=3 

RETURN 
280    CONTINUE 

INUM=INUM»ISIGN 

XNUM=FLOAT(INUM) 

IF  ( (ITYPE.EQ.l.AND.IDEC.LE.l).OR.(ITYPE.EQ.2))  GO  TO  290 

IEOF=3 

RETURN 
290    CONTINUF 

IF  (IDEC.NE.O)  XNUM=XNUM/10.0*»(IDEC-1) 

TF     UTYPF.F0.1)      TVAI    (N\/AI  S7)  rTFTX  (XNIIM) 
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IF  ( 1  I YHt.fc.U.3)  AVAL (NVALbZ) =XNUM 

ICAPD=ILAST*1 

GO  TO  150 
300    CONTINUE 

ILAST=TVALSZ»10 
310    CONTINUF 
C      INPUT  IS  ALPHANUMERIC 

DO  320  J=1,IVALSZ 

IVAL(J)=10H 
330    CONTINUF 

L=ILAST-IFIPST*1 

NVALSZ=(L-1)/10*1 

IF  (L.GT.10»IVALSZ)  GO  TO  270 

ENCODE  (L»330.IVAL<1>  )  { ICHAR ( I )  » I=IFIRST, ILAST) 

ICARD=I1 

RETURN 
C 
330    FORMAT  (80A1) 

END 

SUBROUTINE  ICONV  ( JCHAR , 1 1 , IEOF) 
C      THIS  FUNCTION  CONVERTS  ALPHA  NUMBERS  TO  NUMBER  NUMBEkS 

OIMEMSION  JNUM(IO) 

DATA  JNUM/1H0.1H1 . 1H2 » 1H3 ♦ 1H4 , 1 H5 , 1H6 » 1H7 * 1H8 » 1H9/ 

DO  100  J=l,10 

IF  (JCHAP.NE.JNU^(J) )  GO  TO  100 

I I=J-1 

RETURN 
100    CONTINUE 
C      INVALID  INPUT  NOT  A  NUMBER 

IEOF=3 

RETURN 

END 

SUBROUTINE  ACCPT  (ISAMP) 
C 

C  »»•  INPUT  ANTHROPOMETRIC  DATA 
C  •••  MEANS,  STD  DEV,  AND  CORRELATION  MATRIX 
C 

C0MM0N/mCSM/C0RR(7B) ,ANTM(2»12) *V(25) *NMEAS»KN» 
1     SOROOT(78)  *X(25) ♦ IU,NTmEAS» I  COUNT* 
?     INAME(^) »ITMEAS.NIMEAS»NTEMP 

COMMON/MISC/IOIN.IOOUT.ICAPD 

DIMENSION  3UFF(3) *TEMP(12) 

DATA  BUFF/3*0/ 

ISAMP=0 

ICOUNT=0 

INLINE=0 

I0IN=5 
C 

C  •♦•READ  MODE  -  1=INTERACTIVE  2=FILE 
C 
100    CONTINUE 

WPITE  (6,610) 

CALL  INPT  (1,1 tNOINPT.IMODE.XINPT. IEOF) 

GOTO  (130,110.120.110*490,110),  IEOF*l 
C  INVALID  RESPONSE 
110    WRITE  (6.650) 

GO  TO  100 
C  HELP 
120    CALL  HELP  <<►) 

GO  TO  100 
130    IF  (IMODE.EQ.l)  GO  TO  200 

IF  (IM0DE.NE.2)  GO  TO  110 
C 
C  •••  FILE  TYPE:  l=FORMATTED   2=UNFQRMATTE0 

c 

1<*0    CONTTNIIF 
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WPITE  (6.620) 

CALL  INPT  (1  ,1,NOINPT.IFFORM,XINPT»IEOF) 

GO  TO  (170,150,160,150,490,150),  IEOF*l 
150    WRITE  (6.650) 

GO  TO  140 
160    CALL  HELP  (5) 

GO  TO  140 
170    IF  (IFFOPM.EO.l)  GO  TO  190 

IF  (IFFORM.NE.2)  GO  TO  150 
C 
C      ***  UNFORMATTED  READ 

REWIND  NTM£AS 
C   *»»R(tao  THE  NUMBER  OF  ANTHROPOMETRIC  MEASUREMENTS. 
C   «**CORRELATION  COMPONENTS 

READ  (NTMEAS)  NMEAS.KN. INAME 

IF  (EOF(NTMEAS) .NE.l)  GO  TO  180 

WPITE  (6.670) 

GO  TO  490 
C   »*»RPAD  THF  mean  AND  STANDARD  DEVIATION  VALUES 
180    RFAD  (NTMEAS)  ( < ANTM ( 10. JO) . 10=1 .2) . JO=l .NMEAS) 
C   *«»RPAD  THE  CORRELATION  MATRIX  VALUES 

READ  (MTMEAS)  (COPR ( 10) , 10=1 ,KN) 

REWIND  NTMEAS 

WPITE  (6.690) 

ISAMP=1 

RFTU»N 
C  **•  SET  INPUT  LOGICAL  UNIT  TO  1  FOR  FILE  FMEAS 
190    IOIN=l 
200    CONTINUE 
C 

C   *»«ACCEPT  OPERATOR  SAMPLE  FILE  NAME 
C 

CALL  NAMEFL  (IERROR) 

ICOUMT=ICOUNT*IERROR 

INLINE=INLINE*1 
C      *•»  INPUTS  MEANS  AND  STANOARD  DEVIATIONS  INTO  ARRAY  ANTm 

IF  (IMOPE.E0.1)  WRITE  (6.500) 

IF  (IMOOE.NE.l)  WPITE  (6.510)  INAME 

MAXINPT=2 

DO  320  IPOW=l. NMEAS 
210    CONTINUE 

IF  (.IMODE.EO.l)  WRITE  (6.530)  IROW 
220    NAMT=0 
230    CONTINUE 

IT0TIN=«AXINPT-NAMT 

ICOL=NAMT*l 

CALL  INPT  (2.IT0TIN.N0INPT,IVAL.ANTM(IC0L.IP0W) ,ItOF) 

IF  (IFOF.EQ.5.0P.IEOF.EO.0)  NAMT=NAMT *NOINPT 

IF  (IMonF.EO.l)  GO  TO  240 

INLINE=INLINE*1 

IF  (NAMT.EQ.MAXINPT)  GO  TO  290 

IF  (IEOF.E0.1)  GO  TO  480 

WPITE  (6.540)  IPOW 

IF  (IEnF.E0.5)  GO  TO  270 

WRITE  (*.700)  INLINE 

IC0UNT=IC0UNT*1 

GO  TO  3?0 
C*«*  INTERACTIVE  IEOF  PROCESSING 
2U0         CONTINUF 

GOTO  (280.210,250.260.490.270).  IEOF*l 
250    CONTINUE 

CALL  HELP  (6) 

GO  TO  210 
260    CONTINUE 

WPITE  (6.650) 

GO. TO  210 
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270    CONTINUE 

IOLIST  =  NOI\|PT-ITOTIN 

W°ITE  (6.600)  IOLIST 

IF  (IMODE.NE.l)  ICOUNT=ICOUNT»l 

GO  TO  290 
2b0         CONTINUE 

IF  .  (NAMT.GE.2)  GO  TO  290 

WPITE  (6.550) 

GO  TO  230 
290    CONTINUE 

IF  (NAMT.GT.2)  NAMf=2 

JRERP=0 

DO  300  ICOL=l.NAMT 

CALL  SEASON  ( 1  .  ICOL  »  IROW t ANTM ( 1 . IROW) . ANTM (2. IROW) * IFATAL) 

IRERR=IRERR»IFATAL 
300    CONTINUE 

IF  (IRERR.EO.O)  GO  TO  320 

IF  (IMODE.NE.l)  GO  TO  310 

WPITE  (6.630)  IROW 

GO  TO  ??0 
310    ICOUNT=ICOUNT*IRERR 
320    CONTINUE 
C 

C      •••  ANTHROPOMETRIC  CORRELATION  MATRIX  VALUES  ARE  mCCERTED 
C      *•*  AND  PROCESSED. 
C 

IF  (IMOPE.EO.l)  WRITE  (6.560) 

IF  (IMODE.NE.l)  W»ITE  (6. 520) 

DO  460  IR0W=1»NMEAS 
330    CONTINUE 

NAMT=0 

MAXINPT=NMEAS*l-IROW 
340    CONTINUE 

IF  (iMonE.EQ. LAND. IROW. EO.NMEAS)  WRITE  (6.580)  1*0-  .MAXINPT  .NME^S 

IF  (IMOOE.EQ. LAND. IROW. NE.NMEAS)  WRITE  (6.570)  I *Ow .MAXINPT . IPO*. 
1NMEAS 
350    CONTINUE 

IT0TIN=M4XINPT-NAMT 

ISTART=IROW»NAMT 

CALL  INPT  (2»IT0TIN.N0INPT,IVAL»TEMP(ISTART) , IEOF) 

IF  (IEOF.EQ.0.OR.IEOF.EQ.5)  NAMT=NAMT*NOINPT 

IF  (IMOnE.EQ.l)  GO  TO  360 

INLINE=INLINF»1 

IF  (MAMT.EQ. MAXINPT)  GO  TO  410 

IF  (IEOF.EQ.O)  GO  TO  350 

IF  (IEOF.EO.l)  WRITE  (6.710) 

IF  (IEOF.EO.l)  GO  TO  470 

IF  (IEOF.E0.5)  GO  TO  390 

ICOUNT=ICOUNT*l 

WPITF  (6*700)  INLINE 

GO  TO  4*0 
C  »•♦  INTERACTIVE  INPUT  ERROR  PROCESSING 
360    CONTINUE 

GO  TO  (400,380.370.380.490.390).  IEOF*l 
370    CONTINUE 

CALL  HELP  (7) 

GO  TO  340 
380    CONTINUE 

WPITE  (6.650) 

GO  TO  340 
390    CONTINUE 

IOLIST=NOINPT-ITOTIN 

WRITE  (6.600)  IOLIST 

IF  (IMODE.NE.l)  ICOUNT=ICOUNT*l 

NAMT=MAXINPT 

GO  TO  410 
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400    CONTINUE 

IF  (NOTNPT.EQ.ITOTIN)  GO  TO  410 

ITOTIN=MAXINPT-NAMT 

WPITE  (6*590)  ITOTIN.IPOW 

GO  TO  350 
410    CONTINUE 

ILAST=IROW*NAMT-l 

IRERR=0 

00  420  ICO|_  =  IROW.lLAST 

CALL  PEASON  (2.IC0L»IR0W.TEMP(IC0L) .O.IFATAL) 

IRERR=IPEPP»IEATAL 
420    CONTINUE 

IE  (IRERR.EO.O)  GO  TO  440 

IE  (IMODE.EQ.2)  GO  TO  430 

WRITE  (6»640)  IPOW 

GO  TO  3^0 
430    ICOUNT=ICOUNT*IRERR 

GO  TO  4*0 
440    CONTINUE 

00  450  TCOL=IROW»NMEAS 

N=N0X(IR0W,IC0L) 

CORR(N)=TEMP(TCOL) 
450    CONTINUE 
4*0    CONTINUE 
470    CONTINUE 

IE  (ICOUNT.EQ.O)  ISAMP=1 

IE  (IMOOE.NE.1)  WPITE  (6*680)  ICOUNT 

GO  TO  490 
480    WPITE  (6.660) 

RETURN 
490    CONTINUE 

RETURN 
C 

500    EORMAT  (••  ENTER  PAIRS  OE  MEANS  AND  STANOARD  DEVIATIONS  — ") 
510    EORMAT  (//"  DIAGNOSTICS  OE  FORMATTED  HEAD  '•/••  OPE*  DESC:  "»4A10/> 
b20    FORMAT  (//'«  CORRELATION/ANTHROPO"ETRIC  MEASUREMENTS:") 

530    FORMAT  ('•  MEASUREMENT  -  ".12."  ? ") 

5*0    FORMAT  ("  MEASUREMENT  -".12) 
550    FORMAT  ('•  ENTER  STANDARD  DEVIATION  — ") 

560    FORMAT  ('•  ENTER  ANTHROPOMETRIC  CORRELATION  MATRIX  OATA  IN  A  ROWwIS 
IE  UPPER  '•/"  TRIANGULAR  FORM."/"  NOTE:  ALL  DIAGONAL  ELEMENTS  MUST"* 

♦  "  EQUAL  1.000") 

570    FORMAT  (5H  POW  .12. "(".12."  ELEMENTS.  COLUMNS  ". 12."-". 12.") ") 

ShO    FORMAT  (5H  POW  .12. "(".12."  ELEMENT.  COLUMN  ".12.")") 

5S>0    FORMAT  <2lH  ENTER  THE  REMAINING  .I2.1HH  ELEMENTS  f  OR  HOW  .12) 

600    FORMAT  ("   THE  LAST  ".I?."  INPUTS  HAVE  BEEN  IGNORED") 

610    FORMAT  ("  INPUT  MODE ( 1=INTERACTI Vt ;2=FILE ) — ") 

620    FORMAT  ("  READ  MODE < 1 =FORMATTED: 2=UNF0RmaTTED> — "> 

630    FORMAT  ("   RE-ENTER  OATA  FOR  ".12) 

640    FORMAT  ("  RE-ENTER  DATA  FOR  ROW  ".13) 

650    FORMAT  ("  INVALID  RESPONSE") 

660    FORMAT  ("  UNEXPECTED  EOF  ON  FMEAS  FOR  MEAS  ".12) 

670    FORMAT  ("  FILE  UMEAS  IS  NOT  ATTACHED  TO  ThE  PROGRAM"/ 

♦  "  UNFORMATTED  READ  ON  THIS  FILE  IS  ABORTED  ".//) 

6M0    FORMAT  ("  TOTAL  ERRORS  ON  INPUT  JATA  ".15) 

690    FORMAT  ("   MEAN.STO  DFv  *  CORRELATION  DATA  READ") 
700    FORMAT  ("  INVALID  DATA  AT  LINE  "»I5> 
710    FORMAT  ('•  UNEXPECTED  END  OF  FILE") 

END 

SUbPOUTINE  EDIT 

COMM0N/MCSM/C0RP(78) .ANTM(2.12) .V(25) .NMEAS.KN. 

1  SQROOT (78) .X (25) . IU.NTmE AS. ICOUNT. 

2  INAME(4) .ITMEAS.NIMEAS.NTEMP 
COMMON/MISC/IOIN.IOOUT.ICARD 
DIMENSION  9UFF(3) 

OATA  (Jl l*T'T/.1*.0 / 
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c 

C      o*»  MODIFICATIONS  ON  MEANS , STANDARD  DEVIATIONS* 

C      *••  AND  ANTHROPOMETRIC  CORRELATIONS 

C 

I0IN=5 
100    CONTINUE 

WRITE  (6.350) 

CALL  INPT  (1  ,1,NVA|_SZ,I0PT,XVAL»IE0F) 

GO  TO  (130,100,110,120,340,120) ,  IEOF»l 
110    CONTINUE 

CALL  HELP  (8) 

GO  TO  100 
120    CONTINUE 

WRITE  (6*360) 

GO  TO  100 
130    CONTINUF 

IF  (I0PT.LT.1.0R.I0PT.GT.5)  GO  TO  120 

GO  TO  (140,140,2<»0,330,34Q)  ,  IOPT 
C 

C      **»  ICOL=l  FOR  MEAN  PROCESSING 
C      »»«  ICOL  =  2  FC-^  STANDARD  DEVIATION  PROCESSING 
C 
1«»0    CONTINUE 

ICOL=IOPT 
150    CONTINUF 

MAXINPT=2 

NAMT=0 
160    CONTINUF 

IF  (ICOL.EQ.l)  WRITE  (6,380) 

IF  (ICOL.E0.2)  WRITE  (6,390) 
170    CONTINUF 

IVALSZ=MAXINPT-NAMT 

J=NAMT*1 

CALL    INPT     (2,IVALSZ»NVALSZ,IVAL»8UFF( J) ,IEOF> 

GOTO    (220, lflO, 200, 190, 100, 210) ,     IEOF-1 
180         CONTINUE 

WPITE     (6,370) 

GO  TO  1*0 
190    CONTINUE 

WPITE  (6,360) 

GO  TO  160 
200    CONTINUF 

IF  (ICOL.EQ.l)  CALL  HELP  (11) 

IF  (ICOL.E0.2)  CALL  HELP  (12) 

GO  TO  160 
210    CONTINUE 

IOLlST=NVALSZ-IVALSZ 

WPITE  (6,400)  IOLIST 
220    CONTINUE 

IROW=RMFF(l) 

NAMT=NVALSZ»NAMT 

IF  (NAMT.GE.MAXINPT)  GO  TO  230 

IF  (ICOL.EQ.l)  WRITE  (6,420)  IROW 

IF  (ICOI..EQ.2)  WRITE  (6,430)  IROW 

GO  TO  170 
23P    CONTINUE 

CALL  REASON  ( 1 ,  ICOL  ♦  IROW ,8UFF ( 2) , 0 , IFATAL) 

IF  (IFATAL. EQ.l)  GO  TO  150 

ANTM(ICOL»IROW)=BUFF(2) 

GO  TO  100 
C 

c      ••»  PROCESSING  FOR  ANTHROPOMETRIC  CORRELATIONS 
C 
240    CONTINUE 

MAXINPT=3 

NAMT=0 
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2S0 
260 

270 

2*0 
290 
300 

310 


320 


C 
C 
C 

330 


CONTINUE 

WPITE  (6.410) 

CONTINUE 

IVALSZ=MAXINPT«-NAMT 

J=NAMT*] 

CALL  INPT  (2,IVALS7»NVALSZ»IVAL.3UFF(J) ,IEOF) 

GO  TO  (310,270*290,280, 100,300) ♦  IEOF*l 

CONTINUE 

WPITE  (ft, 370) 

GO  TO  250 

CONTINUE 

WRITF  (6.350) 

GO  TO  250 

CONTINUE 

CALL  HELP  (13) 

GO  TO  250 

CONTINUF 

IOLIST=NVALSZ-IVAL5Z 

WPITF  (6*400)  IOLIST 

GO  TO  320 

CONTINUE 

ICOL=BUFF(i) 

IPOW=BUFF(2) 

NAMT=NAmT*NVALSZ 

IF  (NAMT.GE.MAXINPT) 

IF  (NAMT.EO.l)  WRITE 

IF  (NAMT.EQ.2)  WRITE 

GO  TO  ?60 

CONTINUE 

CALL  REASON  ( 3, ICOL * IROW ,BUFF < 3) . 0 , IFATAL) 

IF  (IFATAL.EQ.l)  GO  TO  240 

CALL  PFASON  (2,ICOL.IPOW.6UFF(3) »0,IFATAL> 

IF  (TFATAL.EO.l )  GO  TO  240 
•••LOAD  THE  VALUE  OF  BUFFO)  INTO  T-»E  APPROPRIATE 
•••POSITION  WITHIN  THE  CORR(N)  ARRAY 

N=NOX(IROW,ICOL> 

CORR(M)=BUFF(3) 

GO  TO  100 

•  ••  PRINT  OPTION 


GO  TO  320 
(6,440)  ICOL 
(6,450)  ICOLtlROW 


C 

340 

C 

350 
360 
370 
380 
390 
400 
410 
4^0 
430 
44  0 
450 


CONTINUE 

CALL  POINT 

GO  TO  100 

•♦♦  ENO  OF  EDIT 

CONTINUE 

RETURN 


MODE 


FORMAT  (••  EDIT  M00E(1-mEANS*2-STl):3-ANTH.C0RR;4-PkinT;5-END> — ") 

FORMAT  ("  INVALID  RESPONSE") 

FORMAT  ("  NO  INFORMATION  ENCOUNTERED") 

FORMAT  ("  INPUT  INDEX  AND  MEAN  VALUE") 

FORMAT  ("  INPUT  INDEX  AND  STANDARD  DEVIATION") 

FORMAT  ('•  THE  LAST", 12,"  INPUTS  HAVE  BEEN  IGNORED") 

FORMAT  (••  INPUT    COLUMN  ,  ROW    AND  CORRELATION") 

FORMAT  {"  INPUT  MEAN  VALUE  FOR  INDEX  ",I3) 

FORMAT  ("  INPUT  STANDARD  DEVIATION  FOR  INDEX  "»I3> 

FORMAT  ("  INPUT  ROW  AND  CORRELATION  FOR  COLUMN  »,I2) 

FORMAT  ("  INPUT  CORRELATION  FOR  ROW  "♦I2»"  COLUMN  ",I2> 

END 

SUBROUTINE  NAMEFL  (IEPROR) 

COMM0N/MCSM/CORR(78) ,ANTM(2,12) »V(25) .NMEAS.KN* 

1  SOPOOT (7fl) ,X(25) »IU,NTMEAS,ICOUNT, 

2  INAME (4) ,ITMEAS»NIMEAS,NTEMP 
COMMON/MI SC/IO I N,IOOUT,ICAPD 


72 


C   »**PROMPTS  USER  FOR  NAMF  OF  SAMPLE  FILE 
C 

DIMENSION  IRUFF(4) ,XVAL<4) 

DATA  IRUFF/4»10H  / 

IFRROR=0 
C 

C   »««INITIALIZE  CONTENTS  OF  IBUFF<4)  WITH 
C   »»#RLANK  CHARACTERS 
C 

DO  100  J=l,4 

IRUFF(J)=10H 
100    CONTINUE 
C 
110    CONTINUE 

IF  (I0IN.E0.5)  WRITE  (6,180) 

CALL  INPT  <4,4,NVALSZ,IBUFF,XVALf IEOF) 

IF  (IEOF.LT.2.0R.IFOF.EQ.5)  GO  TO  150 

IF  (I0IN.NE.5)  GO  TO  140 

GO  TO  (120,130,170),  IFOF-1 
1?0    CONTINUE 

CALL  HFLP  (1^' 

GO  TO  110 
130    CONTINUE 

WRITE  (6,190) 

GO  TO  110 
140    CONTINUE 

WRITE  (6,190) 

IFRROR=l 
C   »*»LOAD  CONTENTS  OF  IBUFF  INTO  INAME  ARRAY 
ISO    CONTINUE 

DO  160  J=l,4 

INAMF (J)=I3UFF(J) 
160    CONTINUE 
170    RFTURN 
C 

1«0    FORMAT  ("  OPERATOR  SAMPLE  DESCRIPTION  (40  CHAk  MAX) ••) 

190    FORMAT  ('•  INVALID  RESPONSE") 

END 

SUBROUTINE  OPSAMP  (ISAmP) 
C   »»»GENEPATE  OPERATOR  SAMPLE 

COMMON/MCSM/CORP(78) ,ANTM(2,12) »V(25) ,NMFAS,KN, 
1     SOPOOT (78) ,X(25)  , IU,NTMEAS, I  COUNT, 
?     INAME(4)  ,ITMEAS,NIMEAS,NTEMP 

COMMON/MISC/IOIN,IOOUT,ICARD 

DIMENSION  IRAN(9) 

OAT A   (IRAN(I) ,1  =  1. 9 )/33 333333, 56555555, 7777777 7 ,22222222  If 
1  444444443, 666666669, 88b888b89, 99999999sr,  l  ^3456789/ 

C 

C   »o«LINEAP  ARRAY  C0RR(78)  IS  LOADED  INTO  LINEAR 
C   oooniJMMY  4PRAY  SQP00T(78)  BEFORE  MATRIX  TRANSFOR- 
C   »»*MATIONS.   ALL  MATRIX  TRANSFORMATIONS  ARE 
c   »»*pfrf0RmED  ON  SQR00T(78)  TO  PRESERVE  C0RR(78) 
C   »*«AS  IMPLEMENTED  WITHIN  THE  EDITING  PROCESS. 
C 

I0IN=5 

00  100  K=1,KN 

SOROOT(K)=CORR(K) 
100    CONTINUE 
c      *•»  pfpfoRM  MATRIX  TRANSFORMATIONS 

CALL  SOPMX  (NERROR) 

IF  (NF.RROR.EO.l)  RETURN 
C      ••*  ACCEPT  OPERATOR  SAMPLE  SIZE 
110    CONTINUE 

WPITE  (A. 320) 

CALL    INPT     (1.1 ,NVALSZ,NSAMP,XVALf IEOF) 

GO    TO     (150. 1?0. 130. 140. 290. 140) •     IF0F*1 
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120    CONTINUE 

WPITE  (6.330) 

GO  TO  110 
130    CONTINUE 

CALL  HELP  (15) 

GO  TO  110 
140    CONTINUE 

WPITE  (6.340) 

GO  TO  110 
150    CONTINUE 

IE  (NSAmP.GT.O)  GO  TO  160 

W»m  (6.350) 

GO  TO  HO 
C 

C   »»*A.SK  USEP  EOP  MACHINE  TYPE  (32  OR  60  BIT) 
C 
160    CONTINUE 

WPITF  (6.300) 

CALL  INDT  (1  ,1  .NVALSZ.MTYPE»XVAL*IEOE) 

GO  TO  (190,160.170, 180,290.180) ♦  IEOE*l 
170    CONTINUE 

CALL  HELP  (20) 

GO  TO  160 
180    CONTINUE 

WPITE  (6.340) 

GO  TO  1*0 
190    CONTINUE 

IE  (-iTYPE.NE.32.ANn.MTYPE.NE.60)  GO  TO  180 
200    CONTINUE 

WPITE  (ft. 310) 
C      ♦«•   OPTION  EOP  ACCEPTING  A  RANDOM  SEED. 

CALL  INPT  (1  .1  ,NVA|_SZ,IS,XVAL.IEOE> 

GO  TO  (?40, 210, 220. 230, 290. 230) ,  IEOF*l 
210    CONTINUE 

WPITE  (6.330) 

GO  TO  200 
220    CONTINUE 

CALL  HELP  (21) 

GO  TO  200 
230    CONTINUE 

WPITE  (6.340) 

GO  TO  200 
240    CONTINUE 

IE  (IS.EO.O)  GO  TO  250 

IE  (IS.LT.1.0P.IS.GT.9)  GO  TO  230 

IU=IRAN(IS) 

GO  TO  260 
250    CONTINUE 

IU=67111133 
260    CONTINUE 

REWIND  NTEMP 
C   «»»WPITE  40  CHARACTER  NAME  EIELD  INTO 
C   ♦••FILE  SAMPLE  (NTEMP) 

WPITE  (NTEMP)  INAME.NMEAS 
C 

DO  280  I=1,NSAMP 
c      »«*   SIMULATE  A  PSEUDO  SAMPLE  CASE. 

CALL  SIMSU9  (MTYPE) 

DO  270  K=1,NMEAS 

V(K)=V(K)*ANTM(2.K) ♦ANTM(l.K) 
270    CONTINUE 

WPITE  (NTEMP)  (V(IO) ,IO=l.NMEAS) 
2b0    CONTINUE 
C   »»«SET  ISAMP  TO  1  TO  INDICATE  OPERATOR  SAMPLE  HAS  BEtN  GtNEHATED 

ISAMP=1 

END  Ell  E.  NTEMP 
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REWIND  NTEMP 

REWIND  NTMEAS 

WPITE  (6*360) 
290         CONTINUE 

RETURN 
C 

300    FORMAT  ("  ENTER  MACHINE  WORD  SIZE  (32  OR  60  bIT)  — '•) 
310    FORMAT  ('•  ENTER  RANDOM  NO  SEED  (0  IF  NONE  Or?  1-9) — '•) 
3?0    FORMAT  ('•  INPUT  SAMPLE  SIZE  — ") 
330    FORMAT  ("  NO  INFORMATION  ENCOUNTERED") 
3*0    FORMAT  ("  INVALID  RESPONSE") 

350    FORMAT  (•'  SAMPLE  SIZE  MUST  BE  GREATER  THAN  0") 
3*0    FORMAT  ("  SAMPLE  DATA  SAVED  ON  FILE  SAMPLE  OR  TAPt  3.") 

END 

SUBROUTINE  PRINT 
C 

C   »»»OUTPUT  OF  SIMULATED  OPERATORS 
C 

C  »«»NOTE:  ALL  OPERATOR  MEASUREMENTS 

C   ♦•^ARE  STORED  ON  FILE  SAMPLE.   PRESENTLY  EACH  RECORD 
C  *••  CONTAINS  NMEA'"  ANTHROPOMETRIC  MEASUREMENTS 
C 

C0MM0N/MCSM/C0RR(78)  . ANTM(2,12)  ,  \/  (25)  .NMEAS.KN, 
1     SOROOT (78) ,X (25) , IU.NTMEAS, I  COUNT* 
?     INAmE(4) ,ITMEAS»NIMEAS»NTEMP 

COMMON/MISC/IOIN.IOOUT.ICARD 

INTEGER  FMJ 

01  MENS  I  ON  FMT(6) .NUM( 12) .TEMP (12) .INDEX (2) 

DATA  FmT/^»H(     .4H      .4H(6X).4H      .4H(F6.»4H3) )  / 

DATA  NUM/lHl,lH2,lH3,lH4,lH5,lH6,lH7,lH8,lHV,2H10»2Hll»2Hl2/ 

ISTOPE=0 

IOFLAG=0 
C 

C      ♦*♦   MENU  CHOICE  OF  OUTPUTS 
C 
100    CONTINUE 

WDITE  (6,630) 

CALL  INPT  (l,l,NVALSZ»IOPT,XVAL«IEOF> 

GOTO  (120.110.110,110.390,110).  IEOF*l 
110    CONTINUE 

WRITE  (6,650) 

GO  TO  100 
120    CONTINUE 

IF  (IOPT.LT.l.OR.IOPT.GT.^)  GO  TO  110 

IF  (IOPT.EO.l)  GO  TO  130 

IF  (IOPT.EQ.2)  GO  TO  150 

IF  (IOPT.E0.4)  GO  TO  390 

GO  TO  190 
C 

c      »««   PRINTOUT  OF  MEANS/STANDARD  DEVIATIONS 
C 
130    CONTINUE 

WPITF  (6,580) 

DO    1^0    IR0W=1,NMEAS 

WRITE     (6.640)     IPOW, (ANTM(JO,IROW) ,J0=1,2) 
WO         CONTINUE 

WRITE     (6,620) 

GO    TO    100 
C 

c      •«•   PRINTOUT  OF  ANTHROPOMETRIC  CORRELATIONS 
C 
150    CONTINUF 

WPITE  (6,590) 

WPITE  (6,660)  (NUM(IO) .IO=l,NMEAS) 

WPITE  (6,600) 

DO  180  TROW=l .NMFAS 
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DO  160  TCOL=l»NMEAS 

N=NDX(IPOw,ICOL) 

TFMP(ICOL>=CORR(N) 
160    CONTINUE 

IF  (IROW.GT.l)  GO  TO  170 

WP1TE     (6,610)     (TEMPdCOLJ  , ICOL=I*Ow,NMEAS> 

GO    TO    1*0 
170         FMT(2)=NUM(IROW-l) 

FMT(4)=NUM<NMEAS*1-IR0W) 

WRITE  (6,FmT)  (TEMP(ICOL)  »ICOL  =  I;?OW.NMEAS) 
ISO    CONTINUE 

WPITE  (6*620) 

GO  TO  100 
C 

C  »♦•  PRINT  OF  ANTHROPOMETRIC  MEASUREMENTS 
C 
190    CONTINUE 

WPITF  (6,430) 

IF  (IOFLAG.LT.l)  WRITE  (6,440) 

IOFLAG=l 
C 

C   •••IMPUT  CARPIAC"  RETURN  —  ENTIRE  SAMPLE  LIST 
C   •♦•input  SINGLE  NUMBER    —  SINGLE  OPERATOR 
C   •*»inPUT  TWO  NUMBERS      —  SEGMENT  OF  OPERATOR  LIST 
C 

CALL  INOT  (1»2»NVALSZ.INDEX»XVAL*IE0F) 

GO  TO  (230,270,200,210,100,220),  IEOF*l 
200    CONTINUE 

CALL  HFLP  (19) 

IOFLAG=0 

GO  TO  190 
210    CONTINUE 

WRIT?  (6,420) 

GO  TO  190 
220    CONTINUE 

I0LIST=NVALSZ-2 

WRITE  (6,450)  IOLIST 
230    CONTINUE 
C   •••INDEX  VALIDATION 

IF  (INDFX(I).GT.O)  GO  TO  240 

WRITE  (6,480) 

GO  TO  190 
240    CONTINUE 
C 

C   »»«TO  DETERMINE  WHETHER  FILES  WILL  BE  REWOUND 
C 

IF  (ISTORE.GT.O.AND.INDEX(l) .GT.ISTO&E)  GO  TO  250 

N»SKIP=INDEX(1)-1 

GO  TO  ?60 
250    CONTINUE 

NPSKIP=INDEX (l)-ISTORE-l 
260    CONTINUE 

ISTAPT=INDEX(1) 
C   •  ••if  TWO  INPUTS  ENCOUNTERED  -  OUTPUT  LIMITS 
C   «»»TF  ONE  INPUTS  ENCOUNTERED  -  OUTPUT  SINGLE  OPERATOR 

IF  (NVALSZ.EO.l)  GO  TO  280 

IF  (NVALSZ.GE.2)  GO  TO  290 
270    CONTINUE 

ISTORE=0 

ISTART=1 

IEND=999 

NRSKIP=0 

GO  TO  320 
280    CONTINUE 

IEND=ISTART 

CO  TO  -no 
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290    CONTINUE 

IF  < lNneX<l).LE. INDEX (2>)  GO  TO  300 

WPITE  (6*490) 

GO  TO  1Q0 
300    CONTINUE 

IEND=INDEX(2) 
310    CONTINUE 

IF  (INOEX(I) .GT.ISTORE.4NO.ISTORE.NE.0)  GO  TO  330 
320    CONTINUE 
C   «»»QEWIN0  FILE  • 

REWINO  NTEMP 

READ  (NTEMP)  INAME*NMEAS 

IF  (EOF(NTEMP) .EQ.l)  GO  TO  410 
330    CONTINUE 

If  (NPSKIP.LT.1)  GO  TO  350 
C   *»«WINO  FILE  FORWARD 

00  340  T=l ,NRSKIP 

READ  (NTEMP)  (V(K) *K=1,NMEAS> 

IF  (EOF(NTEMP) .EQ.l)  GO  TO  400 
340    CONTINUE 
350    CONTINUE 
C   *»*OUTPUT  HEADER 

WRITE  (6*530)  ( INAME  ( 10)  ,  10=1  ,<►) 

IF  (ISTART.EO.IENO)  WPITE  (6,540)  IENO 

IF  (ISTaRT.NE. IENO. AND. IEND.NE. 999)  WRITE  (6,550)  ISTART»IENO 

WPITF  (6*500) 

WPITE  (6,660)  (NUM(IO) ,I0=1,NMEAS) 

WRITE  (6,600) 
C 

C   »»*0UTPUT  DATA  FROM  FILE 
C 

ISKIP=0 

NFOUNO=0 

DO  360  I=ISTART,IEND 
C   »««SKIP  LINE  FEATURE 

IF  (I5KIP.E0.1)  WRITE  (6,460) 

ISKIP=0 

IF  (ISTART.NE.IEND.ANO.I.NE..IEND.AND.  <  1/5)  *>5.EQ.  I  >  ISKIP=1 

READ  (NTEMO)  (V(K) ,K=1 .NMEAS) 

IF  (EOF(NTEMP) .EQ.l)  GO  TO  370 

NFOUND=NFOUND*l 

ISTORE=I 

W»ITE  (6*510)  (V(K) ,K=1 ,NMEAS) 
360    CONTINUE 

GO  TO  ISO 
370    CONTINUE 

IF  (IENn.NE.999,AND.ISTART.NE.IEND)  WRITE  (6,470)  IENU 

I5TORE=I 

IF  (NFOUNO.LT.l)  WRITE  (6,560) 
380    CONTINUE 

WPITE  (6,520) 

GO  TO  190 
390    CONTINUE 

REWIND  NTEMP 

RFTURN 
400    CONTINUE 

I0LIST=I-2»IST0RE 

WPITE  (6*560) 

GO  TO  190 
410    CONTINUE 

WPITE  (6,570) 

GO  TO  100 
C 

<»20    FORMAT  ("  INVALID  RESPONSE") 
430    FOPMAT  ('•  OUTPUT  MODE: — ") 
t*uO  FORMAT  (••  FOR  AIL  OPERATOR  DATA  (ENTER  CARRIAGE  RtTU^N)",/, 
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1"  FOR  INDIVIDUAL  OPERATOR  DATA  (ENTEK  OPERATOR  NU^BE*) "./ 

2"  FOR  A  SEQUENCE  OF  OPERATORS  (E^TER  RANOE  OF  OREKATOR  NUMBERS)"/ 

♦  "  ENO  (ENTER  $)") 

450    FORMAT  ("  THE  LAST". 14,"  INPUTS  HAVE  BEEN  IGNORED") 

<*60    FORMAT  (1H  ) 

470    FORMAT  ("  ERROR  IN  UPPER  INDEX  —  ".14) 

480    FORMAT  ("  INOEX  NUMBER  CANNOT  BE  LESS  THAN  OR  EQUAL  TO  0") 

490    FORMAT  (•'  FIRST  INDEX  MUST  BE  LESS  THAN  SECOND  INUEX") 

500    FORMAT  (27x, "OPERATOR  SAMPLE'S/* 71 ( 1H-) , ) 

510    FORMAT  (11F6.2,F5.2> 

520    FORMAT  (71 (1H-)) 

530    FORMAT  ( 1 X  ,4A10 ♦ 10X . A10) 

540    FORMAT  (10X, "OPERATOR  NO"»I4) 

550    FORMAT  ( 1  Ox ."OPERATOR  NOS",I4,"  — ".14) 

560    FORMAT  ("  INDEX  NUMBER  EXCEEDS  RECORDS  IN  OPERATOK  FILE") 

570    FORMAT  ("  FILE  SAMPLE  IS  NOT  ATTACHED  TO  THE  PROGkAM"/ 

*  ■•  OPFRATOR  OUTPUT  PEQUEST  APORTED ".//) 

5*0    FORMAT  (/.21X. "MEANS  AND  STANDARD  DEV I  AT  IONS" »/♦ IX ,70 < 1H-) ,/,2X ♦ "M 

1EASUREmENT".6X. "MEANS". 6X."ST0  DEV",/.1X,71 (1H-) ) 
590    FORMAT  (/.7X. "CORRELATION  MATRIX  FOR  ANTHOPOMETRIC  MEASUREMENTS"./ 

1.1X.70UH-)  ) 
600    FORMAT  (lX,70(lH->) 
610    FORMAT  (12F6.3) 
6?0    FORMAT  (1X,70(1H-)) 

*30    FORMAT  ("  PRINT  MODE  ( 1-MEAN/STD52-CORP  MATRXI3-OPER  mEAS'.4-END)  -"> 
640    FORMAT  (I9.5X.2F12.2) ) 
650    FORMAT  ("  INVALID  RESPONSE") 
660    FORMAT  (2X,  11 (A2.4X) »A2) 

END 

SUBROUTINE  RAND2  (RNORM.MTYPE.QUAN, IJ) 
C 

C   •♦•GENERATES  A  NORMAL  DISTRIBUTION  *ITH  MEAN  OF  ZERO  AND 
C   •♦•STANDARD  DEVIATION  EQUAL  TO  ONE.... 

C   ••♦METHODOLOGY  MARSAGLIA  -  BRAY 

C 

C0MM0N/MCSM/C0PR(7R)  »ANJTM(2.12)  .«/  (25)  .NMLAS.KN, 

1  SQR00T(78)  .X (25) . IU.NTMEAS » I  COUNT  * 

2  INAME (4) .ITMEAS.NIMEAS.NTEMP 
DIMENSION  QUAN(2) 
COMMON/MISC/IOIN.IOOUT.ICARD 
IJ=U*1 

IF  (IJ.FQ.2)  GO  TO  110 
100    CONTINIJF 

X1=2.0*UNFORM(MTYPE)-1.0 

Y1=2.0»UNFORM(MTYPE)-1.0 

S=X1*X1*Y1»Y1 

IF  (S.GE.1.0)  GO  TO  100 

S=SQPT(-2.0*ALOG(S)/S) 

QUAM(l)sXl«S 

0UAN(2)=Y1»S 
110    CONTINUE 

RNORM=OUAN(IJ) 

IF  (IJ.FQ.2)  IJ=0 

RETURN 

END 

SUBROUTINE  REASON  ( ICALL  .  ICOL  ♦  IROW.TSTORE1  .TSTORE<i»  IFAtAL) 
C 

C      ••*   CHECKS  FOR  REASONABLENESS  OF  INPUT  VALUES 
C 

C0MM0N/MCSM/C0RR(78> .ANTM(2,12) .V(25) .NMEAS.KN. 
1     SQPOOT (78) .X(25)  . IU.NTmEAS, I  COUNT ♦ 
?     INAME(4) .ITMEAS.NIMEAS.NTEMP 

COMMON/MISC/IOIN.IOOUT.ICARD 
C 

C   ICALL  =  1   MEANS  AND  STANDARD  DEVIATIONS 
C.       TC.AU   =  ?       COP.RFI  ATTON  MATPT.X.  VERIFICATION 
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C   ICALL  =  3   INDEX  VERIFICATION 
C 

IFATAL=0 

IAFLAG=0 

IF  (ICALL. LT.O)  IAFLAG=1 

ICALL=IABS(ICALL) 

IF  (ICALL. GT.l)  GO  TO  110 
C 

C      *•*   mfaN  VALUES  MUST  8E  WITHIN  A  THREE  STANDARD  DEVIATION 
C      »•«   NEIGH90RH00D  OF  THE  ACCEPTED  MEAN  VALUES 
C 

IWRITE=1 

IF  (IROW.LT.l.OR.IROW.GT.NMEAS)  GO  TO  120 

IF  (ICOL.E0.2)  GO  TO  100 

RETURN 
C 

C      •**   STANDARD  DEVIATIONS  MAY  NOT  BE  NEGATIVE 
C      ••»   ANO  THE  MEANS  MUST  BE  GREATEH  THAN 
C      *•*   THE  STANDARD  DEVIATION. 
C 
100    CONTINUE 

IWRITE=2 

IF  (TSTORE2.LT.0.)  GO  TO  120 

IWRITE=3 

IF  (TST0RE1.LE.TST0RE2)  GO  TO  120 

RETURN 
C 

C      »««   CHECKS  VALIDITY  OF  INDEX  NUMBERS 
C 
110    CONTINUE 

IWRITE=1 

IF  (IROW.LE.O.OR.IROW.GT.NMEAS)  GO  TO  130 

IWRITE=? 

IF  (ICOL.LE.O.OP.ICOL.GT.NMEAS)  GO  TO  130 
C   •••IS  THE  NUMBER  WITHIN  THE  UPPER  TRIANGULAR  REGION? 

IWRITE=3 

IF  (IPOW.GT.ICOL)  GO  TO  130 

IF  (ICALL. EO. 3)  RETURN 
C 

C      •••   ALL  DIAGONAL  ELEMENTS  MUST  EUUAL  1.000 
C      •♦•   MATRIX  ELEMENTS  MAY  NOT  BE  LESS  THAN  (-1)  OR  GREATER 
C      •••   TmaN  ONE. 
C 

IWRITE=4 

IF  (ICOL. EQ.IROW. AND. TSTORE1.NE. 1.000)  GO  TO  130 

IWRITF=5 

IF  (TSTOREl.LT. -1.000)  GO  TO  130 

IWRITE=6 

IF  (TSTORE1. GT.l. 000)  GO  TO  130 

RFTURN 
C 

C   ««*BRANCH  TO  APPROPRIATE  ERROR  MESSAGE 
C 
120    CONTINUE 

IF  (IWRITE.NE.1.AND.I0IN.NE.5)  W-JITE  (6.280)  IROW 

GO  TO  (190,140,150),  IWRITE 
130    CONTINUE 

GOTO  (190,200,210,160,170,180),  IWRITE 
140    CONTINUE 

WPITE  (6.230) 

GO  TO  220 
150    CONTINUE 

WRITE  (6*240) 

GO  TO  220 
160    CONTINUE 

WPTTF  f*.?SO)  .T.ROW.T<~m 
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170 

i«o 

190 
200 
210 
220 


C 

230 
2*0 
250 

260 

270 

2S0 
290 
300 
310 


IROW.ICOL 


IROW, ICOL 


IROW,NMEAS 


ICOL,NMEAS 


IPOW,ICOL 


C   ** 

c   •*• 


c 

100 


GO  TO  220 
CONTINUE 
WPITF  (6*260) 
GO  TO  220 
CONTINUE 
WRITE  (6.270) 
GO  TO  220 
CONTINUE 
WRITE  (6.290) 
GO  TO  220 
CONTINUE 
WPITE  (6.300) 
GO  TO  220 
CONTINUE 
WRITE  (6.310) 
GO  TO  220 
CONTINUE 
IFATAL=1 
RETURN 

FORMAT  ('•  STANDARD  DEVIATIONS  LESS  ThAN  ZERO") 

FORMAT  (••  STa.JDARD  DEVIATION  GREATER  THAN  MEAN  VALUE.") 

FORMAT  (•*  DIAGONAL  ELEMENT  IN  RO*  ".12."  COLUMN  ".12. 

♦  "  IS  NOT  EQUAL  TO  1.0") 

FORMAT  (••  ENTRY  IN  COLUMN  ".12."  ROW  ",I2. 

♦  "  IS  LESS  THAN  -1.0") 

FORMAT  (••  ENTRY  IN  COLUMN  ".12."  HOW  ",I2» 

♦  "  IS  GREATER  THAN  ♦1.0") 
FORMAT  (110) 

FORMAT  ("  ROW  ".12."  IS  LESS  THAN  ZERO  OR  GREATER  THAN  ",I2) 
FORMAT  ("  COLUMN  ".12."  IS  LESS  THAN  ZEPO  OR  GREATER  THAN  ".12) 
FORMAT  ("  ROW  ",I2,"  COLUMN  ".12."  IS  OUTSIDE  OF  UPPER", 

♦  "  TRIANGULAR  REGION  OF  CORRELATION  MATRIX.") 
END 

SUBROUTINE  SAVE 

••*   THIS  ROUTINE  SAVES  INPUT  DATA  ON  PERMANENT  FILE  TAPE2 

COMMON/MCSM/C0RP(78) ,ANTM(2.12) .V(25) .NMEAS.KNt 

1  S0R00T(78) »X(25) » IU.NTMEAS. I COUNT. 

2  INAME(^) .ITMEAS.NIMEAS.NTEMP 
COMMON/MISC/IOIN.IOOUT.ICARD 
REWIND  NTMEAS 

•SAVE  NO*  OF  MEAS.,  ELEMENTS  IN  COkR  MATRIX.  SEX  FLAG 
WRITE  (NTMEAS)  NMEAS.KN. INAME 

♦  SAVE  MEAN/STD.  VECTOR. 

WRITE  (NTMEAS)  ( ( ANTM ( 10, JO) ♦ 10=1 .2) . JO=l .NMEAS) 

SAVE  CORRELATION  MATRIX 

WRITE  (NTMEAS)  (CORR ( 10) . 10=1 ,KN) 

REWIND  NTMEAS 

WRITE  (6.100) 

RETURN 

FORMAT  (»  MEASUREMENT  DATA  SAVED  ON  FILE  UMEAS  OR  TAPE2") 

END 

SUBROUTINE  SIMSUB  (MTYPE) 


•  «« 
««« 
*»« 
•«» 
•«» 


FOP  A  GIVEN  SOUAR  ROOT  MATRIX  R  OF  ORDER  M  IN  UPPER 
TRIANGULAR  FORM  OF  A  CORRELATION  OR  A  COVARIANCE 
MATRIX  OF  M  COLUMNS,  THIS  ROUTINE  COMPUTES  A  VECTOR 
V  OF  PSEUDO  SCORES  OF  LENGTH  M  SUCH  THAT  IT  WILL 
BE  A  QUASI-RANDOM  SAMPLE  FROM  A  NORMAL  POPULATION 
WITH  A  MEAN  VECTOR  OF  ZEROS  AND  THE  SAME  CORRELATION 
OR  COVARIANCE  MATRIX,  I.E.  FROM  AN  N(O.R)  MATRIX. 


COMM0N/MCSM/CORR(7ft)  .ANTM  (2,1  P)  ,V(?S)  .NJMFAS.KN. 
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c 

c 

»»« 

c 

««« 

c 

*«« 

c 

»«« 

c 

1  S«R00T(78>  .X(25)  . IU.NTMEAS . I  COUNT  * 

2  INflME  U>  .ITMEAS.NIMEAS.NTEMP 
COMMON/MISC/IOIN.IOOUT.ICARD 
DIMENSION  Y<45) 

IJ=0 

DO  100  I=1.NMEAS 

V(I)=0.0 

CALL  RAND2  (RNORM.MTYPE.QUANt IJ) 
100    Y(I)=RNORM 

DO  110  I=ltNMEAS 

DO  110  J=1.I 

K=NDX( J, I) 
110    V(I)=V(I)*Y(J)»SQROOT(K) 

RETURN 

END 

SUBROUTINE  SORMX  (NERROR) 

FOP  A  GIVEN  SYMMETRIC  MATRIX  OF  ORDER  N  STORLD  IN  ITS 
UPPER  TRIANGULAR  FORM.  COLUMNWISE  IN  THE  VECTOR  H  OF 
LENGTH  N»(N*l)/2»  THIS  ROUTINE  TRANSFORMS  IT  INTO  ITS 
"SQUARE  ROOT"  IN  THE  SAME  FORM  IN  THE  /LCTOR  S0*00T<78) 

C0MM0N/MCSM/C0RR(7fl) .ANTM(2»12) tV(25) .N^EAStKN* 

1  SQR00T(78> .  x (25)  . IU.NTMEAS. I  COUNT* 

2  INAME(4) ,ITMEAS.NIMEAS.NTEMP 
COMMON/MISC/IOIN.IOOUT.ICARD 
DIMENSION  ZX(7) 

INTEGER  ZX 
C   •••NERROR  FLAG  IS  USED  TO  SIGNIFY  A  CONDITION 
C   ••♦WHICH  WILL  PRODUCE  A  NEGATIVE  SQRUOT(N)  VALUE 
C   •••JUST  BEFORE  IMPLEMENTATION  OF  THE  SQUARE  ROOT 
C   •••FUNCTION. 

NFRRORsf) 

DO  120  I=ltNMEAS 

ZX(4)=I 

N=NDX(I,I) 
C   ••♦TFST  FOP  NEGATIVE  RADICAL  ARGUMENT 

IF  (SQROOT(N) .LT.O)  GO  TO  130 

SQ=SQRT(SOROOT(N) ) 

DO  100  J=I.NMEAS 

N=NDX(T»J) 

SQROOT (N) =SOROOT (N) /SO 
100    CONTINUE 
C   •••NORMAL  EXIT  FROM  SUBROUTINE 

IF  (I.FO.NMEAS)  NERROR=0 

IF  (I.EO.NMEAS)  RETURN 

11=1*1 

ZX(5)=T1 

DO  110  I2=I1.NMEAS 

00  110  J2=I2.NMEAS 

N=NDX (I2»J2) 

L=NDX(I.I2) 

K=NDX(T,J2) 

SOHOOT (M) =SQROOT (N) -SOROOT (L) «SQROOT (K) 

ZX(6)=I2 

ZX(7)=J2 
110    CONTINUE 

ZX(1)=N 

ZX(2)=L 

ZX(3)=K 
120    CONTINUE 
130    CONTINUF 
C   •••ABNORMAL  EXIT  FROM  SUBROUTINE 

NFRRO«=l 

WRITE  (<S»190) 
)40    CONTINUE 
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WPITE  (6.190) 

CALL  INPT  (4,1 .NVALSZ.IANS.XVAL.IEOF) 

GO  TO  < 170,160.150,170.160,170) ,  1E0F*1 
150    CONTINUE 

WPITE  (6,200) 

WPITE  (6»210)  (ZX(I) ,1=1,7) ,SQROOT(N) 
160    CONTINUE 

RFTURN 
170    CONTINUE 

WPITE  (6,220) 

GO  TO  140 
C 

180    FORMAT  (•'  ANTHROPOMETRIC  CORRELATION  VALUES  ARE  INVALIO") 
190    FORMAT  ('•  —  ENTER  CARRIAGE  RETURN  TO  CONTINUE  PROCESSING",/ 

♦  »  —  FNTEP   ?   TO  DISPLAY  SQRMX  VARIABLES") 

200    FORMAT  (<•  AN  INDEFINITE  QUANTITY  HAS  BEEN  DETECTED  BY  SQRmx"/ 

♦  '•  VARIABLE  LIST  WITHIN  SQRMX:"/"   N     L", 

♦  '•     K     I     II    12    J2    SQ-?OOT"»/,10X,44(1H-)  ) 
210    FORMAT  (8X,4I5,1X,3I5,F9.3,/,10X,44(1H-) ,//) 

220    FORMAT  (/'•  INVALID  RESPONSE  ••*••/) 

END 

FUNCTION  NDX (I  ,J) 

IF  (I.GT.J)  GO  TO  100 

NDX=I*(J»J-J)/2 

RFTURN 
100    NDX=J» (I»I-I)/2 

RFTURN 

EWD 

FUNCTION  UNFORM(MTYPE) 
C 
C   »«»UNlFOi*M  RANDOM  NUMBER  GENERATOR  ON  INTERVAL  (0»1) 

C 

COMMON  /MCSM/  CORP(78) »ANTM(2,12) »V(2b) ,NMEAS»KNt 
1S0R00T(78)  ,X (25)  , IU.NTMEAS. I  COUNT, 
?INAMF(4) ,ITMEAS,NIMEAS,NTEMP 

COMMON/MlSC/IOIN,IOOUT,ICAPD 
C       60  BIT  UNIFORM  RANDOM  NUMBER  GENERATOR 

DATA  I x/ 1 677721 3/, C/28 14749767 10655./ 

IF  (MTYPE.EQ.32)  GO  TO  100 

IU=IU*IX 

UNFOPM=iu/C 

RFTUPN 
C       32  BIT  UNIFORM  RANDOM  NUMBER  GENERATOR 
100    IU=IU»65539 

IF  (IU.GE.O)  GO  TO  110 

IU=IU*2147483647«1 
110    UNFOCm-iu 

UNFORM=UNFORM«. 46566 13E-09 

RETURN 

END 

SUBROUTINE  HELP  <I> 

RETURN 

END* 
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APPENDIX  D 

Program  Listing  for  CAPE  Model 
Modified  for  Use  with  FORTRAN  IV' 


1  C0MM0N/MCSM/C0RR178)  ,ANTM(2.12)  ,V(25)  .NMEflS.KN. 

2  1    SQR00T(78)  .  X  (25)  , 1U. NTMEAS, ICOUNT, 

3  2    INAMEUO)  .ITMEAS.NIMERS.MAXSAM.NTEMP 

4  COMMON/MISC/IQIN.IGQUT.ICARD 

5  DflTfl  NUCREW/8/ 

6  DflTfl  MflXFUN/4/ 

7  DflTfl  MflXOPT.MAXANT.MAXflCT/3,7,3/ 

8  IQ0UT*O 

9  ISAMP=0 

10  MMEflS  =  12 

11  100   CONTINUE 

12  I0IN=5 

13  WRITE  (6,290) 

14  CALL  INPT  (1.1.NVALSZ.NMEAS.XVAL.IE0F) 

15  IE0FP1-IE0F+1 

16  GOTO  (130.100,110.120.280.120).  IEOFPl 

17  110   CONTINUE 

18  CALL  HELP  (1) 

19  GO  TO  100 

20  120   CONTINUE 

21  NRITE  (6,340) 

22  GO  TO  100 

23  130  CONTINUE 

24  IF  (NHEAS.LT.1.0R.NMEAS.GT.MMEAS)  GO  TO  120 

25  KN=NMEflS*(NMEAS+l) /2 

26  11NPT=0 

27  140   CONTINUE 

28  I0IN=5 

29  150   CONTINUE 

30  WRITE  (6.300) 

31  CALL  INPT  (l.l.NVALSZ.IOPT.XVAL.IEOF) 

32  IE0FP1=IE0F+1 

33  GOTO  (180,140.160.170.280.170) .IEOFPl 

34  160  CONTINUE 

35  CALL  HELP  (2) 

36  GO  TO  150 

37  170   CONTINUE 

38  WRITE  (6.340) 

39  GO  TO  150 

40  180   CONTINUE 

41  IF  (lOPT.LT.l.OR.IOPT.GT.MAXANT)GO  TO  170 

42  GOTO  (190.200,220.240.260.280.225) .IOPT 

43  190   CONTINUE 

44  CALL  ACCPT  (ISAMP) 

45  IINPT=1 

46  GO  TO  140 

47  200  CONTINUE 

48  IF(IINPT.EQ.l)  GO  TO  210 

49  WRITE  (6.350) 

50  GO. TO  140 
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51 

210 

52 

53 

5i 

220 

55 

56 

222 

57 

58 

225 

59 

60 

61 

62 

230 

63 

61 

65 

240 

66 

67 

68 

69 

250 

70 

71 

72 

260 

73 

71 

75 

76 

270 

77 

78 

79 

280 

80 

81 

82 

290 

83 

300 

81 

85 

310 

86 

320 

87 

88 

330 

69 

340 

90 

350 

91 

92 

360 

93 

91 

CONTINUE 

CALL  EDIT 

GO  TO  mo 

CONTINUE 

IF  (1SAMP.EQ. 1)  GO  TO  230 

WRITE  (6.320) 

GO  TO  110 

CONTINUE 

IF  (ISAMP.NE.l)  GO  TO  222 

CALL  PUNRT 

GO  TO  110 

CONTINUE 

CALL  SAVE 

GO  TO  mo 

CONTINUE 

IF  (IINPT.EQ. 1)  GO  TO  250 

WRITE  (6.360) 

GO  TO  110 

CONTINUE 

CALL  PRINT 

GO  TO  110 

CONTINUE 

IF  (ISAMP.NE.l)  GO  TO  270 

CALL  OPSAMP(ISAMP) 

GO  TO  110 

CONTINUE 

WRITE  (6.330) 

GO  TO  110 

CONTINUE 

WRITE  (6.310) 

STOP 

FORMAT  C  ENTER  NUMBER  OF  MEASUREMENTS  (1-12)?--') 

FORMAT  C  0PTI0NS(1-1NPUT;2-ED1T;3-SAVE:1-PRINT55-GENERATE;  ' 
1  .  '6-END;7-PUNCH)  ') 

FORMAT  (5X, 'END  OF  OPERATOR  SAMPLE*) 

FORMAT  (IX.  'xxxMCSMUO  UNABLE  TO  SAVE  MEASUREMENT  DATA  -  FATAL 
1  ERRORS') 

FORMAT  (10X, '***MCSM165  UNABLE  TO  GENERATE  SAMPLE  -  FATAL  ERRORS') 

FORMAT  (/.  '**hMCSMQ21  INVALID  RESPONSE') 

FORMAT  (10X.  '**«MCSM130  ANTHROPOMETRIC  DATA  MUST  BE  INPUT  BEFORE'. 
1/.13X, 'THE  EDIT  FUNCTION.') 

FORMAT  (10X.  '*««MCSM150  ANTHROPOMETRIC  DATA  MUST  BE  INPUT  BEFORE*. 
1/.13X.'THE  PRINT  FUNCTION.') 

END 
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1  SUBROUTINE  INPT  (I TTPE. 1 VALSZ. NVALSZ. I VAL. XVAL. IEOF) 

2  COMHON/MISC/IOIN.IOOUT.  ICARD 

3  DIMENSION  I  CHAR  (80) , I VAL (I VALSZ)  . XVAL  (1 VALSZ) 

4  C 

5  DIMENSION    NVflL  (20) 

6  LOGICAL    m    1    NCHftR  (80) 

7  EQUIVALENCE     (NVflL  (1)  ,  NCHAR  (1)  ) 

8  C 

9  C 

10  DATA  IBLK.lQUES.IDOL.lQUO.lPOST.IPD.ICOM.lMlN.lPLUS/lH  .  1H?. 

11  MH$.1H',1Hi,1H..1H.,1H-,1H+/ 

12  DATA  14BLK/4H     / 

13  IST0P=O 

14  IF  (IVflLSZ.GE.O)  GO  TO  100 

15  ISTDP=1 

16  IVflLSZ=IABS(lVALSZ) 

17  100      CONTINUE 

18  IQU0TE=O 

19  NVflLSZ-Q 

20  IEOF-0 

21  IF    (100UT.EQ.1)    GO   TO    120 

22  00   110    1=1,30 

23  ICHflR  (11  -J8LK 

24  110      CONTINUE 

25  READ     (IOIN, 330. END=111) ICHflR,    NCHflR 

26  DO   109    J=1.60 

27  IF  HCHflR(J)  .EQ.  1BLK)  GO  TO  109 

28  GO  TO  130 

29  109  CONTINUE 

30  GO  TO  113 

31  111  CONTINUE 

32  REWIND  IOIN 

33  113  IEOF  -  1 

34  RETURN 

35  120   CONTINUE 

36  Il-ICARD 

37  GO  TO  140 

38  130   CONTINUE 

39  11=1 

40  140   CONTINUE 
VI  I2=B0 

42  150   CONTINUE 

43  C  HRITE(6.501)  IOEC. ISTOP, I VALSZ. NVALSZ.  11.12,  IEOF 

44  IF  (ISTOP.  EQ.  LAND.  IVALSZ.EQ.  NVALSZ)  RETURN 

45  IF  (11. IE. 12)  GO  TO  160 

46  IF  (NVALSZ. LE.IVALSZ)  RETURN 

47  IEOF-5 

48  RETURN 

49  160   CONTINUE 

50  DO  170  1=11,12 


85 


51 

IF  (ICHAR(l)  .EQ.1BLK  )  GO  TO  170 

52 

IFIRST=I 

53 

GO  TO  180 

54 

170 

CONTINUE 

55 

IF  INVALSZ.EQ.O)  1E0F-1 

56 

IF  (NVALSZ.GT.1VALSZ)  IE0F=5 

57 

RETURN 

58 

180 

CONTINUE 

59 

IF  (ICHflR(lFIRST) .EQ.IQUES)  IE0F=2 

60 

IF  (ICHAR  (1FIRST) .EQ.1D0L)  IE0F=4 

61 

IF  (IEOF.NE.O)  RETURN 

62 

IF  (ITYPE.EQ.4)  GO  TO  300 

63 

IF  (ITTPE.NE.3)  GO  TO  190 

64 

IF  (1CHARUFIRST)  . NE. 1QU0. AND. ICHAR  (IF IRST)  .NE.IPOST) 

65 

IF1RST-1FIRST+1 

66 

1QU0TE=1 

67 

190 

CONTINUE 

68 

DO  220  I-IF1RST.I2 

69 

IF  (IQUOTE.EQ.O)  GO  TO  200 

70 

IF  (ICHAR(I)  .EQ.IQUO.OR.ICHAR(I)  .EQ.IPGSTJGO  TO  210 

71 

GO  TO  220 

72 

200 

CONTINUE 

73 

IF  (ICHAR(l) .NE.IBLK.AND.ICHAR(I) .NE.ICQM)  GO  TO  220 

74 

210 

CONTINUE 

75 

ILAST-1-1 

76 

GO  TO  230 

77 

220 

CONTINUE 

78 

1LAST=12 

79 

230 

CONTINUE 

80 

NVALSZ=NVALSZ+1 

81 

I1*ILAST*2 

82 

IF  (NVALSZ.GT.IVALSZ)  GO  TO  150 

83 

IF  (ITTPE  .EQ.3)  GO  TO  310 

84 

ISIGN-1 

85 

IF  (ICHAR  (IF IRST)  . NE. 1  Ml N. AND. ICHAR  (IF IRST) .NE.IPLUS) 

86 

IF  (ICHAR(IFIRST) .EQ.IMIN)  IS1GN=-1 

87 

1F1RST=IFIRST+1 

88 

240 

CONTINUE 

89 

INUM-O 

90 

IDEC-0 

91 

IDECF=0 

92 

DO  260  I=IFIRST,1LAST 

93 

IF  (ICHAR  (I) .NE.IPD)  GO  TO  250 

94 

IDECF=IDECF-M 

95 

1DEC=1 

96 

GO  TO  260 

97 

250 

CONTINUE 

98 

CALL  1C0NV  (ICHAR  (1) .1ADD.IE0F) 

99 

IF  (IE0F.EQ.3)  RETURN 

100 

INUM=1NUM«10+1ADD 

GO   TO    190 


GO   TO   240 
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ioi 

501 

FORMAT!/  ■   TEST  \  101  10) 

102 

IF  (IDEC.NE.O)  1DEC=1DEC+1 

103 

260 

CONTINUE 

104 

IF  tIDECF.LT.2)  GO  TO  280 

105 

270 

CONTINUE 

106 

C 

WRITE  (6.501)  10EC.  10ECF 

107 

1E0F-3 

108 

RETURN 

109 

280 

CONTINUE 

110 

INUM=1NUM«ISIGN 

111 

XNUM=FLOAT  (INUM) 

112 

CC 

WRITE  (6.501)  IOEC.  IDECF.  INUM.  1AD0.  ISIGN.  ITTPE 

113 

IF  ((ITTPE. EQ. l.ANO.IDEC.LE.l) .OR.  (ITTPE. EQ. 2))  GOTO 

im 

1EQF-3 

115 

RETURN 

116 

290 

CONTINUE 

117 

IF  (IDEC.NE.O)  XNUM=XNUM/10.0k*  (IDEC-1) 

118 

IF  (ITTPE. EQ.l)  IVAL (NVALSZ) =IF1X  (XNUM) 

119 

IF  (ITTPE. EQ. 2)  XVAL  (NVALSZ) =XNUM 

120 

1CARD-ILAST+1 

121 

C 

WRITE  (6.501)  IOEC. IDECF. INUM. IflDO.  ISI GN. 1TYPE, NVflLSZ. 

122 

C 

1   IVflL  (NVflLSZ)  ,  ICflRD 

123 

GO  TO  150 

124 

300 

CONTINUE 

125 

ILAST=IVALSZ*4 

126 

909 

FORMAT  (  2X.  80A1) 

127 

910 

FORMAT  (2X.  20A4) 

128 

310 

CONTINUE 

129 

DO  320  J-1.IVALSZ 

130 

IVAL  (JJ-I4BLK 

131 

320 

CONTINUE 

132 

L-ILAST-IF1RST+1 

133 

NVALSZ"  (L-l)  /10-H 

134 

IF  (L.GT.4*IVALSZ)  GO  TO  270 

135 

C 

ENCODE  (L. 340, IVAL  (1)  )  (ICHAR  (1 ) . IF IRST. ILAST) 

136 

J  -  0 

137 

904 

FORMAT  C   ENCODE  \  101 10) 

138 

00  311  I  =  IFIRST. ILAST                   , 

139 

J  -  J  ♦  1 

140 

NCHAR(J)  =  NCHAR(I) 

141 

311 

CONTINUE 

142 

DO  313  I  «  l.IVALSZ 

1<43 

IVAL  (I)  ■'  NVAL  (I) 

144 

313 

CONTINUE 

145 

I CARD- 11 

146 

RETURN 

147 

330 

FORMAT  (80A1.  Tl.  80A1) 

148 

340 

FORMAT  (80A1) 

149 

END 

290 
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1  SUBROUTINE  ICONV  UCHflR. 1 1 . IEOF) 

2  DIMENSION  JNUM (10) 

3  DflTR  JNUM/1H0. 1H1, 1H2.1H3, 1HM.1H5.1H6.1H7.1H8.1H9/ 
M          11  FORMAT  (//'   RT.  ICONV   \  R2.2I10  /) 

5  DO  100  >1.  10 

6  IF     UCHflR.  NE.  JNUM  (J)  )     GO    TO    100 

7  I 1*J-1 

8  RETURN 

9  100   CONTINUE 

10  IEOF  =»  3 

11  WRITE (6, 111  JCHflR.  11.  IEOF 

12  RETURN 

13  END 
14 
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1  SUBROUTINE  ACCPT  (ISflMP) 

2  COMMON/MCSM/CORR (78)  ,RNTM  (2. 12) . V  (25) .NMERS.KN. 

3  1    SQROQT  (78) .X  (25)  , 1U, NTMERS, ICOUNT. 

U  2    1NRME  (10)  , ITMERS.N1MEAS.MAXSAM.NTEMP 

5  C0MM0N/MISC/10IN.100UT.ICARD 

6  DIMENSION  BUFF  (3)  .TEMP  (12) 

7  DRTR  BUFF/3*0./ 

8  1SRMP=0 

9  ICOUNT=0 

10  INL1NE=0 

11  I01N=5 

12  100   CONTINUE 

13  WRITE  (6.610) 

1IJ  CALL  1NPT  (1.1.N0INPT.IM00E.XINPT.IE0F) 

15  1E0FP1-IE0F+1 

16  GOTO  (130.110.120.110.1190.110)  .IE0FP1 

17  110   WRITE  (6.650) 

18  GO  TO  100 

19  120   CALL  HELP  (4) 

20  GO  TO  100 

21  130   IF  (1M00E.EQ.1)  GO  TO  200 

22  IF  (IM00E.NE.2)  GO  TO  110 

23  140   CONTINUE 

21  WRITE  (6,620) 

25  CRLL  INPT  (1 . 1 . NOINPT . IFFORM. XINPT . IEOF) 

26  1EDFP1=IEQF+1 

27  GOTO  (170.150.160.150.1490.150).  IEOFPl 

28  150  WRITE  (6.650) 

29  GO  TO  140 

30  160   CALL  HELP  (5) 

31  GO  TO  mo 

32  170   IF  (IFFORM. EQ.l)  GO  TO  190 

33  IF  (IFFORM. NE. 2)  GO  TO  150 
3M  REWIND  NTMERS 

35  RERD  (NTMERS. END=180)  NMF.AS. KN. INRME 

36  RERD    (NTMERS)      (  (RNTM  (10. JO) . 10= 1 . 2)  . J0=1 . NMER5) 

37  RERD    (NTMERS)      (CORR  (10)  . 10= 1 . KN) 

38  REWIND  NTMERS 

39  WRITE  (6.690) 
•10  ISAMP=1 

HI  RETURN 

U2  180  CONTINUE 

■43  WRITE  (6.670) 

44  GO  TO  490 

115  190   I0IN=1 

•16  200   CONTINUE 

•17  CALL  NRMEFL  (IERROR) 

48  1C0UNT=IC0UNT+IERR0R 

119  C     WRITE  (6.901)  1M0DE.  IERROR.  ICOUNT 

50  901  FORMAT!/  '  RCCPT  \  101  10) 
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51  INLINE- INLINE* 1 

52  IF  (IMODE.EQ.l)  WRITE  (6.500) 

53  IF  (1MG0E.NE.1)  WRITE  (6.510)  1NRME 

54  MfiXlNP-2 

55  C  WRITE(6,90I)  IMODE,  IERROR. 1C0UNT. INLINE. MRXINP. 1E0F 

56  DO  320  IR0W=1,NMERS 

57  210  CONTINUE 

58  IF  (IMODE.EQ.l)  WRITE  (6.530)  IROW 

59  220   NRMT=0 

60  230  CONTINUE 

61  1T0TIN=MRXINP-NRMT 

62  IC0L=NRMT*1 

63  CALL  INPT  (2. I TOTIN. NOINPT. 1 VRL, RNTM  (ICOL. IROW)  . IEOF) 

64  IF  (IE0F.EQ.5.0R. IEOF.EQ.O)  NRMT=NRMT*NOINPT 

65  C  HRITE(6,901)  IMOOE.  ITOTIN.  MRXINP. NRMT. IEOF, ICOL. ICOL 

66  IF  (IMODE.EQ.l)  GO  TO  240 

67  1NL1NE=INL1NE+1 

68  IF  (NRMT. EQ. MRXINP)  GO  TO  290 

69  IF  (1E0F.EQ.1)  GO  TO  480 

70  WRITE  (6,540)  IROW 

71  IF  (IE0F.EQ.5)  GO  TO  270 

72  WRITE  (6,700)  INLINE 

73  I  COUNT* 1  COUNT +1 

74  GO  TO  320 

75  240   CONTINUE 

76  IEOFPl=IEOF+l 

77  GOTO  (280.210.250.260.490.270).  IE0FP1 

78  250   CONTINUE 

79  CRLL  HELP  (6) 

80  GO  TO  210 

81  260  CONTINUE 

82  WRITE  (6,650) 

83  GO  TO  210 

84  270   CONTINUE 

85  I0LIST=N0INPT-1T0TIN 

86  WRITE  (6,600)  IOLIST 

87  IF  (IMODE.  NE.l)  1C0UNT=  ICOUNTM 

88  GO  TO  290 

89  280  CONTINUE 

90  IF  (NRMT.GE.2)  GO  TO  290 

91  WRITE  (6,550) 

92  GO  TO  230 

93  290  CONTINUE 

94  IF  (NRMT.GT.2)  NRMT=2 

95  1RERR=0 

96  DO  300  IC0L=1.NRMT 

97  CRLL  REASON  (1 . ICOL, IROW. RNTM  (1 . IROW) . RNTM (2, IROW) . IFRTRL) 

98  1RERR=>IRERR*1FRTRL 

99  300  CONTINUE 

100  IF  (1RERR.EQ.0)  GO  TO  320 
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101 

IF  (IM0DE.NE.1)  GO  TO  310 

102 

WRITE  (6.630)  IROW 

103 

GO  TO  220 

104 

310 

IC0UNT=1C0UNT+IRERR 

105 

320 

CONTINUE 

106 

IF  (IMOOE.EQ. 1)  WRITE  (6.560) 

107 

IF  (IMOOE.NE.l)  WRITE  (6.520) 

108 

DO  460  !ROW=l,NMEfiS 

109 

330 

CONTINUE 

110 

NRMT=>0 

111 

MRXINP=NHERS+1-IR0W 

112 

340 

CONTINUE 

113 

IF  (IMOOE.EQ. 1. AND. IROW. EQ.NMERS)  WRITE  (6.580)  IP 

114 

IF  (IMQDE.EQ.1.RND.IR0W.NE.NMERS)  WRITE  (6.570)  IP 

115 

1NMERS 

116 

350 

CONTINUE 

117 

IT0TIN=MRXINP-NRMT 

118 

1STRRT=1R0W+NRMT 

119 

CALL  INPT  12.IT0TIN.N0INPT. I VRL. TEMP (1STRRT) .IE0F) 

120 

IF  (1E0F.EQ.0.0R.IE0F.EQ.5)  NRMT=NRMT+NOINPT 

121 

IF  (IMOOE.EQ. 1)  GO  TO  360 

122 

1NLINE=INLINE-M 

123 

IF  (NRMT.EQ.MRXINP)  GO  TO  410 

12>1 

IF  (IEOF.EQ.O)  GO  TO  350 

125 

IF  (IEOF.EQ.l)  WRITE  (6.710) 

126 

IF  (IEOF.EQ.l)  GO  TO  470 

127 

IF  (IE0F.EQ.5)  GO  TO  390 

128 

IC0UNT=IC0UNT*1 

129 

WRITE  (6.700)  INLINE 

130 

GO  TO  460 

131 

360 

CONTINUE 

132 

IE0FP1=IE0F*1 

133 

GO  TO  (400.380.370,380.490.390) .IE0FP1 

134 

370 

CONTINUE 

135 

CALL  HELP  (7) 

136 

GO  TO  340 

137 

380 

CONTINUE 

138 

WRITE  (6,650) 

139 

GO  TO  340 

110 

390 

CONTINUE 

141 

I0LIST=N0INPT-IT0T1N 

142 

WRITE  (6.600)  IOLIST 

143 

IF  (IMOOE.NE.l)  IC0UNT=IC0UNT+1 

144 

NRMT=MRXINP 

145 

GO  TO  410 

146 

400 

CONTINUE 

147 

IF  (NOINPT.EQ.ITOTIN)  GO  TO  410 

148 

ITOTIN=MRXINP-NRMT 

149 

WRITE  (6.590)  IT0T1N.IR0W 

150 

GO  TO  350 

IROW.MRXINP.NMERS 
IROW.MRXINP. IROW, 
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151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 


•410 


CC 


902 
1120 

430 

440 


450 
460 
470 


480 

490 

500 
510 

520 
530 
540 
550 
560 


570 
580 
590 

600 

610 
620 
630 
640 


CONTINUE 

1LAST=IR0W+NAMT-1 

1RERR=0 

DO  420  1C0L=1R0W.ILAST 

WRITE  (6.902)  IROW. ICOL. ILAST.  1MGDE.  TEMP(ICOL) 

FORMAT  ('   ACC.  4110,  4F10.3) 

CALL  REASON  (2. ICOL. IROW. TEMP (ICOL) . 0. IFATAL) 

IRERR=1RERR+IFATAL 

WRITE  (6,901)  I  ROW, I  COL. IRERR.lFATAl 

CONTINUE 

IF  (1RERR.EQ.0)  GO  TO  440 

IF  (IM00E.EQ.2)  GO  TO  430 

WRITE  (6.640)  IROW 

GO  TO  330 

1CGUNT=IC0UNT*IRERR 

GO  TO  460 

CONTINUE 

DO  450  1C0L=1R0W.NMEAS 

N=NDX  (IROW.ICOL) 

CORR(N)  -TEMP  (ICOL) 

CONTINUE 

CONTINUE 

CONTINUE 

IF  (ICOUNT.EQ.O)  ISAMP=1 

IF     1IH0DE.NE.11     WRITE     (6,680)     IC0UNT 

GO  TO  490 

WRITE  (6.660) 

RETURN 

CONTINUE 

RETURN 

(5X.  'ENTER  PAIRS  OF  MEANS  AND  STANDARD  DEVIATIONS--') 
(//, '  DIAGNOSTICS  OF  FORMATTED  READ  ',/,'  OPER  DESC: 


FORMAT 

FORMAT 
l./l 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 
USE  UPPER 
2ENTS  MUST 

FORMAT  (5X.4HR0W  . 12. ' ( 
(5X.4HR0W  . 12. ' ( 
(10X,23H*k*ENTER 


10A4 


FORMAT 
FORMAT 

1) 
FORMAT 

1) 
FORMAT 
FORMAT 
FORMAT 
FORMAT 


(//.1X.  'CORRELATION/ANTHROPOMETRIC  MEASUREMENTS:  ') 

('   MEASUREMENT  -  '.12.'  ?---') 

('  MEASUREMENT  -',12) 

(5X, 'ENTER  STANDARO  DEVIATION  --') 

(5X. 'ENTER  ANTHROPOMETRIC  CORRELATION  MATRIX  DATA  IN  A  ROWW 
\/.5X.  'TRIANGULAR  FORM.  './, 5X.  'NOTE:  ALL  DIAGONAL  ELEM 
EQUAL  1.000') 

.12.'  ELEMENTS.  COLUMNS  ' , 12.  '- ' . 12.  ')  ') 

,12.'  ELEMENT.  COLUHN  '.12.')') 

THE  REMAINING  .12.18H  ELEMENTS  FOR  ROW  ,12 


(10X,  '*««ACCPT051  THE  LAST  ',12.'  INPUTS  HAVE  BEEN  IGNORED' 

(5X. 'INPUT  MODE  (1  =  INTERACT  I VE:2  =  F I LE)--') 

(SX.'READ  MODE (1=F0RMATTED;2=UNF0RMATTED) -- ')       \ 

(10X, '««*ACCPT060  RE-ENTER  DATA  FOR  ',12) 
(10X.  '**«ACCPT1300  RE-ENTER  DATA  FOR  ROW  '.13) 
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201 

650 

FORMAT 

202 

660 

FORMAT 

203 

670 

FORMAT 

201 

1RAM'./ 

205 

680 

FORMAT 

206 

690 

FORMAT 

207 

700 

FORMAT 

208 

710 

FORMAT 

209 

END 

UOX.  '"""ACCPT155  INVALID  RESPONSE') 

UOX. '"""ACCPT  UNEXPECTED  EOF  ON  NIMEAS  FOR  MEAS  '.12) 
UOX,  '"""REA0FLM030  FILE  ANTSTA  IS  NOT  ATTACHED  TO  THE  PROG 
'.10X.  '"""UNFORMATTED  READ  ON  THIS  FILE  IS  ABORTED  ---'//) 
UOX,  '"""ACCPT3000  TOTAL  ERRORS  ON  INPUT  DATA  '.15) 
('  MEAN.STD  DEV  ^CORRELATION  DATA  READ') 
(10X. '"""ACCPT150  INVALID  DATA  AT  LINE  '.15) 
(10X. '"""ACCPT1200  UNEXPECTED  END  OF  FILE') 


93 


1 

SUBROUTINE  EDIT 

2 

COMMON/MCSM/CORR (78) .  ANTM  (2. 121  ,V (25)  .NMEAS.KN, 

3 

1    SQROOT  (78) ,X  (25)  . 1U, NTMEflS, I  COUNT, 

il 

2    INAME  (10) , ITMEAS.N1MEAS.MAXSAM.NTEMP 

5 

COMMGN/MISC/IQIN.  IOOUT,  ICARO 

6 

DIMENSION  BUFF  (3) 

7 

OATA  BUFF/3«0./ 

8 

I0IN-5 

9 

100 

CONTINUE 

10 

HR1TE  (6.350) 

11 

CALL  1NPT  (l.l.NVALSZ.IGPT.XVAL.IEOF) 

12 

]EOFPl=IEOF+l 

13 

GOTO  (130.100,110.120,340.120).  IE0FP1 

111 

110 

CONTINUE 

15 

CALL  HELP  (8) 

16 

GO  TO  100 

17 

120 

CONTINUE 

18 

WRITE  (6.360) 

19 

GO  TO  100 

20 

130 

CONTINUE 

21 

IF  (10PT.LT.1.0R.I0PT.GT.5)  GO  TO  120 

22 

GOTO  (140. 140.240.330.340)  .IOPT 

23 

140 

CONTINUE 

2<J 

IC0L=10PT 

25 

150 

CONTINUE 

26 

MAXINP=2 

27 

NAMT=0 

28 

160 

CONTINUE 

29 

IF  (ICOL.EQ.l)  WRITE  (6.380) 

30 

IF  (IC0L.EQ.2)  WRITE  (6.390) 

31 

170 

CONTINUE 

32 

IVALSZ=MAXINP-NAMT 

33 

J=NAMT+1 

34 

CALL  INPT  (2. IVALSZ.NVALSZ. IVAL.BUFF (J) , 1E0F) 

35 

IE0FP1-IE0FP1 

36 

GOTO  (220.180.200.190.100.210).  1E0FP1 

37 

180 

CONTINUE 

38 

WRITE  (6.370) 

39 

GO  TO  160 

40 

190 

CONTINUE 

HI 

WRITE  (6,360) 

42 

GO  TO  160 

43 

200 

CONTINUE 

VJU 

IF  (ICOL.EQ. 1)  CALL  HELP  (11) 

45 

IF  (ICOL.EQ. 2)  CALL  HELP  (12) 

<16 

GO  TO  160 

•17 

210 

CONTINUE 

48 

I0L1ST=NVALSZ-IVALSZ 

49 

WRITE  (6.400)  10LIST 

50 

220 

CONTINUE 
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51 

1R0W=BUFF  (1) 

52 

NAMT=NVALSZ+NAMT 

53 

IF  (NAMT.GE.MAXINP)  GO  TO  230 

54 

IF  (1C0L.EQ.1)  WRITE  (6.1420)  IROW 

55 

IF  (IC0L.EQ.2)  WRITE  (6. 430)  IROW 

56 

GO  TO  170 

57 

230 

CONTINUE 

58 

CALL  REASON  (1 , 1C0L. IROW, BUFF  (2)  , 0,  IFATAL) 

59 

IF  (IFATAL.EQ.l)  GO  TO  150 

60 

ANTMUCOL,  IROW)  =  BUFF  (2) 

61 

GO  TO  100 

62 

240 

CONTINUE 

63 

MAX1NP-3 

64 

NAMT-0 

65 

250 

CONTINUE 

66 

WRITE  (6.1I0) 

67 

260 

CONTINUE 

68 

IVALSZ-MAX1NP+NAMT 

69 

J»NAHTM 

70 

CALL  INPT  (2.IVALSZ.NVALSZ.IVAL.BUFFU)  .IEOF) 

71 

IEOFP1=IE0F-M 

72 

GOTO  (310.270.290.280.100.300).  IE0FP1 

73 

270 

CONTINUE 

74 

WRITE  (6,370) 

75 

GO  TO  250 

76 

280 

CONTINUE 

77 

WRITE  (6.360) 

78 

GO  TO  250 

79 

290 

CONTINUE 

80 

CALL  HELP  (13) 

81 

GO  TO  250 

82 

300 

CONTINUE 

83 

I0LI5T=NVALSZ-IVALSZ 

8«4 

WRITE  (6.1400)  IOLIST 

85 

GO  TO  320 

86 

310 

CONTINUE 

87 

ICOL  =  BUFF  (1) 

88 

1R0W  =  BUFF(2) 

89 

NAMT»NAHT*NVALSZ 

90 

IF  (NAMT.GE.MAXINP)  GO  TO  320 

91 

IF  (NAMT.EQ.l)  WRITE  (6.410)  1C0L 

92 

IF  (NAMT.EQ.2)  WRITE  (6.450)  1C0L.1R0W 

93 

GO  TO  260 

9>4 

320 

CONTINUE 

95 

CALL  REASON  (3. ICOL. IROW. BUFF (3) . 0. IFATAL) 

96 

IF  (IFATAL.EQ.l)  GO  TO  240 

97 

CALL  REASON  (2. ICOL. IROW. BUFF  (3) . 0. IFATAL) 

98 

IF  (IFATAL.EQ.l)  GO  TO  240 

99 

N=NOX  (IROW. ICOL) 

100 

CORR(N)  =BUFF  (3) 

95 


101 

GO  TO  100 

102 

330 

CONTINUE 

103 

CALL  PRINT 

104 

GO  TO  100 

105 

340 

CONTINUE 

106 

RETURN 

107 

350 

FORMAT 

(5X, 

108 

360 

FORMAT 

(10X 

109 

370 

FORMAT 

(10X 

110 

380 

FORMAT 

(5X. 

111 

390 

FORHAT 

(5X. 

112 

400 

FORMAT 

(10X 

113 

410 

FORMAT 

<5X, 

im 

420 

FORMAT 

(5X, 

115 

430 

FORMAT 

(5X, 

116 

440 

FORMAT 

(5X. 

117 

450 

FORMAT 

(5X. 

118 

END 

'EDIT  MODE  (1-MEANS;2-ST0:3-ANTH.C0RR:4-PR1NT:5-END)  — ') 

, '«««EDIT044  INVALID  RESPONSE") 

.  '«*«E0IT095  NO  INFORMATION  ENCOUNTERED') 

'INPUT  INDEX  ANO  MEAN  VALUE') 

'INPUT  INDEX  AND  STANDARD  DEVIATION') 

.  '«««EDIT  042  THE  LAST'. 12.'  INPUTS  HAVE  BEEN  IGNORED') 

'INPUT    COLUMN  .  ROW    AND  CORRELATION') 

'INPUT  MEAN  VALUE  FOR  INDEX  '.13) 

'INPUT  STANDARD  DEVIATION  FOR  INDEX  ',13) 

'INPUT  ROW  AND  CORRELATION  FOR  COLUMN  '.12) 

'INPUT  CORRELATION  FOR  RON  '.12,'  COLUMN  ',12) 
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1 

2 

3 

U 

5 

6 

7 

8 

9 

10 

11 

12 

100 

13 

110 

14 

15 

16 

17 

18 

19 

20 

120 

21 

22 

23 

130 

24 

25 

26 

140 

27 

28 

29 

150 

30 

31 

32 

160 

33 

170 

34 

180 

35 

190 

36 

SUBROUTINE  NAMEFL  (IERROR) 

CQMMGN/HCSM/CORR (78)  .ANTM (2.  12)  .  V (25)  .NMEAS.KN, 

1  SQROOT  (78) ,X (25) . 1U.NTMEAS. IC0UN7. 

2  I  NAME  (10) . ITMEAS.NIMEAS.MAXSAM.NTEMP 
CQMMON/MISC/IOIN. IOOUT.  ICARD 
DIMENSION  IBUFF  (10) .XVfiL  (10) 

DATA  IBUFF/10*i4H    / 

DATA  I4BLK/4H    / 

IERROR=0 

00  100  J-1,4 

IBUFF  U)=I4BLK 

CONTINUE 

CONTINUE 

IF  (10IN.EQ.5)  WRITE  (6.180) 

CALL  1NPT  (4.10.NVALSZ.  IBUFF,  XVAL.  IEOF) 

IF  (IEOF. LT. 2. OR. IEOF. EQ. 5)  GO  TO  150 

IF  (I0IN.NE.5)  GO  TO  140 

1E0FM1=IE0F-1 

GO  TO  (120,130.170) .  IE0FM1 

CONTINUE 

CALL  HELP  (14) 

GO  TO  110 

CONTINUE 

WRITE  (6.190) 

GO  TO  110 

CONTINUE 

WRITE  (6.190) 

IERR0R=1 

CONTINUE 

DO  160  J-l, 10 

INAME  (J) =IBUFF  (J) 

CONTINUE 

RETURN 

FORMAT  (5X.  'OPERATOR  SAMPLE  DESCRIPTION  (40  CHAR  MAX) ') 

FORMAT  (10X.  '>oo<NAMEFL027  INVALID  RESPONSE') 
END 
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0  SUBROUTINE  OPSAMP  USAMP) 

1  CGMMON/MCSM/CORR  (78) .RNTM  (2. 12) .V  (25) .NMERS.KN. 

2  1    SQROOT  (78) ,X  (25)  . IU.NTMERS. ICOUNT, 

3  2    INRME  (10) . ITMEAS.NIMEAS.MAXSAM.NTEMP 

4  C0MM0N/MISC/10IN.I00UT.ICRRD 

5  DIMENSION  IRAN  (9) 

6  ORTR  IRAN  (1)  .IRAN  (2)  ,  IRAN  (3)  .IRAN  (4)  .  IRAN  (5)  .  IRAN  (6)  .IRAN  (7)  . 

7  1IRAN  (8) .IRAN  (9) /33333333, 55555555. 77777777. 222222221 , 

8  1  444444443.  666666669.  888888889,  999999999.  123456789/ 

9  I0IN»5 

10  DO  100  K=1,KN 

11  SQROOT  (K)=CORR(K) 

12  100   CONTINUE 

13  CALL  SQRMX  (NERROR) 

14  IF  (NERROR. EQ.l)  RETURN 

15  110   CONTINUE 

16  WRITE  (6.320) 

17  CALL  1NPT  (1. l.NVALSZ.NSRMP.XVAL. IEOF) 

18  IE0FP1-1E0F-M 

19  GOTO  (150, 120. 130. mO. 290. 140)  .IEOFPl 

20  120   CONTINUE 

21  NRITE  (6.330) 

22  GO  TO  110 

23  130   CONTINUE 

21  CALL  HELP  (15) 

25  GO  TO  110 

26  140   CONTINUE 

27  WRITE  (6.340) 

28  GO  TO  110 

29  150  CONTINUE 

30  IF  (NSRMP.GT.O.RND.NSAMP.LE.MAXSRM)  GO  TO  160 

31  WRITE     (6.350)    MflXSAM 

32  GO   TO    110 

33  160   CONTINUE 

3>4  WRITE  (6.300) 

35  CALL  INPT  (1 . 1 . NVALSZ. MTTPE, XVAL. IEOF) 

36  IE0FP1-IE0F+1 

37  GOTO  (190,160.170.180.290. 180) .IEOFPl 

38  170   CONTINUE 

39  CALL  HELP  (20) 

40  GO  TO  160 

41  180   CONTINUE 

42  WRITE  (6.340) 

43  GO  TO  160 

44  190   CONTINUE 

145  IF  (MTTPE. NE. 32. RND. MTTPE. NE. 60)  GO  TO  180 

MS  200  CONTINUE 

17  WRITE  (6,310) 

48  CALL  INPT  (1,1. NVALSZ. IS. XVAL. IEOF) 

49  1E0FP1=IE0F*1 


98 


50 
51 
52 
53 
51 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
7« 
75 
76 
77 
78 
79 
80 
81 
82 
83 
81 
85 
86 
87 
88 
89 
90 
91 
92 
93 


210 
220 
230 
240 

250 
260 


270 
280 


370 

290 

300 
310 
320 
330 
340 
350 
360 


GOTO  (240.210.220.230.290.230).  IEOFPl 

CONTINUE 

WRITE  (6,330) 

GO  TO  200 

CONTINUE 

CALL  HELP  (21) 

GO  TO  200 

CONTINUE 

WRITE  (6.340) 

GO  TO  200 

CONTINUE 

IF  (1S.EQ.0)  GO  TO  250 

IF  (IS.LT.l.OR. 1S.GT.9)  GO  TO  230 

IU=IRAN(IS) 

GO  TO  260 

CONTINUE 

IU=671 1 1 133 

CONTINUE 

REWIND  NTEMP 

WRITE  (NTEMP)  INAME.NMEflS 

DO  280  I=l.NSflMP 

CALL  S1MSUB  (MTYPE) 

DO  270  K=1.NMEAS 

V  (K)  =V  (K)  «ANTM  (2.K)  +  ANTM  (l.K) 

CONTINUE 

WRITE  (NTEMP)    (V  (10)  .  10=  1 .  NMEAS) 

CONTINUE 

1SAMP»1 

END  FILE  NTEMP 

REWIND  NTEMP 

REWIND  NTMEflS 

WRITE  (6,370) 

FORMAT  (  10X,  '**  SAMPLE  DATA  SAVED  ON  FILE  9') 

WRITE  (6.360) 

CONTINUE 

RETURN 


(32  OR  60  BIT)  — •] 
IF  NONE  OR  1-9)--') 


FORMAT  (5X,  "ENTER  MACHINE  WORD  SIZE 

FORMAT  (5X.  'ENTER  RANDOM  NO  SEED  (0 

FORMAT  (5X, 'INPUT  SAMPLE  SIZE  --') 

FORMAT  (10X, '***0PSAMP215  NO  INFORMATION  ENCOUNTERED') 

FORMAT  (10X,  '***QPSAMP411  INVALID  RESPONSE ') 

FORMAT  (10X. ****0PSAMP009  SAMPLE  SIZE  MUST  BE  BETWEEN  1  AND  ' 

FORMAT  (10X.  'mmxSAMPLE  DATA  SAVED  ON  FILE  SAMPLE  OR  TAPE  3.') 

END 


13) 
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1  SUBROUTINE  PRINT 

2  COMMON/MCSM/CORR  (78) .RNTM (2. 12) ,V (25) .NMERS.KN. 

3  1    SQROOT  (78) ,X(25) . IU.NTMERS. ICOUNT. 

4  2    INRME  (10) . ITMERS.NlMERS.MfiXSRM.NTEMP 

5  CQMMON/MISC/IQIN.IGOUT.ICRRD 

6  INTEGER  FMT.FMT1 

7  OIMENSION  FMT  (6)  .NUMU2)  ,TEMP(12)  ,INDEX(2)  .FMTl  (5) 

8  DRTR  FMTl  /'(1X,\  '     \  '(2X,\  'R2,2\     'X))V 

9  DRTR  FMT  /MIX.*,  '     ',  '(6X)',  '     '.  'F6.3'.')V 

10  DRTR  MUM/1H1.1H2, 1H3, 1H4, 1H5, 1H6, 1H7.1H8. 1H9.2H10.2H1 1 ,2H12/ 

11  FMTl  (2) =NUM  (NMERS) 

12  ISTORE=0 

13  IOFLRG=0 
H4  100   CONTINUE 

15  WRITE  (6,630) 

16  CALL  INPT  (l.l.NVRLSZ.IOPT.XVRL.IEOF) 

17  IEQFP1=1E0F+1 

18  GOTO  (120.110.110. 110,390. 110) . IE0FP1 

19  110   CONTINUE 

20  HRITE  (6.650) 

21  GO  TO  100 

22  120   CONTINUE 

23  IF  (I0PT.LT.1.0R.10PT.GT.41G0  TO  110 
24.           IF  (10PT.EQ.1)  GO  TO  130 

25  IF  (I0PT.EQ.2)  GO  TO  150 

26  IF  (10PT.EQ.4)  GO  TO  390 

27  GO  TO  190 

28  130   CONTINUE 

29  HRITE  (6.580) 

30  DO  140  IR0W=1, NMERS 

31  WRITE  (6.640)  IROW.  (RNTM  (JO. IROW) . J0=1 . 2) 

32  110   CONTINUE 

33  WRITE  (6.620) 

34  GO  TO  100 

35  150   CONTINUE 

36  WRITE  (6.590) 

37  WRITE  (6, FMTl)  (NUM  (10) . 10=1 .NMERSl 

38  WRITE  (6.600) 

39  DO    180    IR0W=1.NMERS 

40  DO    160    IC0L=1.NMERS 
m  N=NDX  (IR0W.1C0U 

1*2  TEMP(ICOL)  =CORR(N) 

143  160      CONTINUE 

<4>4  IF     (IROW.GT.l)    GO   TO    170 

145  WRITE  (6,610)  (TEMP(ICOL)  ,  IC0L=1R0W.  NMERS) 

46  GO  TO  180 

147  170   FMT  (2)»NUM(1R0W-1) 

148  FMT  (4)  »NUM(NMERS*1-IR0W) 

49  WRITE  (6. FMT)  (TEMP  (ICOL) . ICOL= IROW, NMERS) 

50  180   CONTINUE 
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St 

WRITE  (6.620) 

52 

GO  TO  100 

53 

190 

CONTINUE 

54 

WRITE  (6,1430) 

55 

IF  UQFLfiG.LT.  1)  WRITE  (6.440) 

56 

I0FLfiG=l 

57 

CALL  1NPT  (1.2.NVALSZ.INDEX.XVAL.1EQF) 

58 

IE0FP1=1E0F+1 

59 

GOTO  (230.270.200.210.100.220).  IE0FP1 

60 

200 

CONTINUE 

61 

CALL  HELP  (19) 

62 

10FLAG=0 

63 

GO  TO  190 

64 

210 

CONTINUE 

65 

WRITE  (6.420) 

66 

GO  TO  190 

67 

220 

CONTINUE 

68 

I0LlST=NVflLSZ-2 

69 

WRITE  (6.450)  I0L1ST 

70 

230 

CONTINUE 

71 

IF  (INDEX  (1)  .GT.O)  GO  TO  240 

72 

WRITE  (6.480) 

73 

GO  TO  190 

71 

240 

CONTINUE 

75 

IF  (I  STORE. GT.O. AND. INDEX  (1) .GT. I  STORE) 

76 

NRSKIP=IN0EX  (1) -1 

77 

GO  TO  260 

78 

250 

CONTINUE 

79 

NRSKIP=INDEX(1) -ISTORE-1 

80 

260 

CONTINUE 

81 

ISTRRT=>INDEX  (1) 

82 

IF  (NVALSZ.EO.l)  GO  TO  280 

83 

IF  (NVALSZ.GE.2)  GO  TO  290 

8>4 

270 

CONTINUE 

85 

ISTORE=0 

86 

ISTART=1 

87 

IENQ-999 

88 

NRSKIP=0 

89 

GO  TO  320 

90 

280 

CONTINUE 

91 

IENO=ISTART 

92 

GO  TO  310 

93 

290 

CONTINUE 

94 

IF  (INDEX  (1)  .LE. INDEX  (2) I  GO  TO  300 

95 

WRITE  (6.490) 

96 

GO  TO  190 

97 

300 

CONTINUE 

98 

IEND=INDEX(2) 

99 

310 

CONTINUE 

100 

IF  (INDEX  (1)  .GT. I  STORE. AND. 1  STORE. NE.O) 

GO   TO   250 


GO    TO    330 


101 


101 

320 

102 

103 

101 

330 

105 

106 

107 

108 

340 

109 

350 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

12*4 

125 

126 

360 

127 

128 

370 

129 

130 

131 

132 

380 

133 

134 

135 

390 

136 

137 

138 

•400 

139 

mo 

mi 

142 

410 

143 

mi 

145 

■120 

146 

430 

147 

440 

148 

149 

150 

CONTINUE 

BEHIND  NTEMP 

READ  (NTEMP, END=410)  INAME.NMEAS 

CONTINUE 

IF  (NRSKIP.LT.1)  GO  TO  350 

DO  340  I=1,NRSK1P 

READ  (NTEMP, END=400)  (V  (K)  , K=l .NMEAS) 

CONTINUE 

CONTINUE 

WRITE  (6,530)  (INAME  (10)  .  10=1 ,  10) 

IF  (1START.EQ.IEN0)  WRITE  (6,540)  IENO 

IF  (ISTRRT.NE. IENO. AND. IEND. NE. 999)  WRITE  (6.550)  ISTART.1END 

WRITE  (6.500) 

WRITE  (6.FMT1)  (NUH  (10)  ,  10=  1  .NMEAS) 

WRITE  (6.600) 

ISKIP=0 

NF0UND=0 

DO  360  1=1STRRT.IEND 

IF     IISKIP.EQ. II    WRITE    (6.460) 

1SKIP-0 

IF  (ISTRRT.NE. IENO. AND. l.NE.IEND.RNO.  (1/5) *5. EQ. I)  ISKIP=1 

RERD  (NTEMP, END=370)  (V  (K)  , K=l , NMEflS) 

NFQUND=NF0UND-"1 

ISTORE=I 

WRITE  (6,510)  (V  (K)  ,K=1, NMEAS) 

CONTINUE 

GO  TO  380 

CONTINUE 

IF  (IEND. NE. 999. AND. ISTRRT.NE. IEND)  WRITE  (6.470)  1END 

1ST0RE-1 

IF  (NFOUND.LT.l)  WRITE  (6.560) 

CONTINUE 

WRITE  (6,520) 

GO  TO  190 

CONTINUE 

REWINO  NTEMP 

RETURN 

CONTINUE 

I0LIST=I-2+IST0RE 

WRITE  (6,560) 

GO  TO  190 

CONTINUE 

WRITE  (6,570) 

GO  TO  100 

FORMAT  (10X,  **><kPRINT135  INVALID  RESPONSE') 

FORMAT  (5X, 'OUTPUT  MODE:--') 

FORMAT  (lOX.'FOR  RLL  OPERATOR  DATA  (ENTER  CARRIAGE  RETURN)  ',/. 10X. 
1'FOR  INDIVIDUAL  OPERATOR  DATA  (ENTER  OPERATOR  NUMBER)  ',/, 10X,  'FOR 
2A  SEQUENCE  OF  OPERATORS  (ENTER  RANGE  OF  OPERATOR  NUMBERS)  ',/, 10X,  ' 
3EN0  (ENTER  $)  ') 
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151 
152 
153 
ISM 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
16? 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 


150 

460 
470 
480 

490 

500 
510 
520 
530 
540 
550 
560 

570 

580 

590 

600 
610 
620 
630 

640 
650 


FORMAT 
FORMAT 
FORMAT 
FORMAT 

10  ZERO 
FORMAT 

IX') 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 


(10X.  '***PRINT119  THE  LAST '. 14,  '  INPUTS  HAVE  BEEN  IGNORED') 

(1H  ) 

UOX.  '***PRINT211  ERROR  IN  UPPER  INDEX  --'.141 

(10X. '***PRINT120  INDEX  NO.  LESS  THAN  OR  EQUAL ',/. ****T 

IS  NOT  VALID*) 

(10X, '***PRINT160  FIRST  INDEX  MUST  BE  LESS  THAN  SECOND  INDE 


--'.14) 

NUMBER  EXCEEDS 


RECORDS  IN  OPERATOR 


(27X,  'OPERATOR  SAMPLE'./, 71  (1 H-J ) 
(11F6.2.F5.2) 
(71  (1H-)  ) 
(1X.10A4. 10X.A4) 
UOX.  'OPERATOR  NO*.  14) 
(10X,  'OPERATOR  NOS',14,  ' 
UOX,  '***PR1NT1500  INDEX 
1  FILE***') 
FORMAT  UOX. '***PRINTl600  FILE  SAMPLE  IS  NOT  ATTACHED  TO  THE  PROGR 
1AM',/,10X. '***OPERATOR  OUTPUT  REQUEST  ABORTEO  ---',//) 

FORMAT  1/.21X,  'MEANS  AND  STANDARD  DEVIATIONS'./, IX, 70 tlH-1 ,/,2X, 'M 
1EASUREMENT\6X.  'MEANS ',6X,  'STD  DEV  ' , /, 1 X.71  (1H-) ) 

FORMAT  (/,7X, 'CORRELATION  MATRIX  FOR  ANTHROPOMETRIC  MEASUREMENTS'. 
1/.1X.70UH-)  ) 
FORMAT  (1X.70UH-)) 
(IX,  12F6.3) 
(1X,70(1H-)  ) 
(5X, 'PRINT  MOOE (1-MEAN/STD:2-C0RR  MATRX;  ' .  '3-OPER  MEAS:4-EN 


FORMAT 
FORMAT 
FORMAT 
ID)  —  ') 
FORMAT 
FORMAT 
END 


(7X,I2.5X,2(5X,F7.2)) 

(10X, '***PR1NT103  INVALID  RESPONSE') 
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1  SUBROUTINE  RRND2  (RNORM. MTTPE. QURN, 1J) 

2  C0MM0N/MCSM/C0RR(78)  . RNTM(2.12) ,  V  (25) .NMERS.KN, 

3  1    SQROOT  (78)  ,X(25)  , IU.NTMERS, ICOUNT, 

4  2    1NRME  (10)  ,  ITMERS.NIMERS.MRXSRM.NTEMP 

5  DIMENSION  QURN  (2) 

6  COMMON/MISC/IOIN.IOOUT. ICRRD 

7  1J=IJ-»1 

8  IF  (1J.EQ.2)  GO  TO  110 

9  100   CONTINUE 

10  X1=2.0*UNF0RM(MTTPE) -1.0 

11  Y1=2.0*UNF0RM  (MTTPE) -1.0 

12  S=X1«X1+Y1«Y1 

13  IF  (S.GE.1.0)  GO  TO  100 

14  S=SORT  (-2.0*RL0G(S) /S) 

15  QURNU)  =XI«S 

16  QURN(2)=T1*S 

17  110   CONTINUE 

18  RN0RM=QURN(IJ) 

19  IF  (IJ.EQ.2)  IJ=0 

20  RETURN 

21  END 
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0  SUBROUTINE  REASON  (ICALL. ICOL. IROW. TSTOR1 . TSTQR2. IFATAL) 

1  COMMON/MCSM/CORR (78) ,ANTM(2. 12) .V (25) .NMEAS.KN. 

2  1    SQR00T(78) ,X(25) . IU. NTMEAS. ICOUNT. 

3  2    1NAME  (10) . ITMEAS.N1MEAS.MAXSAM.NTEMP 

4  ,  COMMON/MISC/IOIN.IOGUT.ICARO 

5  IFflTRL=0 

6  IAFLAG-0 

7  C  NR1TE16.901)  ICALL. ICOL, IROW. IFATAL. TST0R1 .TST0R2 

8  901  FORMAT  ['   REASON  \  MHO.  4F10.3) 

9  IF  (ICALL.LT.O)  IAFLAG=1 

10  ICALL-IA8S  (ICALL) 

11  IF  (1CALL.GT.1)  GO  TO  110 

12  IHRITE-1 

13  IF  (IROW.LT.l.OR.IRGW.GT.NMEAS)  GO  TO  120 

14  IF  (IC0L.EQ.2)  GO  TO  100 

15  RETURN 

16  100   CONTINUE 

17  IWRITE-2 

18  IF  (TST0R2.LT.0.)  GO  TO  120 

19  IWR1TE-3 

20  IF  (TST0R1.LE.TST0R2)  GO  TO  120 

21  RETURN 

22  110   CONTINUE 

23  IHR1TE-1 

24  IF  (IROW. LE.O. OR. IROW. GT.NMEAS)  GO  TO  130 

25  1WRITE-2 

26  IF  (ICOL. LE.O. OR. ICOL. GT.NMEAS)  GO  TO  130 

27  1WR1TE*3 

28  IF  (1ROW.GT.IC0L)  GO  TO  130 

29  IF  (IC0L.EQ.3)  RETURN 

30  1WR1TE-4 

31  IF  (1C0L.EQ. IROW. AND. TSTORl.NE. 1.000)  GO  TO  130 

32  IWRITE*5 

33  IF  (TSTGR1.LT. -1.000)  GO  TO  130 
3U  IWR1TE-6 

35  IF  (TST0R1.GT. 1.000)  GO  TO  130 

36  RETURN 

37  120  CONTINUE 

38  IF  (IHRITE.NE.1.AND.10IN.NE.5)  WRITE  (6.280)  IROH 

39  GO  TO  (190.  mO. 150)  ,  IHRITE 
•10  130   CONTINUE 

VI  GOTO  (190,200.210.160.170.180).  IHRITE 

42  140  CONTINUE 

43  WRITE  (6.230) 
U4  GO  TO  220 

45  150   CONTINUE 

46  WRITE  (6.240) 

47  GO  TO  220 

48  160   CONTINUE 

49  WRITE  (6,250)  IR0W.1C0L 
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50  GO  TO  220 

51  170   CONTINUE 

52  WRITE  (6.260)  1R0W.IC0L 

53  GO  TO  220 
51  180   CONTINUE 

55  HR1TE  (6,270)  1R0W.ICGL 

56  GO  TO  220 

57  190   CONTINUE 

58  WRITE  (6.290)  1R0E.NMEAS 

59  GO  TO  220 

60  200   CONTINUE 

61  WRITE  (6.300)  ICOL.NMEAS 

62  GO  TO  220 

63  210   CONTINUE 

61  WRITE  (6.310)  IROW.ICOL 

65  GO  TO  220 

66  220   CONTINUE 

67  IFATAL=1 

68  RETURN 

69  230   FORMAT  (10X.  '*««REAS0N1 150  STANOARD  DEVIATIONS  LESS  THAN  ZERO') 

70  210   FORMAT  (10X. '*mhREAS0N1200  5T0.  GREATER  THAN  MEAN  VALUE.') 

71  250   FORMAT  (10X.  '««»<REASaNl250  DIAGONAL  ELEMENT  IN  ROW  ',  12./.10X,  'wkm 

72  1C0LUMN  '.12.'  IS  NOT  EQUAL  TO  1.0*1 

73  260   FORMAT  (10X. '«h*REASQN1300  ENTRT  IN  COLUMN  ',12.'  ROW  \12,/.10X, ' 
7«1  1*h«IS  LESS  THAN  -1.0') 

75  270   FORMAT  (10X. '«««REAS0N1350  ENTRT  IN  COLUMN  '.12.'  ROW  '.I2./.10X.' 

76  1**«IS  GREATER  THAN  ♦1.0*1 

77  280   FORMAT  (8X.I2) 

78  290   FORMAT  (10X.  '««*REAS0N1500  ROW  '.12.'  IS  LESS  THAN  ZERO  OR  GREATER 

79  1  THAN  ',12) 

80  300   FORMAT  (10X. *«««REAS0N1 600  COLUMN  '.12.'  IS  LESS  THAN  ZERO  OR  GREA 

81  ITER  THAN  '.121 

82  310   FORMAT  (10X.  '*k*REASQN1700  ROW  ',12. 'COLUMN  '.12.'  IS  OUTSIDE  OF 

83  1UPPER',/. 10X,  '"""TRIANGULAR  REGION  OF  CORRELATION  MATRIX.*) 
81  END 
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1  SUBROUTINE  SAVE 

2  COMMON/MCSH/CORR  (78) ,RNTM(2. 12)  .V  (25) . NMERS. KN. 

3  1    5QR00T  (78)  ,X(25) . IU.NTMERS. ICOUNT, 

14  2    INRME  (10)  .  ITMERS.NIMERS.MfiXSflM.NTEMP 

5  COMMON/MISC/IOIN.IOGUT.ICRRO 

6  REMIND  NTMERS 

7  WRITE  (NTMERS)  NMERS. KN. INRME 

8  WRITE  (NTMERS)  (  (RNTM  (10, JO)  , 10=1 . 2)  . J0= 1 . NMERS) 

9  WRITE  (NTMERS)   (CORR (10) . 10= 1 . KN) 

10  REWIND  NTMERS 

11  WRITE  (6,100) 

12  RETURN 

13  100  FORMAT  ('  MERSEUREMENT  DRTR  SAVED  ON  FILE  ANTSTfl  OR  TAPE  2') 

14  END 
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1  SUBROUTINE  SIMSUB  (MTYPE) 

2  C0MM0N/MCSM/CQRR(78) .RNTM (2. 12) ,V (25)  .NMERS.KN. 

3  1    SQROOT  (78) ,X  (25)  . 1U. NTMERS, 1COUNT. 

■I  2    1NRME  (10) , ITMERS.NIMERS.MRXSAM.NTEMP 

5  C0MM0N/MISC/I01N.I00UT.ICRRD 

6  DIMENSION  T (45) 

7  IJ»0 

8  DO  100  I=1.NMERS 

9  V(I)=0.0 

10  CALL  RRND2  (RNORM. MTYPE. QURN, I  J) 

11  100   Y(1)=RN0RM 

12  DO  110  I=1.NMERS 

13  DO  1 10  J=l. I 

14  K  =  NDX(J,1) 

15  110   Vm-V(I)+Y  (J)*SQROOT  (K) 

16  RETURN 

17  END 
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0  SUBROUTINE  SORMX  (NERROR) 

1  COMMON/MCSM/CORR (78) ,ANTM (2. 12) , V (25) , NMEflS. KN. 

2  1   SQROOT  (78) ,X(25) , IU.NTMEflS. ICOUNT. 

3  2    INAME  (10) .ITMEflS.NIMEAS.MflXSflM.NTEMP 

4  COMMON/MISC/IOIN.IOOUT.ICARO 

5  DIMENSION  ZX  (7) 

6  INTEGER  ZX 

7  NERROR=0 

8  00  120  1=1.NMEAS 

9  zxm)=i 

10  N=NDX(I,1) 

11  IF  (SOROOT  (N) .LT.O)  GO  TO  130 

12  SO=SQRT  (SQROOT (N) ) 

13  00  100  J=I,NMEAS 
1W  N  =  NDX(I,J) 

15  SQROOT  (N) =SQROOT  (N)/SQ 

16  100   CONTINUE 

17  IF  (l.EQ.NMEAS)  NERROR=0 

18  IF  (l.EQ.NMEAS)  RETURN 

19  11=1-1 

20  ZX  (5)  -1 1 

21  00  110  12=11. NMEAS 

22  DO  110  J2=12. NMEflS 

23  N=NDX(I2.J2) 

24  L-NDXU.I2) 

25  K=NDX(I,J2) 

26  SQROOT  (N) =SQROOT (N) -SQROOT (L) *SQROOT (K) 

27  ZX(6)=I2 

28  ZX(7)=J2 

29  110   CONTINUE 

30  ZX(1)=N 

31  ZX(2)-L 

32  ZX(3)=K 

33  120   CONTINUE 

34  130   CONTINUE 

35  NERR0R=1 

36  HRITE  (6.180) 

37  1U0   CONTINUE 

38  HRITE  (6,190) 

39  CALL  INPT  (4. 1 . NVflLSZ. IflNS. XVflL. IEOF) 
"40  1EOFP1=IEOF+1 

41  GOTO  (170.160.150.170.160.170).  IE0FP1 

42  150   CONTINUE 

43  WRITE  (6.200) 

44  HRITE  (6.210)   (ZX  (I)  .  1  =  1 .  7)  .  SQROOT  (N) 

45  160   CONTINUE 
48  RETURN 

47  170   CONTINUE 

48  WRITE  (6,220) 

49  GO  TO  140 
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50 

180 

51 

52 

190 

53 

5>4 

200 

55 

56 

57 

210 

58 

220 

59 

FORMAT  (10X.  '««*SQRMX051  ANTHROPOMETRIC  CORRELATION  VALUES  ARE  INV 
1AL10') 

FORMAT  dOX,'--  ENTER  CARRIAGE  RETURN  TO  CONTINUE  PROCESSING  \  /. 10 
IX,*--  ENTER   ?   TO  DISPLAT  SORMX  VARIABLES.') 

FORMAT  (10X,  '*«*SQRMX062  AN  INDEFINITE  QUANTITY  HAS  BEEN  DETECTED 
1BT\/.  10X.  'kxwSQRMX.  VARIABLE  LIST  WITHIN  SQRMX:  ',  /.  10X '  N  L 
2     K     I     11    12    J2    SQRGOT',/. 10X.4H  (1H-) ) 

FORMAT  (8X.14  (3X,  12)  ,4X,3  (I2.3X)  .F6.3./.  10X.144  (1H-)  ,//) 

FORMAT  (10X.  '««*SQRMX060  INVALID  RESPONSE**** ') 

ENO 
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1  FUNCTION  NDX  (I.J) 

2  IF  U.GT.J)  GO  TO  100 

3  NDX=1+  UxJ-J)  /2 
U  RETURN 

5  100   NDX-J*  (1*1-11/2 

6  RETURN 

7  END 
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1  FUNCTION  UNFORM  (MTTPE) 

2  COMMON  /MCSM/  CORR  (78) . RNTM (2. 12) , V (25) .NMERS. KN. 

3  1SQROOT  (78) ,X  (25) . 1U.NTMERS. ICOUNT. 

V  2INRME  (10) .ITMERS.NlMERS.MRXSflM.NTEMP 

5  COMMON/MISC/IOIN.IOOUT.ICRRO 

6  DRTfl    !X/167772l3/.C/28m74976710655./ 

7  IF    (MTTPE. EQ. 32)    GO   TO    100 

8  1U=IU«1X 

9  UNF0RM=1U/C 

10  RETURN 

11  100      1U=IU*65539 

12  IF    (IU.GE.O)    GO   TO    110 

13  IU-1U+21171J836147-M 
IV                 110      UNF0RM-1U 

15  UNF0RM=UNF0RM«.4656613E-09 

16  RETURN 

17  END 
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1  SUBROUTINE  HELP  (I) 

2  RETURN 

3  END 
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99999 

100000 

100001 

100002 

100003 

100004 

100005 

100006 

100007 

100008 

100009 

100010 

100011 


1 


BLOCK  DATA 

COHMON/MCSM/CORR  (78) , ANTM (2. 12) .V (25) .NMERS.KN. 

1  SQROOT  (78) ,X (25)  . IU.NTMERS. ICOUNT. 

2  1NRME  (10)  . ITMERS.NIMERS.MRXSRM.NTEMP 
C0MM0N/M1SC/10IN. 100UT. 1CRRD 

DRTR  I0lN,I00UT,NTEMP.NTMEflS.NlMEflS/5.6,3.4.7/ 
DRTfl  MRXSRM/400/ 
DRTR   NTMERS.NIMERS/2.4/ 
l.ICOUNT/0/ 
DRTR  RNTM/2i4*0./,C0fiR/78«0./.  lNRME/U«yH     / 
END 


1 
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0  SUBROUTINE  PUNRT 

1  COMMON/MCSM/CORR (78) ,ONTM  (2. 12) . V (25) . NMEAS.  KN. 

2  1   SQRG0T(78) ,X(25) , IU.NTMEAS. ICOUNT, 

3  2   1NAME  (10) ,ITMEAS,NIMEAS.MAXSAM,NTEHP 
•1  COMMON/MISC/IOlN.IOOUT.ICflRO 

5  REMIND  NTEMP 

6  REMIND  9 

7  C 

8  2  CONTINUE 

9  READ  (NTEMP. END=140)  INAME.NMEAS 

10  WRITE  (9,271)  INAME.NMEAS 

11  C 

12  4  CONTINUE 

13  READ  (NTEMP. EN0=370)  (V  (K) . K=l ,NMEAS) 
IV  WRITE  (9,271)  (V  (K)  ,K=1 , NMEAS) 

15  271  FORMAT  (20A4) 

16  GO  TO  14 

17  370  CONTINUE 

18  ENDFILE  9 

19  WRITE  (6,272) 

20  272  FORMAT  (10X.  '*«**  DATA  SAVED  ON  FILE  9') 

21  REWIND  NTEMP 

22  RETURN 

23  140  CONTINUE 

24  WRITE  (6, 141) 

25  141  FORMAT  (10X,  ■  ERROR  DATA  SET  D0E3  NOT  EXIST') 

26  RETURN 

27  END 
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